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OPERATIONAL ALIGNMENT AND 
CALIBRATION OP THE ISU FOR PHASE 2 
•OF THE* V/S TOL PROGRAM 

Volume I Calibration • 


By William H. Fincke 

DYNAMICS RESEARCH CORPORATION 
60 Concord Street 
Wilmington, Massachusetts 01887 

1. SUMMARY 

This volume of the final report presents a detailed design and 
analysis of a calibration scheme to be used for the strapdown inertial 
•system for Phase 2 of the V/STOL Program. The calibration scheme is 
capable of being implemented under hanger conditions and runs can be 
completed within one 8-hour day. The precision of calibration is expected 
to be sufficient to support flight operations. 

A total of 11 terms can be determined for each of the 3 gyros and 
8 terms for each of the 3 accelerometers. These terms provide the matrix 
of misalignments between the ISU body axes and the input axes of each 
gyro and accelerometer. Second and third order scale factor error 
coefficients are also estimated as well as bias and gyro g sensitive coef- 
ficients. A second level of calibration can also provide estimates of cross 

•j* 

compliance terms' 1 . A detailed specification of the test procedure is 

If all terms are not to be calibrated, the process can be completed within 
5 to 8 hours. Calibration of the gyro scale factor at a variety of rates 
will take longer. 
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provided, although details of the ISU alignment before each run have yet 
to be finalized. 

The expected la precision of calibration for the accelerometers 

is approximately 2 sec for misalignments, 4/ig for bias, 3 to 13 ng for 

2 

scale factor, and 6 to 12 g gig for compliance terms. For the gyros, 

the expected la precision is approximately 2 sec for misalignments, 

8 mdh for bias, 20 mdh for scale factor at 2°/sec (3 ppm), 6 to 10 mdh/g 

2 

for mass unbalance terms, 11 to 14 mdh/g for major compliance terms 
and 15 to 20 mdh/g for cross compliance terms. 

The data reduction calculations have been made as simple as 
possible since estimation of the calibration terms may be done manually 
using only a desk calculator. Flexibility is provided in the least squares 
estimation of the scale factor error coefficients in that any number of 
input conditions can be accommodated. All calculations required to 
determine the calibration terms are provided. 

A detailed error analysis of the calibration system was made that 
provided not only expected calibration estimation precision, but also formed 
the basis for specifying key performance requirements of the support 
equipment involved. Consideration was given to designing the system so 
as to minimize the complexity and number of test equipments required, 
thereby reducing the cost of implementing the system. Each major 
support equipment is described functionally. 


One of the terms for one of the accelerometers is expected to be 
30 mdh/g 2 (la). 
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2. INTRODUCTION 


2. 1 GENERAL 

The purpose of these studies is to develop operational pre-flight 
alignment and calibration procedures for the strapdown inertial navigator 
to be used during Phase 2 of the V/STOL Program. The calibration pro- 
cedure is designed to be one which can be implemented using the limited 
test equipment which will be available in the hangar. The alignment pro- 
cedure is designed to be one which can be accomplished in the helicopter 
using the fixed point, on-board flight computer. Since the calibration and 
alignment procedure developments can be discussed separately and since 
many personnel at NASA-ERC are primarily interested in only one of 
these developments, the developments are discussed and documented in 
separate volumes: Volume 1 for the calibration procedure and Volume 
2 for the alignment procedure. 

Error analyses are presented which indicate the expected accuracies 
of the developed procedures. These error analyses, of course, are de- 
pendent on the assumed input error models. Developing models for the 
inertial sensors was one of the tasks performed as part of the overall 
study. 


To the extent that ‘they were known, input parameters for the study 
were supplied by NASA-ERC. These parameters describe the following 
types of inputs: 

• specifications on the performance of the accelerometers 
and gyros. 
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6 a description of the dynamic environment in which the 
alignment and calibration must be performed, 

• a description of the test equipment that will be available 
for calibration in the hangar, and 

« a description of the coordinate frames and the mounting 
of the inertial sensors in the strapdown system. 

In those cases where the specifications and descriptions were incomplete, 
the associated error analyses must be considered preliminary in nature. 

In those cases where assumptions had to be made, the resultant error 
analysis provides information on how to set the performance specifications. 

Throughout this report, references are made to a study made by 
UNIVAC [3, 4, 5] . Whereas that study was concerned with calibration and 
alignment in the laboratory, this study is directed towards calibration in 
the field. 


2. 2 OBJECTIVES AND SCOPE OF CALIBRATION STUDY 


The objectives of the ISU calibration procedure are as follows: 


• estimate the inertial sensor calibration terms listed in 
Tables 4 _ 3 and 4~5 of Section 4, with precision goals 
as indicated in the tables 

• perform the calibration under field conditions (viz. , 
in a hangar) 


References are indicated by numbers in brackets and listed in Section 13. 
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® perform the calibration in a minimum time period, with 
8 hours as a goal. Equipment setup time and complexity 
is to be minimized, as well as pre- calibration activities 
and calculations 

® data reduction is to be as simple as possible since it 
may be done manually 


It is anticipated that the calibration procedure may be performed routinely 
once a week and more often if necessary, (e. g. , if an inertial sensor is 
changed or parameter changes are suspected). 

The scope of the study is defined by the following tasks, as 
interpreted by certain ground rules discussed in the next section: 


® develop inertial sensor mo.dels 

® develop operational procedures to obtain the desired 
calibration terms, including an estimate of the time 
required to calibrate 

• develop equations required to estimate the desired 
calibration terms 

• specify critical functional and performance requirements 
of the associated test equipment 

• identify the significant error sources and estimate the 
expected calibration precision 


The study is considered to be a part of the total development effort 
required to implement the ISU calibration function, in that further 
trade-off studies are required, as well as detailed specification of the 
precalibration activities. The study should be considered as specifying 
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a base line system that can be developed and modified as necessary to 
arrive at a final design. In this sense many engineering judgements 

were exercised with the anticipation that design details may be changed, 
as required to provide a "balanced design". 

2. 3 GROUND RULES OF CALIBRATION STUDY 

A set of ground rules were established at the beginning of the 
calibration study, to be used as a guide in developing the calibration 
system. These were discussed in Ref. 2 and further established as the 
study progressed. The ground rules of the study are as follows: 


1. The ISU is to contain three Honeywell GG 334A gyros 
and three Kearfott 2401 pendulous accelerometers, 
using time-modulated pulsed torque-to-rebalance 
loops operating at 128 and 256 kHz respectively, 

and a maximum data sampling rate of 1 kHz. 

2. Polarities of the calibration terms shall be as 
defined in Tables 4~2 and 4“4. 

3. The ISU is to be mounted in a Two- Axis Fixture (TAF) 
as shown in Fig. 5-3 (this is a change from the conven- 
tion used in the UNIVAC study (Refs. 3, 4 and 5). 

4. Accelerometer random noise is assumed to be 
negligible; gyro random noise is to be as specified 
in the UNIVAC study; torque-to-rebalance loop noise 
is no greater than the loop quantization, for both the 
accelerometers and the gyros. 

5. The test stand base motion and effects on the gyro and 
accelerometer outputs are assumed to be the same as 
that presented in the UNIVAC study. 

6. A Single-Axis Test Stand (SATS) with the rotary axis 
nominally vertical is to be used (in order to minimize 
costs). Data is to be taken only as a function of full 
table rotations. 
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7. A Two- Axis Fixture (TAF) is to be used to support the ISU. 
Accurate readouts of the TAF gimbal angles will be 
provided. 

8. To reduce costs and since data recording is to be done 
manually, the data collection equipment is to be minimal. 

A minimum number of pulse counters should be used in 
trade-off with the objective to calibrate in one 8 hour • 
day. 

9. Since data reduction may be performed manually, the 
estimation equations should be as simple as possible 
(e. g. , post run data filtering is to be avoided). 

10. Least squares filtering of redundant data is to be minimized 
and used only when necessary to provide adequate sensitivity 
or flexibility to allow, the inclusion of extra data when 
desired. 

11. The calibration procedure is to be patterned after that used 
in the UNIVAC study, modified as necessary to reflect 

a single-axis test stand (rather than two axes) and pendulous 
accelerometers (rather than vibrating string accelerometers). 
It is desirable to use the same symbology, definitions, 
zero positions, etc. 

12. The calibration procedure is to be based on using only 
test table rates and the gravity vector. Gyro scale 
factor error coefficients to be estimated using input 
rates between + 1 and + 60° / sec. 

13. Error sensitivity equations are to be identified in support 
of possible tradeoffs that may be made as the design 
develops. The estimate of expected calibration precision 
is to be based on engineering judgements of acceptable 
error tolerances in the design in order to arrive at the 
base-line system. 

14. Tradeoffs are to be minimized in this study in preference 
to the development of workable calibration procedures 
and specification of critical calibration system hardware 
requirements . 

15. Details of the pre- calibration activities are to be developed 
separate from this study. 


2-5 



3. SYMBOLS 


The symbols used in the ISU Calibration study are intended to be 
compatible to the ma ximum extent possible with those used in the UNIVAC 
Study [3, 4, 5] and in previous reports on this project [1, 2]. However, 
some deviations and redefinitions were necessary to avoid ambiguities, 
classify meanings, and provide consistency. Changes in these areas are 
noted specifically. The detailed lists of symbols are contained in the last 
Appendix, G, to facilitate rapid and convenient access whenever required. 

The information contained in the Glossary, Appendix G, is intended 
to provide a central location for understanding all nomenclature used in 
the report. Four categories are defined, as follows: 

• Formation of symbols 

• List of prime symbols and abbreviations 

• List of subscripts 

• List of superscripts 

Coordinate systems and related transformations are defined, where used, 
in Sections 5, 9 and Appendix B. 

t 

Positive quantities are defined in the direction of positive axes and 
positive angles are defined according to the right hand rule. Signs asso- 
ciated with nominally cardinal values of dot and cross products of vectors 
are defined as above and in terms of the coordinates in which the right 
hand term in the product is expressed. Positive values of the calibration 
terms are defined in Tables 4-2 and 4-4 in the Math Models Section (4). 



4. ISU AND INERTIAL SENSORS MATH MODELS 


The purpose of this section is to present the math models that 
will be used to characterize' the ISU in general and the accelerometers 
and gyros in particular. The inertial sensors are discussed in Sections 
4. 1 and 4. 2, whereas the ISU is considered below. 


The structural alignment tolerances within the ISU are assumed 
to be as listed in Table 4-1. The ISU porro prism misalignment relative 
to the ISU optical cube is not considered in this study since the ISU cube 
is considered to be the prime reference during calibration. 


Definition of Misalignment 


Value 


® Inertial sensor mounting error 6 sec (max) 

e Inertial sensor OA and SA (or PA) 10 min (max) 

about IA 

• ISU cube and inertial sensor 6 min (max) 

mounting pads 

• ISU cube and ISU mounting 6 min (max) 

surface 


Table 4 - l ISU Internal Alignment Errors 
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4. 1 ACCELEROMETER MATH MODEL 


The accelerometer math model is derived in Appendix A and 
repeated here, as follows: 

t t • 

N N N 

AV C - A$ A E C k - J> (a.) dt + .J [A b+ D a. + D a* + D a* 

1 t J t 

o o 


-C^a + C Ti;) a.a + C T _a.a - MOA a" 

Pp IP 1 p lOio p 

+ MPA a* + Q w w . + J . w + D, ] dt 
o IP p x A o b 

+ e + e n (4-1) 


where 

a' = acceleration in direction of ISU body axis to which 
^ accelerometer PA is nominally parallel 

a 1 = acceleration in direction of ISU body axis to which 
accelerometer OA is nominally parallel 

A 

and in which e includes the effects of random accelerometer noise 

A n a 

(e ) and dynamics of the accelerometer rebalance loop (e^ ). The 

term is the dynamic bias due to the nominal vibration environment. 

The torque -to -rebalance loop operates in the same fashion as the gyro, 

as explained in Refs. 1 and 9. 
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The Cp term in Eq. (4-1), which is proportional to mass unbalance 
along the accelerometer IA, is indistinguishable from MOA and so it 
will be dropped. In addition, the Q^p and terms may be deleted 
since accelerometer outputs are not recorded when the test table is 
rotating and the rates that do exist (earth rate components) are small. 

Other error effects considered negligible are: 


a 


9 

23Q.J& 2 2 2 

rotational cross coupling: 6a = — rr~ a.(w - w. - w ) 


« compliance: 6a = T (a., a , a )/mi (not including C T7 -, and C„) 

c l o p IF lO 


• cross -coupling of accelerations along OA 


where mi is the pendulosity of the accelerometer and K is the equivalent 
stiffness of the rebalance loop. 


The effect of scale factor error is modeled as a third order 
polynomial in terms of acceleration along the IA, as shown in Eq. (4-1). 

Sign conventions associated with the accelerometer calibration 
terms are contained in Table 4-2, and Table 4“3 lists numerical values 
of the calibration terms in Eq. (4-1). Table 4-3 also indicates the 
calibration terms to be estimated and the associated precision goals. 
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Calibration 

Term 

Symbol 

Sign of Output Acceleration for 
Positive Value of Calibration Term 

Bias 

A b 

positive 

Compliance 

C IP’ c io 

same as sign of products 
of accelerations along (IA) (PA) 
and (IA) (OA), respectively. 

Misalignments 

MOA, MPA 

opposite to acceleration along 
nominal PA and same as 
acceleration along nominal OA, 
respectively (corresponding to 
positive rotations of IA about OA 
and PA, respectively). 

Accelerometer 

Scale Factor < 

f 13 !' D 3 

same as acceleration along IA 

Error Coefficients 

L D 2 

positive for either positive or 
negative acceleration along IA 


Table 4-2 Sign Conventions Associated with Accelerometer 

Calibration Terms 


4-4 




Name 

Units 

Nominal 
• Valued ^ 

f ang <3) 

(max)^ 

Est. Free. 
Goal (max)^ 

AO 

Scale factor (SF) 

fps/ pulse I 

. 0025 


- 

A* 

Bias 

M g 

50 


< 25 

D i * 

SF error (SFE) 

Mg/g 

50 


•< 25 

d 2 * 

2nd order SFE 

hg/g 2 



TBD 

D ** 
o 

3rd order SFE 

Mg/ g 3 



TBD 

/-i ■j, 

(, T‘'l' 

IP 

Vibropendulous 

Coeff. 

Mg/g 2 


~ 10 

TBD 

c ** 

io 

Compliance of 
pendulous axis 

Mg/g 2 


~ 10 

TBD 

MOA* 

Misalignment of 
IA about OA 

sec 

+8 

+ 30 

< 5 

| 

MPA* 

Misalignment of 
IA about PA 

sec 

+8 

+ 30 

< 5 

. 


Notes: (1) A single asterisk refers to the basic calibration terms to be estimated. 
The double asterisk refers to additional terms to be estimated to 
provide a full calibration. 

(2) Value at beginning of calibration run, for a normally operating ISU 
(average magnitude) . 

(3) Maximum value at beginning of calibration run, to be used for worst 
case design purposes. 

(4) Ma xim um values are interpreted as being equivalent to a 3cr value. 

TBD indicates "to be determined". 


Table 4~3 Numerical Tabulation of Accelerometer Calibration Parameters 
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GYRO MATH MODELS 


The gyro math model was derived in Ref. 1 and is repeated here. 


as follows: 


2 , c k - f '’Vi dt + f [ V B lV B S V B oV G n a i 


+ C„~a +C T _a.a +C T „a.a + CL 0 a a 
SS s IS i s IO i o OS o s 


+ w.w +Q w.+ Q-W. - MOA w' 
IS i s l 2 l s 


+ MSA w 1 + J w ] dt + e G + e G 
o G o q n 


where 


w 1 = inertial rate about ISU body axis to which gyro SA 
is nominally parallel 

w 1 = inertial rate about ISU body axis to which gyro OA is 
° nominally parallel 


G G 

and in which e includes the effects of random gyro noise (e ) and . 
n q nr 

dynamics of the gyro gimbal/ rebalance loop ( e nre ] C) )* The and 

J terms may be deleted in the calibration process since they contri- 
G 

bute a negligible effect to the gyro output. 
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The effect of scale factor error is modeled as a second order 
polynomial in terms of inertial rate about the IA, as shown in Eq. (4-2). 

Sign conventions associated with the gyro calibration terms are 
contained in Table 4-4, and Table 4-5 lists numerical values of the 
calibration terms in Eq. (4-1). Table 4-5 also indicates the calibration 
terms to be estimated and the associated precision goals. The math 

Q 

model for internal random gyro drift rate (e ) is contained in Table 

nr 

4-6 and is the one used in the XJNIYAC study* as originally specified 
in Ref. 7. 


^Appendix A (pp. A- 8, 9) of Ref. 4 
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Calibration Sign of Output Rate for Positive 


Term 

Symbol 

Value of Calibration Term 

Bias 

R b 

positive (corresponding to 
a negative torque on float 
about OA) 

Mass . < 

Unbalance 

r 

MUIA 

MUSA 

B r B s’ B o 

opposite to sign of acceleration 
along SA and IA, respectively 
(corresponding to mass unbalances 
along positive LA and SA) 

same as sign of acceleration 
along IA, SA and OA, respectively 

Major 

Compliance 

c n' c ss 

positive for positive or negative 
acceleration along LA and SA, 
respectively 

Cross 

Compliance 

C IS' C IO' c os 

same as sign of products of 
accelerations along (IA)(SA), 
(IA)(OA) and (OA)(SA), 
respectively 

Misalignments 

MOA, MSA 

opposite to rate about nominal 
SA and same as rate about 
nominal OA, respectively 
(corresponding to positive 
rotations of LA about OA and 
SA, respectively) 

Scale 

f Q 1 

same as rate about LA 

Factor Error 
Coefficients 

1 s 

positive for either positive 
or negative rate about IA 


Table 4-4 Sign Conventions Associated with Gyro Calibration Terms 
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Symbol 

(1) 

Name 

Units 

Nominal 
Value (2) 

Range 
(max) (3) 

Estimation 
Precision 
Goal (max) (4) 


Scale Factor (SF) 

sec /pulse 

. 844 

(0 to 30°/sec) 
1. 688 

(30 to 60°/sec; 

+ 10% 
± 10% 


"D vO 

b 

Bias 

°/hr 

+ 0. 100 

0. 200 

< 0. 020 


SF error (SFE) 

ppm 

+ 200 


< 50 

O * 

z 

2nd order SFE 

ppm/rad/sec 

50 


TBD 

TD 4* 

B i ~ 

Mass Unbalance 
along SA 

°/hr/g 

+ 0. 300 

+ 1. 5 

< 0. 020 

D a. 

B s ” 

Mass Unbalance 
along LA 

°/hr/g 

+ 0. 300 

+ 1. 5 

< 0. 020 

"D 4,4, 

■a- 

Output axis accel- 
eration sensitivity 

°! hr/g 

. 040 


< 0. 020 

O 

H 

* 

Compliance Coeff. 

°/hr/g 2 

. 020 


< 0. 025 

C 

ss 

Compliance Coeff. 

°/hr/g 2 

. 020 


< 0. 025 

MOA* 

Misalignment of 
IA about OA 

sec 

± 8 

± 30 

< 5 

MSA - 

Misalignment of 
I A about SA 

sec 

+ 8 

+ 30 

< 5 

c ** 

IS 

Compliance Coeff. 

°/hr/g 2 



< 0. 025 

r* ** 
IO 

Compliance Coeff. 

°/hr/g 2 


• 

< 0. 025 

(, *K- 1 ' 

OS 

Compliance Coeff. 

°/hr/g 2 



< 0. 025 


Notes: (1) A single asterisk refers to the basic calibration terms to be estimated. 
The double asterisk refers to additional terms to be estimated to 
provide a full calibration. 

(2) Value at beginning of calibration run, for a normally operating ISU 
(average magnitude). 

(3) Maximum value at beginning of calibration run, to be used for worst 
case design purposes. 

(4) Maximum values are interpreted as being equivalent to a 3cr value. 
TBD indicates "to be determined". 

Table 4-5 Numerical Tabulation of Gyro Calibration Parameters 
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* G 

The gyro random drift rate (e^J is modeled in terms of power 
spectral density, as specified in Ref. 7 and Appendix A of Ref. 4. 

f(Hz) PSD (e G )*, <deg/hr) 2 /Hz 

nr & 


0 to (100 hrs) _1 

3. 02 x 10 2 

(100 hrs) 1 to 10 3 

49-9 

2. 33 x 10 °f 


-7 

10 to 8 

2. 33 x 10 
-6 

8 to 24 

1x10 

-4 -2 

above 24 

6x10 f 


• G - 1 

The rms value of PSD (e ) from f = (14x3600 sec/cycle) to 00 is 

5 mdh. 


Table 4 _ 6 Gyro Internal Random Drift Rate Math Model 


A plot of the PSD has a large discontinuity at f = (100x3600 sec/cycle) 

• Q 

however, the effect of e in the calibration process is small. 

nr 
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5. ISU AND TEST EQUIPMENT CONFIGURATIONS 


The configurations of the Inertial Sensing Unit (ISU) and required 
calibration test equipments are discussed in this section. In order to 
support subsequent analyses and understand the system operation, 
considerations are restricted here to geometrical relationships within 
the major equipments and definition of all necessary coordinate frames. 
Hardware considerations and functional descriptions are discussed in the 
next section (6). 

The layout of the ISU is defined in [6] and repeated in Fig. 5-1. The 
change from the ARMA D4E vibrating string accelerometer to the Kearfott 
2401 pendulum type of accelerometer is reflected in the figure, primarily 
as a reversal of the LAs and definition of OAs and PAs. The addition of 
a porro prism azimuth reference is also shown. In this study, the 
or M faces of the ISU cube will be used as the azimuth reference and it 

O 

is assumed that the alignment between the cube and the porro prism is 
determined separately. Note that the normal mounting position of the ISU 
in the aircraft is with the black cover down such that +z is down, x is 
forward and y is right. 

The ISU cube and inertial sensor relative orientations indicated in 
Fig. 5-1 are repeated schematically in Fig. 5-2 for clarity, and the 
entire ISU has been rotated to show its orientation relative to earth coordi- 
nates for the zero positions of the Single-Axis Test Stand (SATS) and the 
Two-Axis Fixture (TAF). Detailed definitions of the coordinate frames 
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Notes: (1) When installed in aircraft, Z is down, X is forward and Y is right. 

(2) The Body Axes are defined by the ISU Cube faces, nominally 

in the directions of the gyro IAs. 

(3) The Mirror Axes (jVL MM) are defined by the visible surfaces of the 

ISU Cube. J 


Figure 5-1 Layout of Strapdown Inertial Sensing Unit (ISU) 
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N 



Gyro 

Unit 

Cube 


f 


) 

v 


IA = OA x SA 
G = 0 G iS G 


~2 % 



Notes: (1) The ISU orientation relative to earth coordinates is shown for zero 
positions of SATS and TAF. 

(2) All coordinate sets are defined by unit vectors. 

(3) The Body Axes are defined by the orthogonal ISU cube faces, 

nominally in the directions of the gyro IAs (i. e. , 33^ = M-] , Bg s M 2 
and Bg = -Mg). All other coordinate sets are generally misaligned 
with respect to each other. 

(4) The nominal orientation of the Bubble Level coordinate frame 
(L^L^Lg), which is fixed to the TAF inner gimbal, is shown for 
reference purposes. 

Figure 5-2 ISU Cube and Inertial Sensor Orientations for ISU 
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shown in Fig. 5-2 are contained in Table 5-1. The positive ISU cube 
(mirror)axes were chosen to correspond to the visible faces of the cube. 
Except for M„, the Mirror and Body axes coordinate frames are the 

u 

same as in the Univac report [3], as are all of the gyro coordinates. The 
accelerometer system (IA) coordinates are also the same, except for A . 
Finally, the zero position of the ISU has been changed from that used in 
the Univac study. This was done primarily to facilitate mounting some 
electronics, along with the ISU, within the TAF. 

The configuration of the major alignment equipments associated 
with the ISU is shown in Fig. 5~3, and Table 5-2 contains detailed defini- 
tions of the coordinate frames shown and how they are related. Unit 
vectors are used to define all coordinate axes. Unlike the Univac study, 
a test table trunnion axis is not provided, in order to reduce costs. The 
SATS and TAF are used to support the ISU in various orientations relative 
to the rotary axis and the Earth's coordinate frame. The autocollimator 
is used to measure the orientation of the ISU cube relative to vertical, 
about the S axis, and a set of at least three Bubble Levels (BL) are 

O 

used to represent the inner gimbal axes of the TAF. The ways that these 
equipments are used is discussed in Sections 6 and 8. The ISU cube is 
shown in Fig. 5 - 3 for reference purposes to show the nominal relationships 
between the various cube and TAF coordinate frames. 
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Figure 5 3 ISU Test Equipment Coordinate Frames 



Notes: (1) The SATS and TAF are shown in their M zero n positions (i. e. , $ 2 ~ ^ 3 = ^ ' 

(2) The F and Earth's coordinate frames are assumed to be misaligned 
by 60 ^, about ISL 

(3) The Bubble Level coordinate frame is nominally defined by the inner 
gimbal coordinate frame. 

(4) Nonorthogonalities of axes are indicated by defining the O and R coordinate 
frames as small rotations of the O' and R' frames about the O' and R' axes, 
respectively. 
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M is defined by the visible surfaces 
of the ISU cube and is assumed to 
be orthogonal and right-handed in 
the order listed. 

B is defined by the ISU cube faces, 
nominally in the direction of the 
gyro I As. 

The elements of G correspond to 
the gyro IAs and are generally 
not orthogonal. 

The gyro right-handed coordinate 
system is defined relative to the gyro 
case and is generally not orthogonal. 

The elements of A correspond to 
the accelerometer IAs and are 
generally not orthogonal. 

The accelerometer right-handed 
coordinate system is defined 
relative to the accelerometer case 
and is generally not orthogonal. 


E is defined by j? and w 
as an orthogonal right- 
) xU > handed coordinate system. 


and Earth Coordinate Frames 



All coordinate frames are defined by. a triad of unit vectors, unless specified 
otherwise. 


Optical Alignment Coordinate Frame 


• Auto collimator 
(S) 


_S £ line of sight 


S is defined by_g_ and the optical line 
of sight. 


Single-Axis Test Stand (SATS) 

• SATS Frame 
(F) 

F F F 
-1-2-3 

• Test Table 
(T) 


* Rotary Frame' 
(R 1 ) 



F is fixed to the frame of the SATS 

T defines the rotary axis, but does not 
rotate. 

R contains the rotary axis and rotates 
the test table 


Two-Axis Fixture (TAF) 


o Rotary frame R R R 
(R) 1 1 6 


0 Outer Gimbal' O' O' O' 
(O') 1 ^ 6 


R defines the outer gimbal axis and is 
misaligned with R' by MR3', about R^. 
R = R' for MR3' =0 and does not rotate 
with the outer gimbal 

O' contains the outer gimbal axis and 
rotates with the outer gimbal 


» Outer Gimbal P-.O 9 .Q 0 
( 0 ) 1 <4 o 


• Inner Gimbal 

(I) 




O defines the inner gimbal axis and is 
misaligned with O' by M03', about O' . 

0 = 0' for M03' = zero and does not 
rotate with the inner, gimbal. 

1 contains the inner gimbal axis and rotates 
with the inner gimbal 


Bubble Level L.L„L„ 

(L) -1-2-3 

-Level indication marks 



Bubble 

Level 


A unit vector in the L frame is defined by 
the line through the ends of the level indi- 
cation marks that points in the same 
nominal direction as the corresponding 
inner gimbal unit vector (see diagram). 
Elements of L are generally not orthogonal. 


Table 5-2 Definitions of Test Equipment Coordinate Frames. 
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6. DESCRIPTION OF SYSTEM AND 
TEST EQUIPMENT REQUIREMENTS 

The purpose of this section is to describe the total calibration 
system from the equipment viewpoint, considering both functional and 
performance requirements. 

Not only will this help to understand the calibration process, but 
it is the means by which the key test equipment requirements will be 
established. After a brief description of the system, the following major 
test equipments are considered separately: 

® Single -Axis Test Stand (SATS) 

• Two -Axis Fixture (TAF) 

• Optical Alignment Equipment 

® Data Collection System 

• Data Processing System 

In conjunction with specifying the calibration procedure and the data pro- 
cessing equations, the error analysis of Section 10 and the test equipment 
requirements of this section were developed simultaneously in an attempt 
to achieve a balanced design. As mentioned earlier, however, the design 
described here is considered to be an initial or interim one, particularly 
in the area of the pre- calibration alignment equipments. In addition, the 
data collection system was not considered in detail in this study since no 
significant problems were anticipated. Instead, emphasis was placed 
upon determining procedures, estimation algorithms, expected accuracy. 
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and requirements of the critical test equipments. The performance 
requirements specified herein therefore should be considered preliminary. 


6.1 DESCRIPTION OF SYSTEM 


A brief functional description of the major system equipments and 
their interrelationships is contained in this section. There are three 
primary functions performed by the equipments provided within the cali- 
bration system, as indicated in the functional schematic of Fig. 6-1. The 
first is the ISU support and alignment function. The second is the data 
collection (and recording) function, and the third is the function of processing 
the data obtained to determine the required calibration terms. Each of 
these functions is described next. 

The ISU support and alignment function is provided by the Single- 
Axis Test Stand (SATS), Two-Axis Fixture (TAF) and an Autocollimator. 

The ISU to be calibrated is mounted within the TAF, which in turn is 
supported by- the SATS. The geometry of the assembly is illustrated 
schematically in Fig. 5-3, and Fig. 5-1 is an illustration of the ISU. The 
autocollimator is referenced to the vertical (indicated by g in Fig. 6-1) 
and with the bubble levels on the inner girnbal of the TAF nulled at three 
different orientations of the ISU, the transformation matrix between the 
ISU cube (mirror coordinates) and bubble levels (T'^"') is measured and 
recorded. The various orientations of the ISU are achieved by adjusting 
the TAF inner and outer gimbals about their respective axes and by rotating 
the entire TAF about the SATS rotary axis, which is nominally vertical. 
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Figure 6-1 Functional Schematic of Calibration System 
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Data Processing 
( ML System 


Pre- 

Calib Inertial 
Align Sensor 
Info. Info. 



Environ, 
and Test 
Equipment 
Info. 


*4*3*2 w ' 


i0 2 °r4t|J 

Preset j 


.Record data for each 
test condition. (m) 


Notes: (1) Not shown are power supplies, ISU electrical support equipment or all equipments required 
to install the Single -Axis Test Stand (SATS) and initially align and calibrate the Two- Axis 
1 Fixture (TAF-), 














Continuous precise readouts of the TAF inner and outer gimbal angles 
(0 and 0 ) are provided. The test stand is also capable of driving the 

Q O ij-i 

TAF and ISU at a rate (w ) about the rotary axis and full revolutions of 
the test table are indicated by a pulse output. Before each run, the ISU 
is nominally oriented in the desired position and the Bubble Levels nulled 
by adjusting 0 and/or 0 . The gimbal angles obtained are labeled (0 ) 
and (0^) and are recorded. Next, the gimbals are readjusted to indicate 
(0g) and {0^ , which are equal to the nominal values desired plus the 
small corrections determined when the test equipment was originally zeroed. 

The next major part of the calibration system is the Data Collection 
System. In order tominimize the cost of test equipment, only four pulse 
counters are provided to count the number of pulses generated within the 
ISU accelerometer and gyro rebalance loops. It is therefore necessary 
to select the inertial sensor outputs, as shown in Fig. 6-1. A timing 
reference is used to control the start and stop times over which pulses are 
counted, and it indicates the corresponding elapsed time. It is synchronized 
with the Ik Hz reference in the ISU such that counts are started and ended 
only for full cycles and when the lk Hz reference is zero. Initiation of 
the counts is controlled by one of two independent variables which must be 
pre-selected. When the table is non- rotating, time is selected as the 
independent variable and when the start command is given, pulse counting 
is initiated as soon as the next full lk Hz cycle is started, as described 
above. Counting continues until the elapsed time is equal to a value preset 
into the time reference. The actual stop time is determined as explained 
above. When the table is rotating, test table angle (0^) is selected as the 
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independent variable. The SATS is designed to provide a pulse for each 
full revolution of the test table. The first table pulse to be received after 
the start command is given is used to initiate pulse counting as soon as 
the next full lk Hz cycle is started, as described above. Counting con- 
tinues until the number of table pulses received is equal to a value preset 
into the timing reference. The actual stop time corresponds to the end 
of the full 1kHz cycle in progress at the time the table pulse is received. 

As stated in the ground rules initially established for this study, 

data recording and processing is to be done manually. Therefore, pulse 

G A 

counts from the gyros and accelerometers (P and P ) and the timing 
Atm $m 

reference (P and P ) are read from the counters and recorded 
after each run (m) for the particular inertial sensors (j) selected. Data 
processing is performed manually using a desk calculator. However, this 
ground rule is currently under review and more automatic data processing 
{and possibly data collection) may be adopted. Use of the ISU flight 
computer is being considered since it would provide many desirable 
features. 

6. 2 SINGLE-AXIS TEST STAND (SATS) 

The SATS must be capable of supporting the weight of the ISU and 

TAF, and yet retain definition of its rotary axis within several arc 

T 

seconds (lcr). This is only a requirement for w = + 2°/ sec and when 
T 

w =0, since the alignment between the bubble level coordinate frame 
and the rotary axis must be measured. Leveling screws are required to 
adjust the rotary axis with respect to the local vertical with a sensitivity 
of approximately . 2 min. 
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The table rate drive must be adjustable and accurate to . 5% and 
respond to changes in rate and stabilize within about . 5 minute. The 
nonuniformity of the table rate (NUWT) must be. such that 

360°/w T 

NUWT = 15 J sin(w t+0)dt (for w in deg/ sec) 

° ( 6 - 1 ) 

— T 

is less than 4. 3 sec (lcr) (per revolution of the table) for w = + 2°/ sec 

and any given value of (ft. Assuming a maximum positive rate error over 

half a revolution and a maximum negative rate error over the other half, 

T 

the percentage of maximum rate error to w = 2°/sec is .25%. These 
values are to be considered maximum permissible values and it would be 
very desirable to reduce them by a factor of 4 or more (see the error 
analysis in Section 10. 5). 

The test table pulse output that indicates full table rotations should 
be repeatable within 6 sec (la) as a maximum permissible value. A 1 to 
2 sec error would be desirable, particularly for the higher table rates for 
which a sufficient number of table revolutions can be achieved in a shorter 
time, thereby reducing the test time. 

The test stand should also be provided with a means to repeatably 

position the test table within approximately . 5 min about the rotary axis, 

at a zero position determined during initial installation. It would also be 

useful to provide table angle indicators every 180° or better yet 90° in 

order to speed up the process of setting the table rate, particularly for 
T 

w = + 1, 2 and 4° /sec. To achieve . 5%, the angle indications would need 
be accurate only to approximately . 5%. 
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Slip rings must be provided to provide power to the ISU and signal 
outputs to the support equipment. 

6. 3 TWO-AXIS FIXTURE (TAF) 

The TAF is required to provide to the operator precise readouts 
of the inner and outer gimbal angles (0 and 0 , respectively) to an 

O 

accuracy of several arc seconds (Ict). Each readout is only required 
within approximately + 2 min of the cardinal and intercardinal values 
and may be of any convenient form (electrical, mechanical scales, optical, 
etc. ), consistent with rapid readout. In addition, the gimbals must be 
adjustable and capable of being damped quickly and secured within the 
accuracy requirement. 

The inner gimbal shall permanently support 3 to 4 bubble levels 
mounted in an orthogonal fashion so as to represent the inner gimbal 
coordinate frame (see Fig. 5-3 and Table 5-2) within an accuracy of 
approximately 1 min . The directions of the tops of the vials are to be 
defined, as well as the need for the fourth bubble level. The repeatability 
of nulling the bubbles shall be no worse than 1 sec (lcr). Whenever any 
bubble level is in a position that it could be nulled, there shall be a clear 
view of the bubble. 

The inner gimbal to outer gimbal nonorthogonality should be less 
than several arc seconds. However, it is possible to accomodate much 
larger misalignments (up to approximately 20 sec.) and compensate the 
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measured data for the resulting effects. A similar condition applies to 
the outer gimbal to table rotary axis nonorthogonality. 

The gimbal axes must remain defined and stable within the above 
accuracy requirements, which implies that the entire gimbal structure must 
remain stable. Also implied by this is the requirement that mounting and 
bolting down an ISU in the TAF does not distort the bubble level coordinate 
frame relative to the inner gimbal and rotary axis coordinate frames. ISU 
mounting repeatability shall be better than 6 min (max) and the misalign- 
ment between the ISU mounting surface and inner gimbal coordinate frame 
should be less than 6 min (max). 

The TAF structure shall not obstruct a nominally horizontal line 
of sight betweenan externally mounted auto collimator and the ISU cube 
normals that are nominally horizontal, for any of the cardinal or inter- 
cardinal values of 0 and/or 0 . Means shall be provided for aligning 

O tr 

the outer gimbal axis optically with true East, within an accuracy of 
approximately 20 sec (Ict), by adjusting the test table angle (0 ). 

Ci 

6. 4 OPTICAL ALIGNMENT EQUIPMENT 

The autocollimator shall be capable of measuring the angle in the 
nominally vertical plane between horizontal and a normal to the ISU cube 
face that is nominally horizontal, over a range of + 20 min. The align- 
ment accuracy to horizontal shall be within 1 sec (Ict) and the instrument 
readout accuracy shall be within X sec (Ict). Operability shall be sufficient 
to allow accurate measurements to be made within approximately 6 
minutes. 
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DATA COLLECTION SYSTEM 


6. 5 


The primary performance requirement of the Data Collection 

System is to measure the elapsed time (At) of each run within an accuracy 

of 100 p sec (lcr), and preferably 10 to 20 p sec, in order to reduce the 

calibration time at the higher rates. Furthermore, pulses are to be 

counted at the maximum rate of 128 k Hz for the gyros and 256 k Hz for the 

T 

accelerometers. The timing requirements when the table rate w is to be 
determined are somewhat more involved due to the delay between the time 
a table pulse is received and a full 1 k Hz cycle starts. This delay can be 
as large as 1 millisec, which is unacceptable for estimation of the scale 
factor error terms. A possible solution to the problem is to measure and 
indicate the delays both at the start and stop of pulse counting and compute 

the average table rate separately using a different At than for the gyro 
pulse count. 


For the cases when the inertial sensor IAs are up (or down), a 
large number of pulses will be generated. For w = 60°/ sec, the total 
gyro count after 6 minutes of running will be approximately 46x10 . 

Therefore, an 8 digit counter is required. For the accelerometer, the 

* fi 

total count after 10 minutes will be approximately 15, 4 x 10 and a 

8 digit counter should be adequate. When the IAs are both at 45° to the 

horizontal, the maximum counts will require two, 7 digit counters to be 

used. The fourth counter may have a greatly reduced range since the 

input rates and accelerations are considerably smaller. 


■v 


as specified in the operational procedure (Section 8) 


6-9 



7. DEVELOPMENT OF CALIBRATION EQUATIONS 


The purpose of this section is to apply the accelerometer and gyro 
math models presented above (in Section 4) to the ISU assembly, as supported 
by the Two-Axis Fixture (TAF) and the Single-Axis Test Stand (SATS). The 
linear accelerations and angular motions imposed upon each of the inertial 
sensors are related to the gravitational attraction of the earth and the 
rotation of the SATS and the earth in terms of the SATS and TAF relative 
orientations. Various positions (m) of the ISU are chosen as discussed in 
the next section (8), and the resulting equations are then simplified by 
deleting terms that contribute negligible effects. Both the gyro and accel- 
erometer equations are arranged in the following form: 

9 

= H m . Y, + 6 m“ (7-1) 

3 3 3 3 3 

where superscript m refers to a particular ISU test condition (orientation), 
j to a particular inertial sensor, and 

= adjusted measurement inertial sensor 
i for run m. 

At m = t - t s time over which measurement M is 
N o 

obtained, for run m 

• xxi 

]$L = average rate of M over time At 

Y. = column vector of calibration terms for 
“3 

instrument j (defined below). 

E* 1 = row vector of coefficients that relate Y. to M 1 ? 1 
~ ^ for a given run m and instrument j ^ ** 

SM m = error in measurement M 111 
2 2 
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The measurement 3VL is the indicated change in the output of inertial 
sensor j over the time At, corrected(or adjusted)for known and significant • 
inputs (such as earth rate and/or table rate and/or gravity). Errors in 
the measurement and correction process, as well as neglected terms, 
constitute SM 111 . Therefore, since H* 1 . Y . represents the effect of the 

J -J-J 

terms to be calibrated, Eq. (7-1) can be used as a basis for estimating 

the calibration terms, as described below, by neglecting 6 m. . 

J 

In essence, id*. Y . is a linear combination of the various calibration 
1 ~J 

terms to be estimated. Different combinations are obtained by changing 

the test conditions (m) - (viz, table rate and ISU orientation). When the 

number of linearly independent combinations (or measurements M) is 

equal to or greater than the number of calibration terms (unknowns), 

sufficient information is available to solve the set of equations , thereby 

making estimates of the various calibration terms. The .accuracy with 

which this can be done is a function of the 6M* 11 terms and the coefficients 

3 

that constitute the solution to the equations. Section 10 contains a detailed 
error analysis using Eq. (7-1) as a starting point. 

In the case of the scale factor coefficients, Eq. (7-1) is used only 
to estimate the effect of various scale factor error terms. A separate 
regression equation is used to estimate the particular terms in order to 
simplify the estimation processes, while simultaneously providing flexi- 
bility in the magnitude and number of instrument input rates and accelerations 

Assuming no a priori, statistical information concerning the calibration 
terms and/or error sources is available. 
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that can be included in the regression analysis. This is discussed further 
in Section 7. 3 below. 


7. 1 ACCELEROMETER EQUATIONS 


In the calibration test procedure presented in Section 8 it will be 
noted that the accelerometers are always oriented in one of these general 
positions (viz . , IA vertical, IA horizontal (with either OA or PA nom- 
inally vertical or at 45° to the vertical) or IA at 45° to the vertical (with 
either OA or PA nominally horizontal). The effect of the calibration 
terms on the outputs of the accelerometers for each of these three 
positions are determined from Eq. (4-1) as follows: 


For IA Vertical 


N 

'V“ = J <aPdt + At m {6 V +1 , 6V 1 } m +e Am 

J 3 t C C 3 




where 


t__ 

J* (a.) m dt= integral of specific forces sensed by 
t 1 J 

o accelerometer j for test condition m 

A j / A tt \HL 

S At (A/ U) 
g —J ~n 


J 


= pulse count out of accelerometer j for test 

condition m 
N 


(7-2) 


(7-3) 


(7-4) 
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Am _ 
e . = 

3 


all accelerometer j errors in position m, other 
than those shown 


’ +1 ' -i 

and Sv and 6V are the average acceleration errors due to scale 
c c 

factor errors when the IA is up and -down, respectively. It will be noted 

' +1 ■ -i 

that the accelerometer bias (A^) is included in 6v^ and 6V c , for two 
reasons. First, it was learned in the development of the estimation 
equations (in Section 9) that more accurate estimates of the other calibra- 
tion terms would result by doing this. Second, the terms when 

estimated in this form can be used directly in the regression analysis 
described below (Section 7. 3) for estimating the various scale factor 
error terms. Only the nominal value of the dot product term will be used, 

as identified by the subscript n. Errors due to this assumption are defined 
-Ani 

as part of e_. and are considered in the error analysis. 


For IA horizontal 


,.A n Am 

P. 

3 3 



Am 
+ e. 

3 


(7-5) 


where P (a.), dt = 0 
t 

o 

A! 3 accelerometer bias as modeled in accelerometer 
b 

calibration equation 

and NOA and NPA are the nonorthogonalities with respect to Vertical 
of the IA about OA and PA respectively. Note that 
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(PV jjf)^ is nominally equal to g when the PA is nominally vertical, and 
to gl when PA is nominally 45° from the vertical. Also the product 


NOA(P.* _g) is of second order when PA is nominally horizontal, and in 
this case it is neglected. Similar reasoning applies to the term 

NPAtO^f 1 . 


For IA at 45° to the Vertical 


. A-rtAm 


ASfp""* = J* (a.)“dt + At m [{ -NOA(P.- gf a + NPA(0 A - £) m } 
3 3 t 1 3 3 n 3 n 




+ {6V + ‘ ? , 6V~' 7 }] + e 7 ^ 11 
c c j j 


where f (a.) m dt = At m (A‘g) in 
13 —i n 


•+,7 ' 7 

and 6v ‘ and 6V " are the effects of scale factor error when the IA 
c c 

is 45° above and below horizontal, respectively (in which case 


lAj’Sl n = • 707 g)- As in the case for the IA vertical, A fa is included in 
" +. 7 ' - 7 

6V * and 6V * . Since Eq. (7-6) applies only for the IA at 45° to 
C . C A 

vertical and either Por O horizontal, one of the two complicance terms 


(Cjp, Cj.q) nan always be neglected. Again, nominal values of the dot 
product terms will be used, as identified by the subscript n. 
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The above accelerometer calibration equations (7-2, 7-5 and 7-6) 
may be combined into a general form, consistent with that of Eq. (7-1), 
as follows: 



“3 ~J 3 


(7-1) 


where 


M 111 = A<i> A (P A /At) m - 
3 3 3 


g(A'U) 

J 


m 

n 


(7-7) 


Y. - [A^ j (-NOA)g | (NPAJg | 6V* 1 | Sv" 1 ; 6 V*' 7 j 6 V/ 7 j Cj p g 2 | C^g 2 ] J 


H. [( 0 )*! (P..u) i(0 A U) ! (D* ' (-1)* ' (+.707)=!= j (-.707)* ; 

“ 3 i 3 n i 3 n 1 ' ' ' ' 


(A.-U) (P.-U) i(A.-U) (0 A U) ] m 
“ 3 - n - n t —3 — n — 3 - n j 


(7-8) 


(7-9) 


and the asterisk (#) indicates those elements of the row matrix are plus unity 

only for the cases when (A.*U) = the value shown in parentheses. In 

""3 — n 

all other cases the elements are zero. As before, nominal values of the 
dot product terms (as identified by the subscript n) are used and the 

• 

second order effects due to this assumption are defined as part of SM* 11 . 

* m ^ 

The term 6 M. is considered in the accelerometer error analysis 
J 

(Section 10. 4). 


Appendix B includes a complete tabulation of the basic calibration 
equations for each of the accelerometers in each of the orientations 
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specified in the calibration test procedure (Section 8). These equations 
are used as a basis for deriving the accelerometer estimation equations 
in Section 9. 


7. 2 GYRO EQUATIONS 


As is the case with the accelerometers, the gyros are always 

oriented in one of three general positions. These are described in 

Section 8 and correspond to the gyro IA parallel to the test table rotary 

axis, IA normal to the table axis, or IA at 45° to the vertical (with 

either OA or SA nominally horizontal). As described earlier, the test 

table axis is nominally vertical and additional inputs are introduced to 

T 

the gyros during calibration runs by rotating the table at a rate w . 

The effect of the calibration terms on the outputs of the gyros for each of 
the three general positions are determined from Eq. (4-2 ), as follows: 


For I A parallel to table axis (and nominally vertical) 

t. 


A* G P 2111 = r (W.) m dt + At m [R' + BjG.-gfV 
3 3 J t ij b I-j-^n 


N 


2 , r ' + c ‘ m, , Gm 
C n g + ] j +e j 


(7-10) 
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where 


N _ 
r (w.). dt 

i 1 j 


p Gm 

j 


integral of significant rate inputs to gyro j, with respect 
to inertial space,for test condition m 

A0^(G.*RI) m + At m w sinL (G.-U-) m (7-11) 

2 — 3 -1 n e — 3 — n 

pulse count out of gyro j for test condition m 
N 

£ C. (7-12) 

, k 


A0 


m 


= change in table angle 0 over time At 


m 


e Gm 

i 



- all gyro j errors for test condition m, other than those 
shown 

= gyro bias, as modelled in gyro calibration equations 


and 60 ' and 60 are the gyro rate errors due to scale factor errors 

^ C rj-i 

when the IA is parallel to the table axis and when the table rate (w ) is 

plus and minus, respectively. The effect of scale factor error when 
T * 

w = 0 is assumed small and is considered as an error source. Because 

of the way these terms are measured and processed, the term Rj^ is 

considered separately (unlike the case with the accelerometer in which 

* +1 * + 7 

A, is included as part of 6V , 6V ’ , etc. ) A detailed discussion of 
b c c c 

this is included in the calibration test procedure development considera- 
tions (Section 8. 1). As in the case of the accelerometers, only nominal 


It is assumed that gyros available for field use have been previously 
calibrated to a point that when wT = 0, the effect of earth rate inputs 
times the scale factor error at such relatively small inputs is negligible. 
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values of the dot product .terms will be used. Errors due to this assumption 
are' defined as part of e. and are considered in the error analysis. 


For IA normal to table axis (and nominally horizontal) 
t. 


N 


A$ G p Gm = J (w.) m dt + &t m [RJ + B (S G _g) m + fWO G _g) m 
2 3 1 2 b n -O'— j *'n 

o 


+ C co (sP* g) 111 exp2 -NOA w m + NSA w m ], + eP* 11 
SS —j ^ n * s o j 3 


where 


U N 


l N 


J (w.) m dt = w J [(G.-N) m cosL -l- (G.- R;) m sinL]dt 
1 .1 ® x J .1 -*• 


At w (G.*R1) sinL 
e -3 —1 


At m w [(G.*R J Xii sinL + (G.- N) iiA cosL] for w -1- = 0 


,m 


m 


T 

for w / 0 


T 


e —3 —1 


“J 


t 

W S 3 


T, G _ n m 
w (S. • R') 
-3 -1 n 


.m 

(w o ) J 


w (O. • RJ 
-3 -In 


At 


m 


h T - t 
N o 


(7-13) 

(7-14) 

(7-15) 

(7-16) 

(7-17) 

(7-18) 

(7-19) 


and NOA and NSA are the nonorthogonalities with respect to the table 

T 

axis (RI) of the IA about OA and SA, respectively. For w f 0, the 
integral of the first term in Eq. (7-14) will be essentially zero, for 
any initial orientation of the IA relative to North, since the gyro output 
is utilized only for full rotations of the test table (as discussed in Section 
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8. 1). The integral of the second term can be significant and so it cannot 

s 1 1 > 

be defined as part of the error e . . For w =0, both terms can be 

J 

significant and so actual values of these dot products must be used instead 
of nominal values. Only nominal values of the other dot product terms, as 
identified by the subscript n, will be used. It will also be noted that sub- 
stituting Eqs. (7-17 and 7-18) into Eq. ("7 — 1 3 ) results in coefficients of 

NOA and NSA equal to -A0^ 1 (S^* R') m and A0** 1 (O^' R’)™, respectively. 

2 — 2 — 1 n 2 — 2 — 1 n 


For IA at 45° to the Vertical (and w =0) 


N 


. - G -^Gm f* . .m,, ,,ui. . _ xxi , _ m _ m _ , in. 

A#. P. = (w.) . dt + At [R * + B T a. + B„ a + EA a + C TT (a. ) exp2 

.11 " i] b I i S s O o II i 


where 


, „ m rn m m , m m, Gm 

+ C co (a )exp2 + C TC ,a. a +C T ~a. a + C^ e a a ] + e. 

SS s IS i s IO i o OS o s j 


f (w.J^dt = At m w [(G.*U) m sinL + (G.* N) m cosL] 

U T “1 p " ,_ l *” 1 ■"'I ““ 


1 1 


e — j 

m TT .m 

a i 

=g<s?-u) m 


~1 


m 


2 — n 


G T , v m 


C = g(o , *U)' 


o 


2 ~ n 


(7-20) 

(7-21) 

(7-22) 

(7-23) 

(7-24) 


As in the case for the IA normal to the table axis, the actual values of the 
dot products in Eq. (7-19) must be used’, whereas only nominal values will 
be used for the other dot product terms, as identified by the subscript n. 
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The above gyro calibration equations (7-10, 7-13 and 7-18) may 
be combined into a general form, consistent with that of Eq. (7-1), as 
follows: 


m“ = H 0 ”- Y. + 6 M m 


(■7-1) 


where 


• m . 

M . = J 
2 N 


[^fp G - A 0 ^(G.‘R') m ]/At m - w (G.-RL)fsinL 
2 2 z ~2 ~ 1 n e —3 —1 - 


A$ G (P./At) m - w [(G.*R»)“sin L + (G.* N) m C os L] 
3 J e — 3 — 1 * 


for w T f 0 


(7-25) 


for w T = 0 


(7-26) 


* = the subscript n for IA nominally vertical. For IA level or at 45°, 
the actual value of the product is used. 

*3 - [ s b ! .B, i g B s j g B 0 : g 2 c n | g 2 c ss ; g 2 c IS ; g 2 c IO ; g 2 c os i hVoA ; 


Ti 1 * + 1 ; - t 
w |nsa;60 c ; 6 #> c r 

J 


(7-27) 


r Gm 


,G 


1 ,„G 


,G 


H = [1 i (G.* U) i(S-U) 1(0. -U) i(G.-U) exp2 i(S • U) ex P 2 i 

-j . -j - m 1 -j -ni“] - n | - n , -j — n 1 


! <2,- S) n (sf ■ U) n : (a,- £> n <o G -E> n i ( 0 °. U) n (S°. 2 > n i 

-K (S?*R!) ;K (0 G -RI) !(1)** ! (-l)**] m 
x 1 x n x i x n 1 * 


(7-28) 


--indicates elements = unity only for case when (G. R’ ) = the value shown 

] 1 & 

in parentheses. Otherwise, elements = zero. 
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(7-29) 


0 for w T = 0 
Kj = \ +1 for w T = plus 

-1 for = minus 

The second order effects due to using nominal values of the dot product 

terms, identified by the subscript n, are defined as part of CM 01 . The 

* m ^ 

term 6IVL is considered in the gyro error analysis (Section 10. 5). 

Appendix B includes a complete tabulation of the basic calibration 
equations for each of the gyros in each of the orientations specified in the 
calibration test procedure (Section 8). These equations are used as a 
basis for deriving the gyro estimation equations in Section 9. 

7. 3 SCALE FACTOR ERROR COEFFICIENT EQUATIONS 


In Section 4, the errors of the accelerometer and gyro outputs 
that are functions of the respective inputs were modeled as polynomials 
of the following form: 



A, + D,a. + D*a 2 + D o a 3 
b 1 i 2 l 3 l 


6 0 = R, + Q t w. + Q„w. 

c b 1 i 2 l 


(for each accelerometer j) 

(7-30) 

(for each gyro j) (7-31) 


where 

6v s error in indicated acceleration from accelerometer, 
c 

due to scale factor errors 

m 

60 ^ = error in indicated rate from gyro due to scale factor errors 
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= acceleration input to accelerometer 
w^ = rate input to gyro 

and and R^ are the accelerometer and gyro biases defined as constants 
in the scale factor error equations.' The distinction between and 
and and R^ is discussed later in this section. The A, R, D and Q 
coefficients are to be estimated for each inertial sensor j as part of the 
ISU calibration process. Although the coefficients could be included as 
part of the accelerometer and gyro calibration equations derived above, 
considerable flexibility and some improvement in accuracy is realized 
by using Eqs. (7~30 and 7 -31). This will be demonstrated as the theory 
is developed below. 

Assuming no a priori knowledge re: the expected values of the 
coefficients, at least as many pairs of input/output measurements as 
terms to be estimated are required to realize a solution. This is done 
by defining the following regression equations, based on Eqs. (7-30 and 
7-31): 

A 

* 2 3 

(SV ) = A + D (a.) + D (a.) + D (a.). (for each accelerometer j) 

(7-32) 

A 

( 6 ^c^k = R b + + Q 2 <w T ) 2 (for each gyro j) 


(7-33) 



A 


where (6Y ), and (5ip- ), are estimates determined from the accelerometer 
c k c k 

and gyro calibration processes (as presented in Sections 8 and 9). The 

T 

inputs (a.) and (w ), are provided by components of gravity and test 

X K 


m 


table rates, respectively, averaged over At and corresponding to the 


various test conditions k. Note that 5V = A,’ when a. =.0 and 
a T c b 1 

dm = R’ when w =0. This emphasizes the definitions of A, and R. 
c b T v b b 

as being based on all a_. and w runs; not just those when the inputs are 
zero. 


When more pairs of data are available than coefficients to be 
estimated for a given inertial sensor, filtering is possible to improve 
the accuracy of the estimates for that sensor. By utilizing the above 
regression equations, rather than incorporating them into the basic 
accelerometer and gyro calibration equations (7-1, 7-7 through 7 -9, and 
7-25 through 7-29), different amounts of data under various conditions 
can be handled without having to modify the original estimation equations. 
The regression equations are also convenient for determining just how 
many coefficients are significant and should be included once test data 
becomes available. 

The detailed derivation of the scale factor estimation equations is 
contained in Section 9 and Appendix C. 
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8. CALIBRATION TEST PROCEDURE 


The purpose of this section is to establish the basic procedures 
required to affect calibration of the ISU, Considerations and trade-offs 
in the development of the accelerometer and gyro calibration processes 
are discussed and the basic procedures are presented in sufficient detail 
to support the preparation of detailed step- by- step procedures. However, 
the procedure described in this section is considered to be an initial or 
interim design since certain trade-offs have yet to be made. This is 
particularly true concerning pre-calibration alignment and compensation 
for the various misalignments measured. Furthermore, it is likely that ■ 
the estimation accuracy of some of the calibration terms can be improved 
and time to calibrate reduced by changing some of the ISU orientations 
and/or specifying different combinations of runs and/or using data from 
all inertial sensors from each run (the case if the ISU flight computer is 
used to collect and reduce the data). Finally, the design may have to 
be changed if analysis of inertial sensor test data indicates that- different 
math models should be used. 

Calibration of the ISU gyros and accelerometers may be performed 
at four different levels, as indicated in Table 8-1. The minimum level of 
calibration is designed to estimate accelerometer and gyro biases, mis- 
alignments and scale factor errors under only one set of inputs. In 
addition, the gyro mass unbalance terms (B , B Q and B_) are determined. 
Six different orientations of the ISU are required and the test table is 
driven at + 2°jkec for some of the runs, compensation of the gyro data for 
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Table 8-1 The Four Levels of ISU Calibration. 
Minimum Calibration. 


Accelerometers 


MOA 

MPA 

6V +1 


T 


BA 


V 


6V _1 | D„ , 

c J 1 ’ 2 


Gyros 


From w = 0 


R b 


B t 


B, 


B 


O 


T 

From w = + 2°/ sec 


MOA 


MSA 


6i) + (at+2°/sec)l R 
c I b, 

6$ (at -2° /sec) \ Q. , Q„ 

T C -'12 


BG 


Partial Calibration {above plus the following) 


5V 


,+. 7 


A b'°l 


6V~* 7 o d 
c J 2’ 3 


II 


'SS 


Full Calibration (above plus the following) 




Gyro Scale Factor Calibration 

' + 1 " 

Additional measurements of dip and dip are made for a variety of table 
rates (w ) between + 1 and + 60° /sec, to refine the 3 gyro scale factor 
error coefficients (R , 0 and Q ). 

OX cj 

Note: A matrix format is used in the estimation of the above calibration 
terms, and the terms are arranged in column vectors in the order shown. 
The gyro terms dependent on table rates constitute the last four elements 
of the 13 element vector. 
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the C and C terms is possible, if necessary. The partial calibration 

level requires an additional three ISU orientations with the sensor IAs at 

45° to the vertical for a total of nine positions. This provides estimates 

of two additional accelerometer scale factor terms and the gyro major 

compliance terms C and C cc . Compensation of the sensor output data 

XX bb 

for the C „, C c and C TT3 cross compliance terms is possible, if necessary. ' 
XO Xb IX - * 

The full level of calibration provides estimates of the cross compliance 
terms and requires an additional three ISU orientations, for a total of 
twelve positions. The fourth level of calibration provides a variety of 
table rates to improve the estimation accuracy of the gyro scale factor 
error terms. 

Presentation of the operational calibration procedure is divided 
into three parts. The first, pre-calibration procedures, is concerned 
with measurements of the relationships between the ISU cube, the bubble 
levels, the test table rotary axis and the earth's coordinate frame. The 
second part ISU orientation and calibration data collection, is concerned 
with the process of obtaining inertial sensor output data under the proper 
conditions so that the calibration terms can be estimated. The third 
part, processing of the data to estimate the calibration terms, is con- 
sidered separately in the following section (9). Fig. 8-1 is an operational 
flow diagram that illustrates the total calibration procedure as described 
above. Details of the figure are discussed in the following sections. 

Other operational procedures associated with the calibration 
process which are considered beyond the scope of this study, include 
turn-on/ shut down procedures, maintenance and checkout of equipments, 
and initial installation and calibration of the test equipment. 
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Figure 8-1 Operational Flow Diagram of Total Calibration 



CALIBRATION RUNS 




i 


Caiib. 

GSF 

only? 


| Cahb. 
GSF° 


i 


set w 


Minimum 

Cahbr. 


Table 8-3 


GH 


Integrated 
Minimum Cahb. 
arid GSF Cahb. 


w 

n 


Calibr. 

GSF? 


Partial 
Calibr. (PC) 


(Table 8-4) 


P. C. Complete’ 


£J 


(Table 8-5) 



Calibr. 


Full 


Calibr, 

n 


GSF ’ 


Cahbr. 


GSF? 

■ 


Continue 

Calibration 


Set w 


Gyro Scale 
Factor (GSF) 
Tests 


-.Am 


„itm 




,Gm 


HRunv/ = + 2°/ sec 


(Table 8-6) 

ESTIMATION OF CALIBRATION TERMS 



Record data for 
each test condition 

























DEVELOPMENT CONSIDERATIONS 


8 . 1 


The calibration test procedure was developed in accordance with 
the objectives and ground rules presented in the Introduction (Section 2). 
Additional considerations that were used in establishing the design are 
presented in this section and the general calibration philosophy is discussed. 
Although the Univac study was used as a model, the procedure evolved is 
different in several areas. Not only are the accelerometers different, but 
the test stand has only one degree of freedom, which is constrained to be 
vertical. This automatically required a change in certain test conditions, 
which necessitated further changes to provide adequate sensitivity in the 
estimation process. 

The general calibration philosophy is to use the earth's gravity 
vector and test table rate capability to exercise the ISU in a variety of 
orientations. Having developed the propagation of these inputs (as a 
function of the calibration terms) to the outputs of the gyros and accelero- 
meters (which can be measured), the terms can be estimated. The 
objective was to maximize the effects on the inertial sensor outputs of 
the calibration terms, in order to improve sensitivity of the estimation 
process, while minimizing the effect of various error sources. Further- 
more, an attempt was made to choose test conditions that would cause 
instruments outputs that were a function of only a few of the calibration 
terms for any one run. This has the advantage of requiring fewer runs 
to estimate any particular calibration term, thereby simplifying the 

and combinations of test conditions 
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data reduction effort and providing flexibility in determining a term without 
a large number of runs. It also has the advantage of minimizing the effect 
of random gyro drift changes between runs. For certain test conditions a 
calibration term not to be estimated affects the instruments outputs to 
an extent that compensation is required. This case occurs in two forms: 
in the first, errors in the knowledge of the calibration term has ah insig- 
nificant effect, whereas in the second, the calibration term error may be 
significant. 

The calibration procedures were designed to be performed in a 
series of groups, such thatobtaining additional data would allow estimation 
of addition calibration terms. In the case of the gyros, the I A to cube 
-misalignments can be determined directly from a pair of runs, as can 
the scale factor error at a given input rate. Once the required data for 
a minimum calibration has been obtained, only three additional runs are 
required to complete a partial calibration (see Table 8-1). A full calibration 
requires four more runs, as explained below. 

The pre-run alignment procedure shown in Fig. 8-1 was developed 
with the following approach as a basis. The effect of gyro and accelero- 
meter iA misalignments is to cause outputs that are proportional to the 
accelerometer IA to jg nonorthogonality and the gyro IA to table rotary 

■ju 

axis nonorthogonality , respectively. To simplify the alignment activities 
before each run, the bubble levels are used to determine the tilt of the 


the effect of rotary axis to vertical errors cause a negligible cross 
coupling of earth rate to the gyro IAs 
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rotary axis from g and then the TAP gimbals are positioned using the 
gixnbal angle readouts so that the bubble levels are aligned with the rotary 
axis. In this way the gyro IAs are referenced to the bubble level coordinate 
frame at the beginning of each run and the accelerometer output data 
can be compensated for changes in the rotary axis to vertical alignment, 
thereby also referencing the accelerometer IAs to the bubble level 
coordinate frame. The transformation matrix between the ISU cube and 
the inertial sensor IAs (which is the end result desired) is then calculated 
separately using the cube to bubble level calibration made at the time the 
ISU is installed in the TAF (as shown in Fig. 8-1 and explained in Section 
8. 2 below). Furthermore, compensation for earth rate coupling to 
the gyro IAs is accomplished using the bubble level to gyro IA transforma- 

"j- ■a" 

tion matrix determined above . 

Most of the calibration runs were chosen with the TAF inner 
gimbal axis nonvertical in order to allow use of the bubble levels in 
determining the rotary axis tilt from vertical. A total of 6 positions with 
IAs either nominally horizontal or vertical are required to perform a 
minimum calibration. This effectively eliminates the effect of cross 
compliance terms, and for the IAs vertical, the effect of IA misalignments 
to the ISU cube. 


this alignment is performed at assembly of the bubble levels, TAF and 
SATS, when the equipment is aligned, calibrated and zeroed. 

•J* vO 

this is required for all gyro runs except those required to measure the 
transformation matrix itself 
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The gyro IA misalignments are measured by driving the test table 
at + 2°/ sec and observing the outputs of those gyros whose IAs are nomi- 
nally normal to the rotary axis (Rp. The process is repeated twice for 
each gyro; once with the OA normal to R^ and once with the SA normal 
to R!^. A higher rate is not recommended since rotation of a gyro about 
its OA at rates greater than approximately 2°/ sec will cause a potentially 
excessive gyroscopic torque on the gyro float about the IA, which may 
cause additional gyro drift rates that are significant. Lower table rates 
would require more run time to achieve the same calibration precision. 

Data from the gyros are taken only for full revolutions of the test 
table, for two reasons. First, the effect of certain components of earth 
rate tend to cancel out over full revolutions and second, the cost of the 
SATS can be reduced by eliminating the requirement for a continuous, 
precise readout of table angle. 


The partial calibration procedure is designed primarily to provide 
a change in accelerometer input accelerations so that the first and third 
order scale factor error terms can be separated. Both plus and minus 
45° orientations to horizontal are chosen for each accelerometer, even 
though only one of these inputs would suffice, in order to minimize the 
potentially excessive error that may occur for inputs having the opposite 
sign from that chosen. The 45° orientation also is used to estimate the 
major gyro compliance terms. 


The remaining four calibration procedures are designed\to obtain 
data to estimate the cross compliance terms for both the gyros and 



accelerometers. It is likely that the combination presented in Section 8. 3 
(Tables 8-3, 8-4 and 8-5), however, is not an optimum one in that the 
last run could probably be eliminated by choosing a different combination 
of test conditions. A total of six combinations were considered, all with 
the constraint of only observing four out of the six inertial sensor outputs 
for any one run. If the flight computer is made available, all six outputs 
could be utilized from each run and either the number of runs could be 
reduced or the precision of calibration term estimation could be improved. 

8. 2 PRE-CALIBRATION OPERATIONAL PROCEDURE 

There are many ways that the precalibration alignment measurements 
can be made. As mentioned earlier, the method described in this section 
should be considered only illustrative of the type of measurements required 
since detailed trade-off studies have not yet been performed. Regardless 
of which specific method is adopted, the requirements of the precalibration 
alignment procedures can be divided into the following four major categories, 
as illustrated in Fig. 8-1: 

• alignment and calibration of the bubble levels (which are 
located on the inner gimbal of the TAF) to the inner gimbal 
axis and to the SATS rotary axis (Rj), including zeroing of 
the TAF readouts 0 and 0 . This is normally done during 

O ”E 

test equipment installation and checked periodically thereafter 
depending upon the stability characteristics of the TAF and 
SATS. 
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• alignment of the SATS rotary axis to vertical and the TAF 
outer gimbal axis to East. The alignment to vertical 
consists of a nominal physical orientation of the rotary 
axis using leveling screws and then a precise measure- 

. ment of the residual misalignment (for the scheme con- 
sidered in this section). Depending upon the particular 
operational procedure finally adopted, it may be necessary 
to physically reduce the residual misalignment to a 
negligible amount. Alignment to East is done by adjusting 
the test table angle (0 ). 

a 

• Calibration of the ISU optical cube to the bubble level 
coordinate frame, in the form of precise measurements 
of the misalignments. This measurement must be made 
every time the ISU is mounted in the TAF. 

• Alignment of the ISU relative to the SATS rotary axis and 
the earth's coordinate system immediately prior to each 
run (m). For the scheme considered in this section, this 
is done using the bubble levels and TAF gimbal readouts 
to physically orient the ISU and measure certain critical 
misalignments. 

In addition to these precalibration alignment procedures, it is necessary 
to set up the data collection system, turn the ISU on, etc. , before an 
actual calibration run can be initiated. The operational, requirements 
of each of the above precalibration procedures is discussed next, and 
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in conjunction with the hardware characteristics presented in Section 6, 
form the basis for the support equipment errors used in the system error 
analysis of Section 10. Again, the particular method described here is 
only illustrative of the measurements and accuracies required since other 
variations are feasible and have yet to be included in the calibration 
trade-off analysis. 


The initial alignment, calibration and zeroing of the bubble levels, 
TAF and SATS are done at the time of test equipment installation. The 
process consists of mechanically adjusting the bubble levels so that they 
represent the TAF inner gimbal coordinate frame (I) (see Fig. 5-3 and 
Table 5-2). This may be done by noting the TAF gimbal angle changes 
as the bubble levels are nulled for various orientations of the TAF. 
Ideally, the bubble levels should be either parallel or perpendicular to 
the inner gimbal axis. The inner gimbal angle (0 ) is zeroed such that 
the bubble level is parallel to the outer gimbal axis. The outer 
gimbal angle (0 ) is zeroed such that the inner gimbal axis is parallel 

u 

to the rotary axis. The values of the TAF gimbal angles when the bubble 

levels are normal to the rotary axis are labeled (0) and (0) and 

3 z 4 z 

are recorded. Nonorthogonality checks between the 0 , 0 and 0 axes 
are recommended to assure acceptable performance, and if compensation 
is to be applied in the data processing equations, the nonorthogonalities 
must be measured and recorded. 


The alignment of the SATS rotary axis to vertical is accomplished 
by using the leveling screws, as indicated by the TAF gimbal angle 
changes required to null the bubble levels when the test table is rotated 180°. 
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Alignment of the TAF outer gimbal axis to East is accomplished using 
the autocollimator and adjusting the test table angle separate 

survey is required to establish East, as a reference for the autocolli- 
mator. The device for indicating this position of the table is secured 
and this angle of the table is labeled (0 O ) . 

Calibration of the ISXJ cube to bubble level coordinate frame 
misalignment is performed each time the ISU is mounted within the TAF. 
The procedure is to place each of the three bubble levels in the nominally 
horizontal position, in sequence. In each instance, the bubble level is 
nulled by adjusting 0 and 0 . With the autocollimator referenced to 

«J ut 

level, the angle the visible ISU cube face makes with respect to level 
is measured. This angle is then a measure of the misalignment between 
the bubble level pointing in the direction of the autocollimator and the 
ISU cube face observed. The process is repeated for each bubble level, 
and assuming the ISU cube faces are all orthogonal, all six elements of 
the bubble level to mirror coordinate frame transformation matrix (T ) 
are established. It is estimated that approximately 1. 5 hours will be 
required to mount the ISU, setup the autocollimator and calibrate the ISU 
cube/bubble level misalignments. This is based on the following time 
allocation. 


l. e. , 


ISU cube 
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• Mount ISU 

• Setup autocollimator and level 

« Perform optical measurements (3 times) 


15 minutes 
20 


45 

® 4 minutes to index ISU 

® 6 minutes to measure 

9 5 minutes to complete records 

15 minutes per ISU position 80 minutes = 1, 33 hrs. 

An additional estimated 15 minutes is required to set up the data collection 
system;, turn the ISU on and prepare to start the calibration runs. 

The alignment of the ISU prior to each calibration run (m) is 
performed using the bubble levels and the TAF gimbal angle readouts. 

However, since the alignment procedure is part of the operational calibra- 
tion procedure, it will be explained in detail in the next section (8. 3). In 
Fig. 8-1, the procedure is shown as a precalibration activity, but separate 
from the other alignment procedures because it is repeated for each test 
condition (m). 

8. 3 CALIBRATION OPERATIONAL PROCEDURE 

An operational flow diagram of the calibration procedure is illustrated 
in Fig. 8~1. The various options are indicated, and the heavy line repre- 
sents the one that would normally be followed if all four levels of ISU 
calibration (as listed in Table 8-1) were to be performed. Although the 
pre-run alignment procedure is included as a precalibration activity, it 
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is required before each run (m), except for the Gyro Scale Factor (GSF) 
tests. 


The pre-run alignment procedure is designed to properly orientate 
the ISU relative to the SATS rotary axis and the earth's coordinate system, 
immediately prior to each run. The five step procedure is outlined in 
Table 8-2 and provides a way to compensate for short-term tilts of the 
rotary axis from vertical. If the physical tilt is excessive (say greater 
than . 3 min), the releveling procedure is initiated, as described in Section 
8. 2, in order to minimize errors in the gyro output data. At the end of 
the procedure, the bubble levels are normal to the rotary axis and the 
misalignment of the rotary axis to vertical is known. The latter is 
necessary for compensation of the accelerometer output data. 

The detailed operational test procedures to perform the four levels 
of ISU calibration, as listed in Table 8-1, are contained in Tables 8 -3 
through 8-6. The tables are organized to indicate the ISU orientation 
required for each test condition (m); both in terms of directions of each 
inertial sensor coordinate frame and the corresponding SATS and TAF 
gimbal angle indications. The time estimated to set up each run is 
tabulated. Also tabulated are the data to be collected and the elapsed 

HI 

time (At ) of each run. The runs within any one table may be made in 

T 

any order; however, data taken with w = + 2°/sec provides misalignment 
information that may be helpful in subsequent runs. 


The G.. and A. column headings refer to the particular gyro and 
accelerometer pulse outputs to be counted, A0 refers to the test table 


pulse output, and TAF refers to 0 


3 a 


0 gimbaf angle readouts. 

rr 
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Procedure 


® 2 >z 


1. Position TAF to pre-run align orientation 


W 2>z <*3>” 


2. Null the two bubble levels that are 
horizontal, by adjusting 0o and 0^. 
Record the gimbal readouts as 

(0 3^ and (0 4 } ? 


, ( 2 ) 


3. Check for excessive tilt of rotary axis 

4. If either (S0;g) m or (60 N ) m is greater 
than . 3 min, relevel SATS, as described 
in Section 8. 2, and repeat step 2 above. 


, . , .m , 

^ 2 z ^3 z ^4 z 


5. Index ISU by positioning the TAF gimbals 
as follows: 

(0 o ) m = (0 o ) ni + value of 0 from zeroing 

3 z 3 n 3 

process 

(0 ,) m = (0„) m + value of 0 , from zeroing 

4 z 3 n 4 

process 

where subscript n refers to the nominal 
value. 


Notes: (1) The ISU should not be rotated about any axis at a rate greater 
than 60°/ sec, during any of the positioning procedures, in 
order to prevent gyro loop saturation. 

(2) Check for tilt (60 and 60,,) outrof-tolerance by computing 

-hi N 

<^E> m * <Vz ' ( *3>“ 

' <*4>z + 


Table 8-2 Pre-Run Alignment and ISU Indexing Procedures 
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Perform the following after 
completing the pre-calibration procedure; 



Subtotal 


10 


81 minutes 


6 ) =90 » actual value of 0, to be recorded from Pre-run Alignment using 
3r r 3 

Bubble Levels when (fij __ = 90“ 


^ I 1.73 I !• 35j hours 
Minimum Calib. 3.' 08 hours 

Pre-Calib. 1..58 

Total 4. 57 hours 


) =90 s 90“ + value of from zeroing process. 

w 3z z 3 

Similarly for (0^, and other nominal values of « 3 and« 4 

Table 8-3 Operational Test Procedure for Minimum Calibration 
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Perform a Minimum Calibration plus the following 


m 
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(0 ) “90 - actual value of 0 to be recorded from 
3 r r o 

Pre-run Alignment using Bubble Levels, when (0^) = 

(0„) = 90 = 90® + value of 0_ from zeroing process. 

3 z z 3 

Similarly for ( 0 , ) , (0 J and other nominal values of 0 _ and 0 . 

4 r 4 z 3 4 

Table 8-4 Operational Test Procedure for Partial Calibration 
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Pre-run Alignment using Bubble Levels when (0 3 ) n - 90® 
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Similarly, for(0.) , ( 0 .) and other nominal values of 0„ and 0 . 

4 r, 4 Z o 4 


Partial Calibration 
Subtotal 

Pull Calibration 
Pre- Calibration 
Total 


60 70 

140 111 


f 200 181 

L 3.33 3. 02 


minutes 

hours 


S. 35 hours 

1.58 

7. 9 hours 


Table 8-5 Operational Test Procedure for Pull Calibration 
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continue for w^'=±8, ±20 and ±60°/ sec (m=45 + , 45“ 46 + , 4S“, 47* 47“> 



continue for w T =±4, ±8, ±20 and ±60‘7sec (m=64 + , 64“, 65+. 65", 66*, 66", 67* 67") 


8-6b. Calculation of Time Required to Calibrate 


If w'^'=±2*/sec has not been 
run (as part of other calibra- 
tion tests), must add an ad- 
ditional . 4 hr per gyro. 


Table 8-6 Operational Test Procedure for Gyro Scale Factor Error Calibration 



Table entries are times required to perform the operations indicated. 
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The "Pre-run align" and -"-Index ISU" procedures referred to in the 

tables are as discussed above and summarized in Table 8-2, The "Set 
T 

w " procedure is the process of adjusting the test table rate drive to the 
T 

w value indicated. 

This is done by checking the time required for the table to rotate 
through a known angle (90°, 180° or 360°, depending upon the table readouts 
provided). If the rate is not within + . 5% of the desired value, the rate 
drive is readjusted and the timing process repeated. The data collection 
system is described in Sections 6. 1 and 6. 5. 

The operational test procedure for performing a minimum calibra- 
tion is specified in Table 8~3. A total of 12 runs for 6 different orientations 
of the ISU are required to collect the necessary data. This includes 
operating the test table at plus and minus 2° /sec for each of 3 of the ISU 
positions. An estimated 3. 1 hours is required to complete the runs, which 
in addition to 1. 58 hours to perform the precalibration activities', requires 
a total of 4. 7 hours to calibrate. 

To perform a partial calibration, an additional 3 runs for 3 ISU 
orientations at 45° to the vertical are required, as indicated in Table 8-4. 
This consumes an additional 1. 1 hours which when combined with the 4. 7 
hours for the minimum calibration, requires a total of 5. 8 hours to 
complete a partial calibration of the ISU. 


as determined in Section 8. 2 
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The additional procedures to perform a full calibration are 
specified in Table 8-5. The new ISU positions are also at 45° to the vertical 
since the cross compliance terms are to be estimated. Each of the first 
3 runs are repeated twice since data from all inertial sensors are required; 
yet only 4 outputs can be processed at any one time. The additional data 
is required to avoid excessive sensitivity in the estimation of the calibration 
terms. Test condition No. 13 is required to provide additional data from 
accelerometer No. 2 in order to estimate one of the cross compliance 
terms. It is likely that a different combination of test conditions could 
eliminate the need for the extra run, particularly if data from each inertial 
sensor could be collected for each run (the case, for example, if the flight 
computer were to be used to collect data). As shown in Table 8-5, an 
additional 2. 1 hours is required to collect the data. Combining this 
with the 5. 8 hours for a partial calibration, a total of 7. 9 hours is 
required to complete a full calibration of the ISU. 

Table 8- 6 specifies the procedures and times required to operate 
the test table at various rates between + l°/sec and + 60°/sec, in order 
to provide data for the gyro scale factor error regression analysis dis- 
cussed in Sections 7. 3 and 9. 3. The rates listed are only suggestions 
and more or less data may be taken depending upon time available and 
number of samples desired to be included. A total of 4. 8 hours is required 

to run 5 pairs of plus and minus rates for each gyro. If only gyro scale 

T 

factor calibration is desired, w - + 2°/ sec is also required (at . 4 hr. 
per gyro) and the total calibration time becomes 6. 0 hours. 

If gyro scale factor calibration is to be done in conjunction with 
the other calibration runs, it is recommended that they be combined, as 
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indicated in Fig. 8-1. In this way, the first 2 steps in each of the test 
condition runs m = 4, 5 and 6 in Table 8-6 can be eliminated, since they 
are already performed as part of the minimum calibration (see Table 8~3). 
Therefore, only 4. 4 hours is required in addition to the times listed in 
Tables 8-3, 8-4 or 8-5 to perform the calibrations indicated plus a full 
range gyro scale factor calibration. 
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9. ESTIMATION OF CALIBRATION TERMS 

The purpose of this section is to derive the equations to be used in 
estimating the calibration terms. The algorithms presented are based on 
making a least squares fit to the data obtained whenever more measure- 
ments are used compared to the number of terms to be estimated. When 
the number of measurements and terms to be estimated are equal, the 
least squares equations degenerate into a deterministic solution. The 
pre -calibration calibrations are also presented. 

Both the accelerometer and gyro calibration equations are arranged 
in the following form; as shown in Section 7: 

M^/At* 1 s M 111 = EpY. +6M m (7-1) and (9-1) 

3 3 3 — 3 3 

in which the various terms are defined by Eqs. (7-7), (7-8), and (7- 9) for 
the accelerometers and Eqs. (7-25) through (7-29) for the gyros. Tables 
B-2 through B-7 in Appendix B contain the results of applying the calibra- 
tion equations to each test condition (m) for each inertial sensor. The 
equations are organized in the following matrix form; 

M. = H. y. (9-2) 

~ 3 3 3 

where 

H. = a matrix whose rows are the rf 11 of Eq. (9-1), for (9 _ 3) 

J the various test conditions (m) 3 

and 

M . = a column vector whose elements are the M ^ from (9-4) 

3 Eq. (9-1) 3 
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The matrix H. can assume various dimensions, depending upon the mea- 
* .3 

surements (M) used and the calibration terms in the vector Y that 

* 

propagate into M. For the case when H. contains the same or more rows 

rj> .1 

than columns, and (H_. H ' ) is nonsingular, it can be shown that the least 
squares solution of Eq. (9-2) is as follows: 


A q 1 _i q * 1 • 

Y. = (H. H.) H. M. (9-5) 

“3 3 3 3 — 3 

A 

where Y. is the estimated value of Y.. This equation was applied to the 
3 3 

calibration equations in Tables B-2 through B-7, for the te^t conditions 
specified in the various calibration procedures of Tables 8-3, 8~4 and 
8-5. Results are presented in Appendix C and discussed in Sections 9. 2 
and 9.3. 


In the case of the scale factor error coefficients, defined in 
Eqs. (7-32) and 7-33) of Section 7. 3, the general least squares solution 
is as follows: 


(F 4 Z 0 - F 2 Z 2 ,/(KF 4- F 2> 

(9-6) 

< F 6 Z 1 - F 4 Z 3 )/(F 2 F 6 - F I> 

(9-7) 

< K V F 2 Z 0 >/(KF 4- F 2> 

(9-8) 

'Vs' F 4 Z X )/(F 2 F 6- F 4> 

(9-9) 
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where the C's correspond to the respective A^, D, R^ and Q coefficients 

in Eqs. (7 -3 2) and (7-33) and the F and Z terms are functions of the 

inertial sensor inputs and indicated outputs. These terms are defined 

below for the accelerometers and gyros. The term K is the number of 

data samples used in the least squares fit. The above equations are based 

on the assumption that inputs to the inertial sensors occur in equal and 

opposite pairs (e. g. , + g//~2 for the accelerometers and + 4°/sec for the 

gyros, as well as other pairs of inputs). The equations also assume that 

A! and R* are used, corresponding to zero input to the accelerometer and 
b b ' 

gyro, respectively. 


The estimated inertial sensor misalignments transformation 
A BA A B G 

matrices T and T are determined from the estimated accelerometer 

A/ 

and gyro nonorthogonalities, as contained in the Y vectors of the calibra- 

A 

tion estimation equations. However, these Y elements are the misalign- 
ments between the LAs of the inertial sensors and the bubble level 

A A 

coordinate frame (T • and T ), as discussed in Appendix B (following 
Eq. B-7); and not the misalignments between the IAs and the ISU cube 

A a bq 

(mirrors) which are elements of the desired T and T matrices. 
Having measured the elements of T’ M ' L , as described above under the 

A A gQ. 

pre-calibration calculations (Section 9. 1), the T and T matrices 
can - be determined as follows: 


^ba = t bm t ml ( ^al } -i 

and 

£BG = t bm t ML ( ^GLj-1 


(9-10) 

(9-11) 
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, BJV[ _ . . _ . * ^ A T i Q\ > 

where T and T are given m Eq. (B“6) and T and T are 

* * 

composed of the estimated nonorthogonalities of the three accelerometers 

and three gyros with respect to the babble level coordinate frame. The 

various matrices are defined in Table 9-1 and are expanded in the following 

a A B G 

sections to estimate the elements of T and T , the transformation 
matrices ultimately desired. 


Matrix 


T 


ML 


T 


BM 


AL 


T 


GL 


,BA 


J3G 


Description of Transformation Matrix 

misalignment of ISU cube (mirror) coordinate frame 
relative to bubble level coordinate frame; elements 
defined in Eqs. (B~6) and (9-18). 

definition of body axes in terms of ISU Cube (mirror) 
coordinate frame 

misalignment of coordinate frame defined by IAs of 
the accelerometers relative to the bubble level coordinate 
frame; elements of matrix are NOA and NPA for each 
accelerometer 

misalignment of coordinate frame defined by IAs of the gyros 
relative to the bubble level coordinate frame; elements of 
matrix are NOA and NSA for each gyro 

matrix used to transform accelerations sensed along IAs 
of the accelerometers into accelerations in the body axes 
coordinate frame; elements of matrix are MOA and MPA 
for each accelerometer 

matrix used to transform rates sensed about IAs of the 
gyros into rates about the body axes coordinate frame; 
elements of matrix are MOA and MSA for each gyro 


Table 9-1 Transformation Matrices Used to Calculate Inertial 
Sensor to ISU Cube Misalignments 
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PRE-CALIBRATION CALCULATIONS 


9. 1 


Calculations in this category consist of processing data obtained 
before the actual calibration runs are made. The calibration procedure 
was arranged such that only two calculations had to be made before each 
run, to determine the tilt of the test table rotary axis from vertical to 
assure that it was not excessive. The other calculations are involved in 
determining the bubble level to ISU cube transformation matrix (T ML ). 


The calculations to determine the tilt of the rotary axis are 
presented in Table 8~2 and repeated as follows: 


.. .m ,m . .m 

(6 V = ^3 z ' ^3 r 


(9-12) 


< 6 v m = ( V" 


(9-13) 


where 


HI m 

(0 ) = (0„) + value of 0 from zeroing process 

o z on «j 


(9-14) 


(0 ,) m = (0.) :rn + value of 0 from zeroing process (9-15) 

4 z 4 n 4 

HI 

(0 3 )r = TAF outer gimbal readout when bubble levels (9-16) 

are nulled in the pre-run align orientation, 
just prior to run (m) 

(0^)^= TAF inner gimbal readout when bubble levels (9-17) 

1 are nulled in the pre-run align orientation*, just 
prior to run (m) 

* the pre-run align orientation is when 0 = 90°, 0 = 0°, as discussed in 

Section 8. 3 and Table 8 -2. 4 
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and subscript n refers to the nominal values of (0 o ) m and (0 > m for run 

. 3 4 

(m). The values of 0 and 0 from the zeroing process are those values 
that align the bubble levels with respect to the rotary axis. The 60 
and 60^. terms are used in the accelerometer estimation equations, as 
shown below. 


The transformation matrix (T^^) contains 6 elements, 3 of which 
are independent and must be measured separately. This is done using the 
bubble levels and the auto collimator at the time of TSU mounting as 
explained in Section 8. 2. The T^^ matrix is defined in Eq. (B _ 6) and 
Table B~1 and repeated as follows: 



13* —1 

1 

% 

L 

ML 





ppJ.VJ.J_f _ 

“22 — 1 

*"—21—2 

1 



-1 

-M • 

M ■ 

L 



— 33 —2 

— 32 



(B-6) and (9-18) 

in which the elements are defined in the form = nonorthogonality 

of M a about with respect to L measured about the positive L ^ 
axis in accordance with the right-hand rule. The matrix elements are 
determined from the following transformation equation: 


t ml = t ms t se (t le } -i 


( 9 - 19 ) 


where the transformations are defined in Table 9-2. The bubble levels 
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Matrix Description of Transformation Matrix 

misalignment of ISU cube (mirror) coordinate 
frame relative to bubble level coordinate 
frame; elements defined in Eqs. (B _ 6) and (9~18). 

matrix of measurements of ISU cube, expressed 
in autocollimator coordinate frame 

SE 

T matrix relating autocollimator coordinate frame 

relative to earth coordinates (specifically to 
vertical); assumed equal to the identity matrix 
when the autocollimator is leveled. 


LE 

T matrix representing bubble level coordinate frame 

relative to earth coordinates (specifically to vertical); 
assumed equal to the identity matrix when the bubble 
levels are nulled 


Table 9~2 Transformation Matrices Used in Alignment of ISU 
Cube and Bubble Level Coordinate Frame 


9-7 



SE 

and autocollimator are adjusted relative to vertical such that T and 
LE 

T are essentially equal to identity matrices. Therefore, the auto- 
collimator to ISU cube measurements provide elements of T^^ directly. 
As mentioned above, only 3 measurements are required since the coor- 
dinate frames are essentially orthogonal and therefore the following 
relationships exist: 


21* —2 

= “11—3 

(9-20) 

13* — 1 

= ~— 33* —2 

(9-21) 

—32 —3 

" —22—1 

(9-22) 


Thus all elements of T' ML< in Eq. (9-18) are determined. 

9. 2 ACCELEROMETER CALIBRATION ESTIMATION EQUATIONS 

The calibration estimation equation [Eq. (9- 5)] was applied to the 
accelerometer calibration equations in Tables B-2, B~3 and B~4 to 
provide the accelerometer calibration estimation equations. The results 
are tabulated in matrix form in Tables C~1 , C~2 and C-3 of Appendix C 
for the minimum, partial and full levels of calibration respectively. It 
will be noted that the equations do not include compensations to be applied 
to the M measurements in order to account for gravity coupling due to 
significant misalignments and tilt of the rotary test table axis from 
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vertical. These corrections can be specified when the final pre-calibra- 
tion alignment scheme is defined. Compensation for C ^ and C^p in 
the partial calibration estimation equations is shown in Table C~2. 


The accelerometer scale factor error estimation equations were 

A ■ 

determined from Eqs. (9~6) through (9-9), using the estimates from 

a . +1 *P- -1 A - +.7 

the accelerometer calibration equations (viz., <5V , 6V , 6V ’ , 

. c c c 

" * 7 

6 V ) and the term A! . The results are tabulated in Table C~4 of 

b 

Appendix C. 


a 

The elements of the T transformation matrix, which are a 

function of the misalignments of the IAs of the accelerometers relative 

to the ISU cube, are computed using Eq. (9-10) as a basis. The matrix 
a 

T in Eq. (9~10) is composed of the estimated nonorthogonalities 
(in Y_f~, Y and from Tables C~l, C~2 and C~3), as follows: 

Zi tj 


t al 


NPAj 

-1 

-NOA 

A 

A 


n °a 2 

NPA 2 

1 


A 

A 

1 

-NPA 

NOA 


(9-23) 


A BA 

and writing T in terms of the inertial sensor 

A A 

misalignments to be estimated (viz. , MOA and MPA), Eq. (9-10) is' 


a 

Upon inverting T , 


expanded to yield: 
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-NOAg 1 

A A 

-NOA +NPA X 
1 -NOA 2 

(9-24) 


Upon performing the multiplications indicated and equating corresponding 

A BA 

elements in the resulting matrices, the elements of T are computed as 
follows: 


A A A 


MOA 1 

= NOA x 

~ —i i —3 

(9-25) 

A 

A 

A 


MPAj 

= NPA 2 

- M 13 * 

(9-26) 

A 

A 

A 

(9-27) 

MOA 

= NOA 

“ M * L, 

2 

2 

— 22 —1 

A 

A 

A 


MPA 2 

= NPA 

u 

+ n 21 -k 2 

(9-28) 

A 

A 

A 


MOAg 

= NOAg 

+ “32—3 

(9-29) 

A 

A 

A 


MPA 3 

= npa 3 

+ — 33 " —2 

(9-30) 


9-10 



9. 3 GYRO CALIBRATION ESTIMATION EQUATIONS 

The gyro rate errors due to scale factor errors and the gyro 
misalignments are calculated in a deterministic fashion as implied in 
the first four rows of the calibration equations tabulated in Tables B~5, 
B-6 and B~7. In this way a minimum number of runs are required to 
estimate the terms and the calculations are relatively simple. The 
results are tabulated in Table C~5. 

G 

The remaining gyro calibration terms in Y are computed by 
applying the least squares equation [Eq. (9-5)] to the gyro calibration 
equations in Tables B-5, B-6 and B~7, The results are tabulated in 
matrix form in Tables C-6, C-7 and 08 of Appendix C for the minimum, 
partial and full levels of calibration, respectively. It will be noted that 
the equations do not include compensations to be applied to the M mea- 
surements in order to account for earth rate coupling due to NOA and 
NSA, the elements of (rpGM^ML^ These corrections can be specified 
in detail when the final pre-calibration alignment scheme is defined. 
Compensation for C^ and Cgg in the minimum calibration estimation 
equations is shown in Table C~6 and a similar approach is used in Table 
C-7 to compensate for C and C^, in the partial calibration estimation 
equations. 


The gyro scale factor error estimation equations were determined 

A. 

from Eqs. (9~6) through (9~8), using the 60 estimates from Table C-5, 
(viz. , 60 , 60 pairs for various table rates w ) and R, 1 as estimated 

Aq C C b 

in Y . Since a third order term is not modeled for the gyro, the term 
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in Eqs. (9-7) and (9-8) can be allowed to approach infinity (which eliminates 
the third order term) and the resulting estimation equations are as 
tabulated in Table C-9 of Appendix C. If the third order term is desired, 
the F term is retained, and Eq. (9-9) provides the desired estimate 

® A 

(note that in this case would be given by Eq. (9~7)). 

A BG 

The elements of the T transformation matrix, which are a 

function of the misalignments of the IAs of the gyros relative to the ISU 

a ql 

cube, are computed using Eq. (9-10) as a basis. The matrix T in 
Eq. (9-10) is composed of the estimated nonorthogonalities from Table 
C-5, as follows: 


rJ GL = 


N f l 

-NOA r 


NOA. 


-nsa 2 1 

~NSA 3 NOA 3 


a qj i B G 

Upon inverting T , and writing T in terms of the 

A A 

misalignments to be estimated (viz. , MOA and MSA), 


expanded to yield 


(9-31) 


inertial sensor 
Eq. (9-11) is 



A 

A 

* 





A 

A 


1 

-MOAj 

“MSAj 


1 

0 

0 


NSA g 

-NOA 3 

1 

A 


A 





A ML 


A 

A 

MSA 2 

1 

moa 2 

= 

0 

1 

0 


1 

-NOA 

-NSA 1 

A 

A 







A 


A 

MSA 

-MOA 

O 

1 


0 

0 

-1 


NSA 2 

1 

NOA 

o 


— 


,, , 






C-i 


(9-32) 
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Upon performing the multiplications indicated and equating corresponding 

AgQ 

elements in the resulting matrices, the elements of T are computed 
as follows: 


A A A 


MOA 1 

= NOA 1 

-*U' L 3 

(9-33) 

A 

A 

A 

(9-34) 

MSA^ 

= NSA 1 

+ — 13 L 1 

A 

A 

A 


MOA 

u 

= NOA 

+ —22 L J 

(9-35) 

A 

A 

A 


MSA 2 

= NSA 2 

' —21 L 2 

(9-36) 

A 

A 

A 


MOA 3 

= noa 3 

+ — 3 2 L 3 

(9-37) 

A 

A 

A 


MSA 3 

= NSA 3 

+ M S 3 L 2 

(9-38) 
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10. ERROR ANALYSIS 


10.1 GENERAL FORM OF ERROR EQUATIONS 

The error analysis of the calibration process is based on the general 
form of the calibration equation developed in Section 7. Eq. (7-1) defines 
the relationships between the various calibration terms (Y) and adjusted 
measurements (M), based on the average of the inertial sensor outputs 
over the calibration run time (At). The sensor outputs are compensated 
for significant and known inputs to the instruments such that 

M m =■ H m ■ Y. +&M m ( 7-1) & (10-1) 

3 ~3 “ 3 3 

where Siyf 11 is the total of all errors in the quantity M^ 1 . The error analysis 
consists of identifying all of the significant error sources that contribute to 
6M*j a and then determining the effect this has on the precision with which 
the calibration terms can be estimated. Since the terms are estimated 

« tn 

using Eq. (9-5)., the problem is to determine the effect of 6 m . on 
A J 

Y. in Eq. (9-5). 

J 

The errors are all assumed to be random variables with zero 
means. The analysis is generally performed on a statistical basis in 
terms of variances. Following this approach, it can be shown that the 
covariance of errors in the estimation of the calibration terms is given 
by the following matrix equation: 

E[Y‘Y T ], = (H T H). _1 H T E[6M.6M T J. H.(H T Hh 1 (10-2) 

3 3 3 — — 3 3 3 
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where 


6M. 

E[6m-6m T ] 


A A m 

E[Y- Y ]. 

*3 


• xn 

column vector of errors 6M for instrument j 
expected value of the matrix of error correlations 

error covariance matrix of M., composed of 

m " * m 

e elements, the sum squared errors in 6.M. : 
ss j 

assumed to be a diagonal matrix (10-3) 

A 

error covariance matrix of Y.; in which the 

3 

square root of the diagonal elements are the 
standard errors (SE) of each of the estimated 

A 

calibration terms in Y. (10-4) 

— J 


and H. is defined in Eq. (9- 3) as a matrix whose rows are the IT^ 1 of 
3 ~~ 3 

Eq. (9-1), for the various test conditions (m). The error analysis now be- 

m • Ijl 

comes a matter of defining E[6M* 6M ] . for each gyro and each accelerometer. 


The inertial sensor error sources (6M m ) are identified by taking 
partial derivatives of the general calibration equation (7-1), as defined for 

the accelerometers by Eqs. (7-7) thru (7-9) and for the gyros by Eqs. (7-25) 

* 

thru (7-29). Errors in the compensation corrections within the M calcula- 
tions, as given by Tables B-2 thru B-7, C-2, C-6 and C-7, must also be 
considered. The general categories of the error sources are associated 
with inertial sensor output pulse counting, sampling time, IA, OA, SA and 
PA orientation errors, environmental acceleration and rate effects (viz. , 
linear and rotational vibration of the test stand), inertial sensor noise 
(internal and rebalance loop effects), compensation errors, second order 
terms and other neglected error sources. Additional error sources in the 
calibration of the gyros are due to test stand rate errors. 
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A gj^ A gQ 

In the case of determing T and T , using Eqs. (9-24) thru (9-30) 

and (9-32) thru (9-38) respectively, errors in measuring the elements of 
A ML 

T - must be included~to determine the errors in estimating the misalign- 
ments between the inertial sensor IAs and the body axes. 


10. 2 PRE -CALIBRATION ALIGNMENT ERRORS 

One of the errors in this category is associated with the measure- 
ment of T^"*L Since details of the procedures have not yet been defined 
the error analysis of this portion of the calibration procedure must be 
postponed. The effect of the errors, however, is not difficult to incor- 
porate later since they occur at the end of the analysis, in the determina- 

a da A g(q 

tion of T and T as shown in Eqs. (9-25) thru (9-30) and (9-33) thru 
-•(■9-38), respectively. 

ML 

The effect of errors in T on compensation for earth rate coupling 

in the gyro calibration process are of second order and completely negligi- 

“ ML 

ble. The accelerometer calibration process does not use T compensa- 

a BA 

tion, except as noted above in the computation of T 

Pre -calibration alignment errors due to use of the bubble levels 
and indexing of the ISU have not been defined in detail pending final defini- 
tion of the pre -calibration alignment procedure and significant error 
sources. The primary effect is expected to be associated with compensa- 
tion of the accelerometer output for 60_ and 60 changes. ISU indexing 
precision will also affect the accuracy with which the gyro nonorthoganal- 
ily terms can be estimated. 
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10. 3 ENVIRONMENTAL ERROR SOURCES 


The effect of test stand base motion is to cause extraneous rates 
and accelerations to be sensed by the gyros and accelerometers. The 
analysis is patterned after that in the UNIVAC study as pointed out in the 
ground rules (Section 2. 3), and the significant error sources are: 

A0 n = change in 60 ^ and 60^ (the derivative of this (A0^) is also 
considered an error source) 

a = vibration acceleration, assumed to be isotropic in all 
directions 

These errors are specified in terms of power spectral densities 
(PSDs) as shown in Figs. 10-1 and 10-2*. The PSDs are approximated 
for computational purposes as shown in Tables 10-1 and 10-2**. It is 
noted in Table 10-1 that there are several apparent inconsistencies 
between Table 10-1 and Fig. 10-1. Since a ground rule of the analysis was 
to use the results of the UNIVAC study (Ref. 4), this was done as shown in 
Appendices E and F. Generally, a conservative approach was taken. The 
discrepancies should be resolved, however, since the error source is 
significant (for both the gyros and accelerometers) and the expected calibra- 
tion precision is directly influenced. 

Since the effect of A0^ on the accelerometer output is particularly 
critical, as discussed in Section D. 1 of Appendix D, the bubble level com- 
pensation procedure presented in Sections 8. 2 and 8. 3 is specified. The 
objective is to attenuate the low frequency portion of the A0^ spectrum. In 

*These PSD's are as listed in Ref. (4), Figs. 3-3 and 3-4, which in turn 
were reproduced from those presented in Ref. (8), Fig. 1-1. 

**The approximations are as listed in Ref. (4), Tables A-3 and A-2. 
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accordance with the ground rules of Section 2. 3, the UNIVAC analyses was 
used. On page 2-17 of Ref. 4, the following two models were proposed to 
characterize the use of bubble level compensation: 

(1) The rotational noise spectrum, PSD (A 0 ), was reduced to zero 
below a frequency corresponding to a 50 minute period. 

(2) The rotational noise spectrum, PSD (A0 ), was assumed to be 

the squared modulus of a first order transfer function having 

an RMS noise in'A0 of 4. 5 sec and a half-power frequency of 
-2 n 

10 Hz, corresponding to a 1. 7 minute period. 


The results on the accelerometer output using these assumptions are 
plotted in Fig. E-5, as reproduced from Fig. 2-5 on Ref. 4. Model #1 
yields unsatisfactory performance and would imply that compensation is 
required more often than every 50 minutes. Model § 2 is more compatible 
with the calibration precision goals and is used in the error analysis. A 
preliminary analysis has indicated that the true model may be one that 
slopes to zero as the frequency decreases from that corresponding to the 
averaging time. Because of the critical nature of the error source, it is 
recommended that further analyses be conducted. 
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Figure 10-1 Representative Angular Vibration and 
Base Motion of Test Stand 
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Figure 10-2 Representative Linear Vibration of Test Stand 
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Table 10-1 Approximation of Test Stand Angular Vibration PSD 



PSD (A0 n ) 

f (Hz) 

,2 . 

rad /Hz 

sec /Hz 

— R 

0. to 5 x 10 

9. 42 x 10" 4 

4. 0 x 10 7 

5x10 to. 317 

2. 345 x 10~ 4 f” 2 

9. 95 x 10V 2 

. 317 to 3. 17 

2. 345 x 10" 13 

IQ" 2 

3. 17 to 15. 7 

-12 -2 

2.345 x 10 f 

-2 -2 

9. 95 x 10 f 

above 15. 7 

1. 022 x 10“Y 5 

2 -5 

4. 34 x 10 f 


Notes: 

(1) This table was derived from Table A-2 of Ref. (4) by removing 

the earth rate coefficient of PSD (A0 ). 

n 

(2) The second and third lines of the table are as listed in Ref. (4) 
even though the values do not correspond to those in Fig. 10-1, 

This is discussed in the text (Section 10. 3) 


Table 10-2 Approximation of Test Stand Linear Vibration PSD 


f (Hz) 

PSD(a ), g 2 / Hz 

V 

0 to 10” 2 

lo’ 15 

- 2 -1 
10 to 10 

to" 11 

' lo." 1 to 1 

o 

CD 

1 to lo 2 * 

10- 8 

10 2 to 10 3 

to- 11 

3 

above 10 

0 


# 3 

This was listed as 10 in Table A -3 of Ref. (4). 
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10. 4 ACCELEROMETER ERROR ANALYSIS 


10. 4. 1 Error Sensitivities 

The accelerometer error sources were identified as described in 
Section 10. 1. The detailed error sources considered are tabulated in 
Table 10-3, including lcr values of the magnitudes of the errors. It is 
noted in the table that pre -calibration alignment error sources are to be 
determined, as discussed in the previous sections (10. 2 and 10. 3). 

The effect of these error sources on the accelerometer output 
(i. e. , components of &M m ) are considered in detail in Appendix E for each 
test condition (m). Each component of 6M is related to the error sources 
of Table 10-3 as shown in Tables E-l and E-2. The actual propogation of 
the errors is shown in Table E-l, for accelerometer #1, and Table E-2 is 
a summary of the numerical calculations used to determine the accele- 
rometer error covariance matrix elements for each test condition (m). 

Appendix E also contains plots of the effect on 5 M of each potentially 
significant error source (see Figs. E-l thru E-6). Use of these plots is 
discussed in Appendix E. Fig. E-l graphically shows the estimation pre- 
cision goals stated in Table 4-3 and provides a basis for determining 
acceptable levels of the various error sources, as described in Appendix 
E. The critical error source is shown in Table E-2 to be the cross-coupl- 
ing of gravity due to rotation of the test stand base (A0^). The errors listed 
are for the case when bubble levels are used. Detailed discussions of 
errors associated with bubble level compensation are contained in Sections 
10. 2, 10. 3 and D. 1. Fig. E-5 shows the final effect on 6M, as a function 
of run time (At). 
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TABUS 10-3 ACCELEROMETER ERROR SOURCES CONSIDERED 


ERROR. NO. 

ERROR SOURCE 

SYMBOL 

1 a VALUE 

1 

Compensation for C^q 

§C 

IO 

(1) 

2 

Compensation for C^p 

6c 

IP 

(1) 

3 

Compensation for magnitude of 
g(effect of error is correlated) 

6g 

-3jug 

4 

Compensation for ISU cube to 
bubble level misalignment 
(for calculating T BA) 

Compensation for (6 6 0 ^ } 

6(B. L) 

(2) 

5 


(2) 

6 

Error in setting bubble level 




to null 

(L- V) 

1 sec 

7 

TAF positioning error 

6 0 

3 

(2) 

8 

TAF 0^ positioning error 

6 0 

4 

(2) 

9 

Nonorthogonality of 0 and 0 . 
axes ^ 

M03' 

(2) 

10 

Nonorthogonality of 0 g and0^ axes 

MR3 1 

(2) 

11 

Misalignment of IA to body 
axes (cube) 

Mli 3) 

5/3 min. 

12 

Misalignment of OA to body 

A 



axes (cube) 

(O -B) 

10/3 mm. 

13 

Misalignment of PA to body 




axes (cube) 

(P- B) 

10/3 min. 

14 

Misalignment of body axis to 




bubble level 

(B- L) 

10. 4/ 3 min. 

15 

Effect of uncompensated Cjp 

C IP 

io«g 

16 

Effect of uncompensated C-j-q 

O 

i — i 

o 

10 w g 

17 

Sampling time error 

5(At) 

100 u sec 

18 

Quantization 

q A 

e A dt 

-n 

. 0025 a/T fps 

412 

19 

Accelerometer internal random 

negligible 


noise 

111 

20 

Accelerometer rebalance 

A 

e nreb 

included 


loop noise 

with #18 
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TABLE 10-3 
ERROR NO. 


ACCELEROMETER ERROR SOURCES CONSIDERED 
(Cont. ) 

ERROR SOURCE SYMBOL 1 cr VALUE 


21 

22 


Angular base motion of test 
stand., about level axes 

Linear vibration of test 
stand (in all directions) 



See Section 

t> 

10. 3 

a v 

See Figure 
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Negligible Error Sources 

• Qj P w p w i 

• Compliance (other than Cjq and C^p) 

• Rotational cross- coupling 

• Cross -coupling of accelerations along OA 

• Second order terms (except as listed above) 


Notes: (1) These error sources are assumed negligible for this analysis, 
pending results from laboratory tests 

(2) Values for these error sources are to be determined once the 
pre- calibration alignment scheme is specified in detail 

(3) MIA - (MOA 2 + MPA 2 ) 1/2 
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The error covariance matrix elements derived in Table E-2 were 
used in Eq. (10-2) to determine the expected precision (i. e. , standard 
error) with which each accelerometer calibration term can be estimated. 
This is considered in the next section. 


10. 4, 2 Expected Accelerometer Calibration Precision 

The error covariance matrix elements required to perform the 
error analysis using Eq. (10 -2) are derived in Table E-2 and summarized 
in Table 10-4. The H matrix required in Eq. (10-2) is given for each 
accelerometer in Tables B-2, B-3 and B-4 by selecting the appropriate 
rows according to the measurements used in the estimation equation being 
analyzed (see Tables C-l, C~2 and C-3). The results of the accelerometer 
error analysis are summarized in Table 10-4. 


Generally the estimation precision goal is met, except for 6V c ~,7 
for accelerometer #3 which is 32% oyer the goal. It should be mentioned 
that not all error sources have been included, as discussed in Sections 
10. 2 and 10. 3, pending final definition of the pre-calibration alignment 
procedure. 


The effect of 6g (equal to 3 ug) is included in the 6V c terms, taking 
into account the fact that the error is correlated between the various 

A . A . a. 

Sm. The error propogates one-to-one for 6 and and by 

• j. 7 A • _ 7 

11/2 for 6V * and 6V 
c c 

tion terms. 


c c 

The error has no effect on the other calibra- 


The precision of the Minimum and Partial Calibration processes 
are essentially the same as that shown in Table 10-4, which is for the Pull 
Accelerometer Calibration Process. In the case of Partial Calibration, 
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TABLE 10-4 EXPECTED ACCELEROMETER CALIBRATION PRECISION 
Calibration Term Precision Standard Error(SE) of Y in 


Estimated 

Units 

Goal 


Accel 

. # 





#1(2) 

Kvmn 

mmm 

K ISIS! 

A 

K 

b a 

Mg 

8. 3 

3. 96 

3. 76 

3. 98 

2. 41 

-(NOA)g 

ug 

TBD 

6. 43 

5. 35 

3. 72 

2. 91 

(NPA) g 

ug 

TBD 

'3. 59 

3.76 

6. 62 

3. 40 

6Vc +1 

Mg 

8. 3 

3. 02 

3. 02 

3. 02 

3. 02 

^Vc" 1 

A ^ 

wg 

8. 3 

3. 02 

• 3. 02 

3. 02 

3. 02 

«v c + - 7 

A C 

Mg 

8. 3 

3. 28 

3. 28 

11. 00 

2.77 

6 y c -- 7 

Mg 

8. 3 

3. 28 

6.70 

2. 77 

2.77 

9ip 

“ elet 

TBD 

10. 39 

8. 63 

6. 35 

5. 43 

i — i 

o 

Mg/g 

TBD 

7.14 

7. 02 

12. 28 

5. 98 


Square Root of Error Covariance 
Diagonal Elements 


m 

TT 

V K ss 

"~n — 

fm 

V e ss 

#3 
rrrf 
V e ss 

mail) 

v e ss 

4 

. 63 

6. 2 

6. 2 

. 63 

5 

6. 2 

6. 2 

. 63 

6. 2 

6 

6. 2 

. 63 

6. 2 

6. 2 

1 

5. 00 

- 

. 54 

5. 00 

2 

- 

. 54 

5. 00 

5. 00 

3 

. 54 

5. 00 

- 

.54 

7 

3. 56 

- 

3. 56 

3. 56 

8 

3. 56 

3. 56 

- 

3. 56 

9 

- 

3. 56 

3. 56 

5. 00 

10 

3. 56 

- 

3. 56 

3. 56 

11 

3. 56 

3. 56 

5. 00 

3. 56 

12 

5. 0 

3. 56 

3. 56 

5. 0 

13 

- 

3. 56 

- 

- 
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compensation for C and C is required. Errors in this compensation 

A • _j. ^ A . _ ^ 

will only cause additional errors inf>V" c * and &V c ‘ that are equal to 
. 5g 2 6 C IQ and . 5g 2 5 C jp 

Included in Table 10-4 is a column that indicates the expected 

calibration precision assuming all data available is utilized. This .was 

done only for accelerometer #1 to illustrate the potential improvement that 

is possible. However, it does not necessarily follow that more data (i. e. , 

measurements 3$[) will always yield better results. Since the least 

squares estimation algorithm is not optimum, it is possible that an 

additional measurement can have a relatively large error associated with 

it that can propogate through the weighting coefficients {of Appendix C) 

so as to cause certain terms to be estimated with less accuracy. A 

weighted least square algorithm would provide weighting coefficients that 

could avoid this problem. It is possible that the relatively large error 

+.7 

associated with 0V * for accelerometer #3 could be reduced in this 

c 

manner. A study of this approach is recommended. 

a 

The precision with which T can be estimated is dependent upon 

ML A A 

the measurement accuracy of T , as well as SE(NOA) and SE(NPA) as 
indicated by Eqs. (9-2 5) thru (9-30), The error propogation is one-to-one. 
From Table 10-4, 

— A A _ 

.7 sec SE{NOA, NPA}< 1.4 sec (10-5) 

— ML 

which is within the precision goal of 1. 7 sec. To meet this goal, SE(T ) 

must be <1.0 sec. 
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10. 5 GYRO ERROR ANALYSIS 


10. 5. 1 Error Sensitivities 

The gyro error sources were identified as described in Section 
10.1. The detailed error sources considered are tabulated in Table 10-5, 
including la values of the magnitudes of the errors. It is noted in the table 
that pre- calibration alignment error sources are to be determined, as dis- 
cussed in the previous sections (10. 2 and 10. 3). 

The effect of these error sources on the gyro output (i. e. , com- 
ponents of ^M m ) are considered in detailed in Appendix F for each test 
condition (m). Each component of 6M m is related to the error sources of 
Table 10-5 as shown in Tables F-l and F-2. The actual propagation of the 
errors is shown in Table F-l, for gyro #1, and Table F-2 is a summary of 
the numerical calculations used to determine the gyro error covariance 
matrix elements for each test condition (m). 


Appendix F also contains plots of the effect on 6M of each potentially 
significant error source (see Figs. F-l thru F-8). Use of these plots is 
discussed in Appendix F. Figs. F-l thru F-3 graphically show the estima- 
tion precision goals stated in Table 4-5 and provides a basis for determining 
acceptable levels of the various error sources, as described in Appendix F. 
The critical error sources are shown in Table F-2 to be the angular rate of 
vibration of the test stand (A^h) and the nonuniform table rate error(NUWT). 
Very little can be done economically with the former, other than perhaps to 
provide a more massive foundation for the test stand, whereas the later 
error may be reduced by a better quality table rate drive. For the gyro 
scale factor tests, the largest error source is the test table rotary angle 
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TABLE 10-5 GYRO ERROR SOURCES CONSIDERED 


ERROR NO. 

ERROR SOURCE 

SYMBOL 

1 cr VALUE 

1 

Compensation for C^ 

6 c 

II 

(1) 

2 

Compensation for Cgg 

6 c 

SS 

(1) 

3 

Compensation for C^q 

6 c 
TEO 

(1) 

4 

Compensation for C^g 

6 c 
^IS 

(1) 

5 

Compensation for ISU 
cube to bubble level mis- 
alignment (for calculating 
^BG) 

6<B . L) 

(2) 

6 

Compensation for earth 
rate coupling through 

t gm t ml 

S(G. L) 

< .5/ 3 min. 

7 

Bubble level coordinate 
frame to table axis mis-' 
alignment. 

(L • ft') 

(2) 

8 

Nonorthogonality of 0„ and 
0^ axes 

M03 ' 

(2) 

9 

Nonorthogonality of 0„ and 
02 axes 

MR 3 ’ 

(2) 

10 

TAP 0g positioning error 

60 

3 

(2) 

11 

TAF 0£ positioning error 

60 4 

TBD(2) 

12 

Scale factor error at earth 
rate 

«Q e 

200/ 3 ppm 

13 

Azimuth misalignment of 
bubble level coordinate 
frame to East 

(L ♦ N) 

20 sec 

14 

Misalignment of gyro I A 
to body axes(cube) 

mia (3) 

5/3 min. 
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TABLE 10-5 GYRO ERROR' SOURCES CONSIDERED 

(cont. ) 


ERROR NO. 

ERROR SOURCE 

SYMBOL 

1 a VALUE 

15 

Misalignment of gyro SA to 
body axes about I A 

(S G . B) 

10/3 min 

16 

Misalignment of body axes 
to bubble level 

(B • L) 

10. 4/ 3 min. 

17 

Effect of uncompensated 

< B I> 

1. 5°/ hr/ g 

18 

Effect of uncompensated Bg 

B S 

1. 5°/hr/ g 

19 

Tilt of test table axis about 
East 

% 

. 2 mTn. 

20 

Nonuniform table rate 

NUWT 

3. 7 sec 

21 

Tilt of test table axis 
about North 

6 0 

N 

. 2 min 

' 22 

Distance of gyro along 
normal to table axis 

m 

10 inches 

23 

Sampling time error 

6< At) 

100 jusec 

24 

Uncertainty of test table 
angle read out 

5 <A0 2 ) 

6 sec 

25 

Quantization 

G 

q 

See Note (4) 

26 

Gyro random noise 

e G 
nr 

See Table 4-6 

27 

Gyro rebalance loop noise 

e G 

nreb 

Included in #25 

28 

Angular vibration of test 
stand base, about level 

A0 

n 

See Fig. 10-1 


axis 



29 

Angular test stand rate 
about level axis 

M 

n 

See Fig. D-17 

30 

Error in setting table 
rate 

5(w T ) 

. 5% 
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TABLE 10-5 GYRO ERROR SOURCES CONSIDERED 

_ (cont. ) 


Negligible Error Sources 

® Linear environmental vibration 
® Error in setting bubble level to null 
• Second order terms (except as listed above) 


Notes: (1) These error sources are assumed negligible for this analysis, 

pending results from laboratory tests 

(2) Values for these error sources are to be determined once 
the pre- calibration alignment scheme is specified in detail 

(3) MIA = (MOA 2 + MSA 2 ) 1/2 

(4) = , 844 sec for |wj < 30°sec. and 1. 688 sec for 
30 < i w. i < 60°/sec. 

1 i 
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output error (5 (A$2^ with misalignment of the gyro IA parallel to the 
rotary axis as the next largest. Both of these- errors could be reduced 
without much trouble, if desired. Use of the bubble levels has little 
effect on the gyro calibration precision (other than in measuring T^ ). 

The error covariance matrix elements derived in Table F-2 were 
used in Eq. (10-2) to determine the expected precision (i. e. , standard 
error) with which each gyro calibration term can be estimated. This is 
considered in the next section. 

10. 5. 2 Ejected Gyro. Calibration Precision 

The error covariance matrix elements required to perform the 
error analysis using'Eq. (10-2) are derived in Table F-2 and summarized in 
Table 10 -6. The H matrix required in Eq. (10-2) is given for each gyro 
in Tables B r -5, B-6 and B-7 by selecting the appropriate rows according 
to the measurements used in the estimation equation being analyzed (see 
Tables C-r7, C-8 and C-9). The results of the gyro error analysis are 
summarized in Table 10-6. 

Generally, the estimation precision goal is essentially met*. for 
the bias and g sensitive terms, except for gyro #2 which is about 50% 
oyer the goal. For the major compliance terms, the errors are from 30% 
to 62%over but probably are only significant for the term of the Z gyro 
(#1), since the IAs and SAs are all nominally level for the other two gyros 
when in flight. If all the measurements (i. e. , 2 extra samples) are used. 


within approximately 20% of the goal 
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TABLE 10-6 EXPECTED GYRO CALIBRATION PRECISION 

'/i ’ 

Calibration Term Precision Standard Err or(SE) of Y in 


Estimated 

Units 

Goal 


Gyro # 


A 



#l(z) 

fflx) 

■my) 

t 1 (all data) 

R b 

mdh 

7.0 

8. 52 

7.42 

7.77 

7. 03 

A 

B I 

mdh/ g 

7. 0 

6. 30 

10. 33 

6. 02 

2. 78 

A 

B s 

mdh/g 

7.0 

8. 54 

7.45 

7. 29 

8. 54 

A 

B o ■ 

mdh/g 

7.0 

8. 54 

7. 74 

6. 38 

6.96 

A 

mdh/ g 2 

8. 3 

13. 46 

13.93 

10.90 

8.15 . 

c ss 

mdh/ g 2 

8. 3 

11. 66 

10.75 

11. 00 

10.83 

& 

mdh/ g 2 

8. 3 

17. 26 

19. 07 

16. 05 

1-4. 90 

C IO 

A 

mdh/ g 2 

8. 3 

17. 23 - 

15. 21 

14.98 

12. 93 

c os 

mdh/ g 2 

8. 3 

20, 45 

29. 34 , 

20. 82 

. 17. 09 


m 

Square Root of Error Covariance 
Diagonal Elements 

#1 

v ss 

#2 
. / m 

V e ss 

#3 • 
V ss 

#l(all) 
■ V s s 

4 

2. 8 

11.4 

11.4 

2. 8 

5 

11. 4 

11.4 

2. 8 

11. 4 

6 

' 11.4' 

2. 8 

11. 4 

11. 4 

1 

12.4 

1,2. 4 

- - 

12.4 

2 

12. 4 

- 

12. 4 

12.4 

3 ' 

- 

12.4 

12. 4 

2.7 

7 

8. 9 

- 

8. 9 

8. 9 

8 

8. 9 

8. 9 

- 

8. 9 

9 

- 

8. 9 

8. 9 

12. 3 

10 

8. 9 ■ 

12. 3 

8. 9 

8. 9 

11 

8. 9 

8. 9 , 

12. 3 

8. 9 

12 

12. 3 

8. 9 

8. 9 

12. 3 

13 

- 

8. 9 

12. 3 

■ 
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the goal can be met, as shown in Table 10-6. In the case of the cross 
compliance terms, the"errors' are from 80 to 150% over and in the case of 
C oc for the #2 gyro (x), 2 50% over. The estimation precision goals can- 
probably be relaxed for the cross compliance terms since gravity effects 
will yield only second order errors under normal flight conditions, and in 
-the worst case" when the'ISU is tilted 45° from the nominally level positioi 
the resulting rate errors propagates at half the values shown in Table 10-t 
Although not all error sources have been included (viz. , those associated 
with the pre- calibration alignment procedure), they are not anticipated to 


alter the results of Table 10-6 significantly. They are important, however, 


A 

- in tRe estimation precision of T 


BG 

9 


which is discussed below. 


The precision of the Minimum and Partial Calibration processes 
are approximately the same as that shown in Table 10-4, which is for the 
Full Gyro Calibration process. In the case of Minimum Calibration, 


compensation for C_ and C is required. Errors -in this compensation 

Ur bb 

will only cause additional errors in and that are equal to . 56Cgg and 

(. 5§Cgg+6c^), respectively. The total error should be about the same 

as that shown in Table 10-6. The B 0 and B terms .are not affected at all. 

d o 


In the case of Partial Calibration, compensation for C and C is 

IxJ Id 

required. Errors in this compensation propagate at 40% into the estimation 

A. A. 

of the Bj and C terms, and is expected to cause a slight degradation in 


estimation precision compared to that of a Full Calibration. The other 


terms are unaffected. 


Included in Table 10-6 is a column that indicates the expected 
calibration precision assuming all data available is utilized. ^This was 
done only for gyro #1 to illustrate the potential improvement that is 
possible. However, it does not necessarily follow that more data 
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(i. e. , measurements M) will always yield better results. Since the least 
squares estimation algorithm is not optimum., it is possible that an addi- 
tional measurement can have a relatively large error associated with it 
that can propagate through the weighting coefficients (of Appendix C) so 
as to cause certain terms to be estimated with less accuracy. A weighted 
least square algorithm would provide weighting coefficients that could avoid 
this problem. It is possible that the relatively large errors associated with 

with B , C and C for gyro #2 and C for gyro #1 could be reduced in this 
1 11 CJo 11 

manner. A study of this approach is recommended. 


A BG 

The precision with which T can be estimated is dependent upon 

ML A A 

the measurement accuracy of T , as well as SE(NOA) and SE(NSA) as 
indicated by Eqs. (9-33) thru (9-38). The error propogation is one-to-one. 

The precision with which NOA and NSA can be estimated is determined using 

♦ 

the estimation equations in Table C-5 and the SM errors given in Table F-2 

T 

for m = 41 ; 51 and 61 andjw | = 2° /sec. The results are as follows: 


A A 

se (noa, nsa3 


11.13 x 10 
7 r [2 — ( —2 )] / 360 


= .32 sec (10-6) • 

T 

The error in knowing W has an insignificant effect in Eq. (10-6). Consider- 
ing errors in aligning the bubble levels to the test table rotary axis, and 

ML 2 2 1/2 — ' 

errors in the measurement of T , a total of (1. 7 - . 32 ) =1.6 sec 

(let ) can be allocated to these measurements and yet meet the estimation 

precision goal of 1. 7 sec (1 cr ). 

In the case of the gyro scale factor error estimate, the errors pro- 
pagate one-to-one as shown in Table C-5. From Table F-2, for m = 72, 

82, and 92, 
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(10-7) 


SE&J C ) = 6M = 18. 6 mdh for W T = + 2°/sec 

for any of the three gyros. This is within the precision goal of 62 mdh, as 
shown in Fig. F-2. 
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11. CONCLUSIONS 


The following conclusions are drawn from the studies made: 


1. The calibration system and procedure developed in 
this study can be used in the laboratory as well as 
under field conditions. 

2. The calibration can be performed within an 8 hour 
period, providing a minimum number of gyro scale 
factor rates are used. 

3. Accelerometer calibration precision is generally within 
the goals specified. 

4. Gyro calibration precision is generally equal to the 
goals specified, except the B term of gyro #2 (which 
is 13 mdh/g, rather than 7 mdh/g). 

5. Estimation precision of the compliance terms is worse 
than that for the basic terms, but probably of acceptable 
level for the accelerometers and marginal to probably 
acceptable for the gyros. 

6 . The additional runs required to estimate the cross 
compliance terms generally have only a minor affect 
on the precision with which the other terms can be 
estimated. 

7. Errors in the knowledge of | g| affect only the accelerometer 
scale factor terms. 

8 . The data sampling error ( 6 (At)) is critical for the gyros, 
but not the accelerometers. 

9. The angular vibration (A 0 n ) is critical for the accelerom- 
eters and the angular vibration rate (A 0 ) is critic ' 

the gyros. 
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10. Using data from additional runs to improve precision of 
estimates does not always yield the desired results, 
depending upon the errors associated with the additional 
measurements. 

11. The hardware requirements to support the total calibration 
process are considered to be feasible technically and 
economically. The alignment accuracy of the ISU and 
TAF gimbals before each run, however, are yet to 

be specified in final detail. 
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12. RECOMMENDATIONS 


It is recommended that development of the calibration procedure 
continue along the lines currently being followed since most objectives 
are basically being met. However, the design can very likely be improved, 
if deemed necessary, by considering a different combination of test con- 
ditions. The primary results of this would be to reduce test time and 
possibly improve estimation precision. Alternatively, estimation pre- 
cision requirements should be reviewed to determine if any can be relaxed, 
particularly for the gyro cross compliance terms. 

Another way to provide significant improvements in the design is 
to consider using weighted least squares, particularly for estimation of 
the scale factor error coefficients. This is recommended since instances 
of reduced calibration precision for certain terms were experienced in 
this study upon using additional data, due to the relative errors of the 
various measurements. This effect would be eliminated using weighted 
least squares estimation algorithms, and may provide significant im- 
provements in calibration precision and/or reduce test time with little 
effect on computational complexity. 

It is also recommended that the resolution of the test table pulse 
output(for A0 ) be improved from 6 sec (la) to approximately 1 to 2 sec, 
along with reducing the sampling time error from 100 jusec to about 10 p sec. 
If the nonuniformity of table rate (NUWT) can be held better than 4. 3 sec/rev 
(as defined in Eq. (6-1)), gyro calibration precision can be further improved 
sig nif icantly. Reduction much less than ~2 sec/rev is not warranted, how- 
ever, assuming the other errors remain as shown in this report. 
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It is believed that the pre- calibration alignment procedure proposed 
is basically sound and is an improvement over previously used methods. 

It is recommended that development of the method be continued along the 
lines proposed, and that the associated procedures and error analysis be 
updated accordingly. Similarly, final details of the data collection system 
require definition. 

The “'base-line" design proposed herein should be reviewed to de- 
termine what changes are desired, if any, to provide a more "balanced" 
design (i. e. , tradeoffs between the various equipment performance 
requirements, calibration run times, alternate procedures, etc. ). Not 
only is the error analysis useful in doing this, but consideration should be 
given to the revisions associated with use of the flight computer to collect 
and reduce the data, assuming a final decision is made along these lines. 

It is recommended that the error analysis be revised to reflect 
more final consideration of the various error sources. Also analysis of 
the final precision in estimating T-^ and T^A should be performed, as 
well as that of the scale factor error coefficients. This would include 
analysis of errors to be identified once the final pre -calibration alignment 
procedures are finalized. The effect of errors in compensation for com- 
pliance terms in the Minimum and Partial Calibration procedures should 
also be refined. 

It is strongly recommended that the apparent inconsistencies of the 
environmental models and associated error analysis in the Uni vac study be 
■resolved and the analyses in this report corrected accordingly. Not only is 
❖ 

Such as compensation for 60^ and 60 , gimbal nonorthogonalities, 
gimbal angle readout nonlinearities, etc. 
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the environmental angular motion particularly critical, but an initial inves- 
tigation has indicated a better analytical approach that will more realisticall • 
determine the effect of bubble level compensation. 

Finally, it is recommended that the centrifugal acceleration analysis 
be refined to reflect a more precise value of the moment arm (m) of the 
gyros from the SATS rotary axis. This effect on the gyro scale factor 
error coefficients is potentially significant at the higher table rates. 
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APPENDIX A 


DEVELOPMENT OF ACCELEROMETER MATH MODEL 

The purpose of this appendix is to document the salient errors 
associated with a pendulum type of accelerometer under the condition 
of vigorous linear and angular case motions. Detailed derivations are 
not presented., rather the error sources are stated and categorized. 

The performance equation for a pendulous accelerometer 
(Fig. A - l) can be evolved by writing a moment equation about the 
accelerometer output axis, thus, 

1(0 +w)=-D0 - K0 + T + £ T . (A-l) 

o o o o o u OA 

where I , D, and K are the moment of inertia, damping constant and 
spring constant, respectively and 0^ is the pickoff angle on the output 
axis. The spring constant in pulse-rebalanced instruments is the result 
of torquing the float so as to null the float error angle. T denotes the 
uncertainty torques, while £ T represents the sum of all remaining 

xJ-LX 

kinematic and gravity related torques. This last term can be separated 
as follows: 


E T 


OA 


£ T 

. linear 


+ £ 


Using small angle approximations, the component of output-axis torque 
due to linear accelerations in the absence of mass unbalance and gimbal 
deformations is 


A-l 




(A-3) 


mZ (a. + 0a ) 
1 p 


Torques arising from mass unbalance and (kinematically produced) gir 
Reflections are expressed as linear and quadratic terms in the three c. 
accelerations. Here we lump them as one compliance term. 


M (a., a , a ) 
c i o p 


(A-4) 


and breakout the compliance of the pendulous axis as the following separate 
term: 


-0 a = (C T _. a.)a 
p o IO i o 


(A- 5) 


The dominant gyroscopic torque is 


(I ~ I.) w w. 
P i p i 


while that due to rotational acceleration is 


(A- 6) 


2 a/ 2 2 2 . 

mZ 0(w - w. - w ) 
p i o 


(A-7) 


Using these relationships Eq. (A-l) can be rewritten as 


10 + D0 + K0 = m4(a. + 0 a ) - I w + (I - I.) w w. + C T ,-.a.a 
oo o o l op oo pipi IO 1 o 


+ mj& 2 0(w^w/Y/^) + T (a., a , a ) + T 
pio ciop u 


(A-8) 
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We now define the equation of motion of an ideal accelerometer as 


10 + D0 T + K0 T = mi a. 
o I I I 1 


(A-9) 


Denoting the difference in angular response of these two instruments as 


e c = 0 - 
. 6 o I 


(A-10) 


yields 


I e + De + Ke = mi 0a - I w + (I - I.)w w. + C _ a. a 
00 0 o poopipi IO 10 


+ mi^0(w 2 ' w?- w^) + T (a., a ,a ) + T 
pio ciop u 


(A-ll) 


The corresponding steady- state error is 
j mi a. 

e = — (mi ( — ~-) a - I w +(I - I.)w w. + C,„a.a 
0 K K p o o p i p l IOio 


2 2 2 2 

+ mi ( — — ) (w ~ w. - w } + T (a. s a , a ) + T ) 
K p i o c rop u 


(A- 12) 


where we have used the relationship 


* o mi 
0 « 0 T f» ~ ~ a. 
o I K l 


(A-13) 
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Converting the angle error to an acceleration error by multiplying by 
K/mX yields 


e 

a 


"f 


m $, 
K 

I 


a. a 
1 p 


o • 

— » w 
mX o 


(I - I.) 
a 

mX 


w w. 
p x 


(cross- coupling or vibropendulous error) 


(OA rotation error) 


(anisoinertia or dynamic unbalance error) 


+ C _a.a 
IO 1 o 


, mi ,2 2 2, 

+— — — a.(w - w.- w ) 
K x p i o 


T (a , a , a ) 
+ -...C-J- o P 
mX 


+ 


T 

u 

mX 


, ,, mX_ 

+ (6 — ) a. 


(compliance of pendulous axis) 


(rotational cross-coupling error) 


(compliance error - not including Cjq) 


(uncertainty error, including bias term) 


(scale factor error) 


(A-14) 


The last term displays the dominant effect of an error in knowledge of the 
instrument scale factor, mX/K. 



APPENDIX B 

TABULATION OP CALIBRATION EQUATIONS 

In Section 7 3 it was shown that both the accelerometer and gyro 
.c~i.li! ration equations could be arranged in the following form: 

M^At* 1 s M m = H m * Y. + 6M m (7-1 and B-l) 

j 3 ~2 ”] J 

in which the various terms are defined by Eqs. (7-7) (7~8) and (7-9) for 
the accelerometers and Eqs. (7-25) through (7-29) for the gyros. The 
purpose of this appendix is to apply these equations to each inertial sensor 
for each test condition (m) specified in the calibration procedure (Section 
8). These equations are then used as a basis for deriving the calibration 
term estimation equations in Section 9, as well as the error analysis in 
Section 10, 

The first step is to determine the relationships (i. e. , transforma- 
tion matrices) between each inertial sensor coordinate frame and the 

T 

"environment" (i. e. * w and w ). The transfonnation matrices are 

e 

defined in terms of the unit vectors that define the coordinate frames, as 
follows: 


- 


- 


- 

Si 


% 


u 

-2 

= t BE 

^ 2 - 

= t BE 

E 



3 _ 


N 


B-l 



th 

where a multiplication of the k row of the 3x3 transformation matrix 
BE 

T with the E^ column represents B^ expressed in the E^ frame. 
For example, a vector known in the E^ frame can be expressed 
in the B frame by computing the dot product of both sides of Eq. (B~2), 
as follows: 



and 

w ‘B = component of w along B = w sin L(B,‘ U) +0 (B • E) + 

— e — 1 — e & — 1 e —1 — — 1 — 

w cos L(B *N) (B-4) 

e — 1 — 

where (B^ TJ), (B^ E) and (Bp N) are the first row elements of T 
The components of w along B and B are computed similarly. In 

6 A o ijt 

this appendix, the components of _g, and w along the inertial 
sensor axes are determined for each test condition (m), using the above 
approach. 


The test conditions are defined in terms of the TAF gimbal 
angles (0 and 0 ), the SATS rotary angle (0 ), and the test table rate 

rp O 4 2 a 

(w ). In addition, the following four angles are included in the trans- 
formation matrices in anticipation of their possible use in the estimation 
equations : 
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60 ^ vertical tilt error of the rotary axis about North 

80 vertical tilt error of the rotary axis about East 

J i 

MR3' Nonorthogonality between the TAF outer gimbal 

axis and the test table rotary axis, measured 

about the R ' axis 
o 

M03 1 nonorthogonality between the TAF inner and outer 

gimbal axes, measured about the O' axis 

O 


The transformation matrices relating the "environment" and the ISU 
body axes are as follows: 

: 

rpBRm _ notary axes to body axes transformation matrix 


.pBEm _ earth axes to body axes transformation matrix 


,pBEm „ T BRm T R'Em 


(B-5) 
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(B-6) 


T BRm = T BM T ML T LI T IOM T 00' T 0’Rm T RR' 

n 



T R'Tm t TF t FE 


B-4 



and the various transformation matrices are defined in Table B-l and the 
coordinate frames in Tables 5-1 and 5-2. Fig. 5~3 contains an illustration 
of the various coordinate frames. The elements of the matrix 

represent the 1SU optical cube to bubble level misalignments that are 
measured when the ISU is mounted within the TAF. In the derivation 
here, these elements are set equal to zero since the misalignment terms 
in the Y vectors are defined relative to the bubble level coordinate frame 
(L). The way in which these terms are introduced is shown in Section 9. 


The nominal transformations from the body axes to the inertial 
sensor coordinate frames are based on the axes orientations illustrated 
in Fig. 5-2, and are expressed in matrix form, as follows: 


“a 


-l 

0 

o' 



~G 


~1 

0 

o” 

o A 

= 

0 

0 

1 

B 

(B-8) 

o G 

= 

0 

0 

1 

p_ 


0 

1 

0 



s G 


0 

-1 

0 

- — 

i - 

- 


- 

n 


- 

1 

- 


— 


B 


n 


A " 


0 

1 

o ' 



G 


~0 

1 

0 ~ 

o A 

= 

1 

0 

0 

B 

(B-9) 

lo 

Q 

= 

-1 

0 

0 

p 

— — 

2 

0 

0 

-1 

— 

n 




0 

0 

1 

A 


0 

0 

1 



G 


0 

0 

1 

o A 

= 

-1 

0 

0 

B 

(B-10) 

o G 

= 

-1 

0 

0 

p 

- -J 

3 

0 

-1 

0 

n 



3 

0 

-1 

0 


(B-ll) 


(B-12) 


(B-13) 
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Matrix 


T 

T 

T 


FE 

TF 

R’T 


t HR' 

T°'R 




n 


rpML 


t bm 


Description of Angle(s) that Characterizes Transformation 

misalignment of test table rotary axis about North (60^.) 

misalignment of test table rotary axis about East (^0^) 

rotation of test table about rotary axis ( 0 g); zero 
value corresponds to TAF outer gimbal axis (0 axis) 
pointing East 

nonorthogonality between rotary axis and 0 axis (MR3') 

o 


rotation of TAF outer gimbal about 0 3 axis ( 03 ); zero 
value is when TAF inner gimbal axis (0^ axis) and 
rotary axis are parallel to each other 

nonorthogonality between TAF 0 and 0 axes (M03 1 ) 

rotation of TAF inner gimbal about 0^ axis ( 04 ); zero 
value corresponds to the bubble level aligned 
parallel to the 0 axis. 

o 

nominal relationship between bubble level and inner 
gimbal coordinate frames - misalignments assumed 
to be negligible 

misalignment of ISU cube (mirror) coordinate frame 

relative to bubble level coordinate frame; elements 

of the form M a ^* L c = nonorthogonality of M about 

L^with respect to measured about the positive L 

axis that is normal to L, and L 

— b — c 

definition of body axes in terms of ISU cube (mirror) 
coordinate frame 


Table B-l Transformation Matrices Relating Body Axes to Earth ■ 

Coordinates 


B -6 


At this point all the transformations required to relate the "environment" 
coordinate frames to the inertial sensors axes has been defined. 


The next step is to define the "environment" as vectors in an 
appropriate coordinate frame, as follows: 


K 


w 

~e 


T 

w 


- - 

g 

o 

o 

- 

w sinL 
e 

0 

w cos L 
e 

t" 

w 

0 

0 


in the earth coordinate frame (E) 


in the earth coordinate frame (E) 


in the rotary axes coordinate frame (R 1 ) 


(B-14) 


(B-15) 


(B-16) 


The above "environment vectors ", transformation matrices, and inertial 

sensor calibration equations were combined as indicated by Eqs. (B~3) 

and (B-4), using a digital computer. The elements of the H^^and H.^ m 

J J 

row matrices were determined for each inertial sensor, for each of 
the test conditions (m) (assuming 6 m = 0). The results were organized 
in matrix form and are summarized in Tables B~2 through B~4 for the 
accelerometers and Tables B~5 through B-7 for the gyros, in which 
products of small quantities have been neglected. Errors that are not 
compensated for in the computation of M, or that are not negligible, are 
considered in the error analysis (Section 10) as part of the 6M term. 
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In the case of the gyros, the runs for which the test table was 

rotating are tabulated separately, at the bottom, of each table. These 

runs are combined with each other and other runs, as shown at the top 

of each table. This is done to separate out certain calibration terms and 

T 

simplify the estimation calculations. When other table rates (w ) are 

used, as indicated in Table 8-6, the same six equations at the bottom of 

each table are used, substituting the new values of m for those shown. 

T 

Rates other than w = + 2°/ sec are used only in the estimation of the gyro 

A A A 

scale factor error coefficients (R^, and G^), as discussed in Section 
9 and Appendix C. 


T+ T~ 

Although the w and w values for a given test condition m are 

T 

nominally equal and opposite, the errors in setting w are sufficiently 

T 

large to warrant using the actual (i. e. , measured) values of w , as 
indicated in Tables 33~5, B-6, and B~7. The compensation terms that 
are functions of NOA and NSA should be computed after NO A and NSA 
have been estimated. The terms can be estimated from T^^ T^^, 
providing T*^' is assumed known within approximately 15 to 30 sec (lcr). 
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Table B-2 Accelerometer No. 1 Calibration Equations 


A$(P A /At) 4 +g 


♦ 4 
M 


0 0 0 

0 

1 

0 

0 

0 

0 


A$(P A /At) 5 


M 5 


1 0 1 

0 

0 

0 

0 

0 

0 


* AJ 1 

A$(P A /At) 6 


M 6 


1+10 

0 

0 

0 

0 

0 

0 


b 

-4. 85 NOA 

A$(P A /At) X + 0 


M 1 


1 0 -i 

0 

0 

0 

0 

0 

0 


4. 85 NPA 

A$(P A /At> 2 


M 2 


1-10 

0 

0 

0 

0 

0 

0 


6V +1 

A<£(p A /At) 3 ~ g 

= 

M 3 

= 

0 0 0 

1 

0 

0 

0 

0 

0 


c 

ay" 1 

A4>(P A /At) ? - g/V2 


M 7 


0 0.7 

0 

0 

1 

0 

0 

+.5 


C 

• + 7 
6 V 

A$<P A /At) 8 + g/f2 


M 8 ■ 


0 +.7 0 

0 

0 

+1 

0 

-.5 

0 


C 

» - 7 

5 V 1 1 

A$(P A /At) 9 


M 9 


1 -.7 -.7 

0 

0 

0 

0 

0 

0 


c 

_ 2 

A$(P A /At) 10 - gl\f2 


M 10 


0 0 -.7 

0 

0 

1 

0 

0 

-.5 


c ip* 

2 

A^>(p A /At) n +gA(T 


M U 


0 -.7 0 

0 

0 

0 

1 

.5 

0 


_ C I° g 

A«S>(P A /At) 12 


M 12 


1 +.7 -.7 

0 

0 

0 

0 

0 

0 


A$(P A /At) l3 -gA/2 


_m 13 _ 


1 +.7 0 

0 

0 

0 

1 

.5 

0 


2 = 

1 

j 

= 1 








= 1 



Notes: (1) g = magnitude of gravity at location of SATS (to the 4th decimal place) 

(2) Units of M are jug 

(3) Matrix elements shown as . 7 are actually equal to 1 U _ 

(4) Units of NOA and NPA are sel 



Table B- 3 Accelerometer No. 2 Calibration Equations 


A$(P A /At) 4 



* 4 

M 
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0 

+1 

0 

0 

0 

0 

0 

0 


A$(P A /At) 5 



M 5 


l 

-1 

0 

0 

0 

0 

0 

0 

0 


A$(P A /At) 6 

-g 


M 6 


0 

0 

0 

1 

0 

0 

0 

0 

0 


A b 

A<!?(P A / At) 1 

+ 0 


M 1 


1 

+1 

0 

0 

0 

0 

0 

0 

0 


-4. 35 NOA 

A$(P A /At) 2 

+ g 


M 2 


0 

0 

0 

0 

+x 

0 

0 

0 

0 


4. 85 NPA 

A$(P A /At) 3 

*4" 0 


M 3 


1 

0 

-1 

0 

0 

0 

0 

0 

0 


SY +1 

c 

A^(P A /At) ? 

4* 0 

s 

M 7 

= 

1 

-.7 

-7 

0 

0 

0 

0 

0 

0 


o> 

<!• 
O i 

A<S>(P A / At) 8 

- g/V 2 


M 8 


0 

0 

+.7 

0 

0 

+1 

0 

0 

+.5 


» -i- 7 

6V ‘ 

c 

A$(P A /At) 9 

+ g/V2 


M 9 


0 +.7 

0 

0 

0 

0 

+1 

-.5 

0 


* - 7 
6 V 

c 

A 1 0 

A$(F /At) + 0 


M 10 


1 

+.7 

-.7 

0 

0 

0 

0 

0 

0 


_ 2 
C IP g 

A$(P A /At) n + gl{2 


M U 


0 

0 

+.7 

0 

0 

0 

+1 

0 

-.5 


c io gZ 

A^(P A /At) 12 - g//2 


M 12 


0 +.7 

0 

0 

0 

1 

0 

+.5 

0 



A^>(P A /At) 13 - g//2 


M 13 


0 

0 

-.7 

0 

0 

1 

0 

0 

-5 

_ 



3 = 2 

3 = 

2 









j 

= 2 



Notes: (X) g = magnitude of gravity at location of SATS (to the 4th decimal place) 

(2) Units of M are /^g 

(3) Matrix elements shown as . 7 are actually equal to 1//2 - 
-(4) Units of NO A and NPA are S ' SC . 



Table 33-4 Accelerometer No. 3 Calibration Equations 


A$<P A /At) 4 




10-1 

0 

0 

0 

0 

0 

0 


A$(P A /At) 5 - g 


M 5 


0 0 0 

1 

0 

0 

0 

0 

0 


A$(P A /At) 6 


M 6 


1-10 

0 

0 

0 

0 

0 

0 


A b 

A$(P A /At) 1 + g 


M 1 


0 0 0 

0 

+1 

0 

0 

0 

0 


-4.85 NOA 

A$(P A /At) 2 + 0 


M 2 


110 

0 

0 

0 

0 

0 

0 


4. 85 NPA 

A$(P A /At) 3 + 0 


M 3 


1 0 1 

0 

0 

0 

0 

0 

0 


e <+l 
6 V 

c 

A$(P A /At) 7 - g/\T2 


* 7 
■M 


0 0 .7 

0 

0 

1 

0 

0 

.5 


6V 1 
c 

A4>(P A /At) 8 + 0 

= 

M 8 

= > 

1 -.7 -.7 

0 

0 

1 

0 

0 

0 


•+. 7 
SV 

c 

A<i>( P A /At) 9 + gA/2 


M 9 


0 .7 0 

0 

0 

0 

1 

-.5 

0 


6V"’ 7 

c 

A^(P A /At) 10 + g/J 2 


<10 

M 


0 0 .7 

0 

0 

0 

1 

0 

-.5 


c ipg 2 

A$(P A /At) n + 0 


M U 


1 .7 -.7 

0 

0 

0 

0 

0 

0 

- 

_ c io g2 

A$(P A /At) 12 + g/v/2 


M 12 


o 

c~ 

o 

0 

0 

0 

1 

.5 

0 


3= 

A$(P A /At) 13 


M 13 


1 -7 .7 

0 

0 

0 

0 

0 

0 



j " 3 3 = 3 j = 3 


Notes:’ (1) g = magnitude of gravity at location of SATS (to the 4th decimal place) 

(2) Units of M are /ug 

(3) Matrix elements shown as .7 are actually equal to 1//2 

(4) Units of NOA and NPA are s'g'c 



,*,42+ *4 

M - M 


,.•,72+ 

M 











1 

1 

O 

0 

M i2 '- M 4 


.'72- 

M 





( 







----1-9 


UVt 52 - M 52 ')/2 


M 51 





V 


) 




1 0 
1 

c 52 , 

o 

(M 62 - M 62 ")/2 

•+ * 

M 6 - 1 











1C 62 ' 

t 

_°_ J _ 

A$(P G /A t) 4 - w sinL 


M 4 


i 

i 

0 

0 

1 

0 

0 

0 

0 

1 



[{M 52+ + M 52 ')-'c 52+ NSA]/2 


*50 

M 


i 

0 

0 

l 

0 

0 

0 

0 

0 

1 

! 



[(M 62+ 4 M 62 ’)-c 62+ NOA]/2 


M 60 


i 

0 

-1 

0 

0 

1 

0 

0 

0 

1 

1 



i^(P G /At) 1 + w e cosL 


M 1 


i 

0 

0 

-1 

0 

0 

0 

0 

0 

1 



A$(P^/At)^ + w cosLi 
e 


« 2 

M 


i 

0 

1 

0 

0 

1 

0 

0 

0 

1 

1 



A<HP G /At) 3 +w sinL 
e 


M 3 


i 

-1 

0 

0 

1 

0 

0 

0 

0 

t 

r\ 


A$(P G /At) 7 + w (sinL + COsD/VT 


M 7 


i 

-.7 

0 

.7 

.5 

0 

0 

-.5 

0 

i 

1 



A$(P G /At) 8 + w (cosL - sinL)/V2 
e 


* 8 


i 

.7 

-.7 

0 

.5 

.5 

-.5 

0 

0 

1 



A$(P G /At) 9 - w^cosL 


M 9 


i 

0 

.7 

-.7 

0 

.5 

0 

0 

-.5 

1 



A$(P G /At) 19 - w^CcosL- sinL) //if 


M 10 


i 

-.7 

0 

-.7 

.5 

0 

0 

.5 

0 

1 



A$(P G /At)^ 1 - w (cosL - sinL) 1^2 


M U 


i 

.7 

.7 

0 

.5 

.5 

.5 

0 

0 

1 

1 



A<i*{P G /At) 12 - w cosL 


'12 

M 


i 

0 

-.7 

-.7 

0 

.5 

0 

0 

.5 

i 

1 



A$(P G /At) 13 + w (sinL+cosL)/Vif 


,■,13 

M 


i 

-.7 

-.7 

0 

,5 

.5 

.5 

0 

0 

1 


— 

3 = 1 

3 = 

i 












3=1 

[<A$P G - 36OOS 0 P 0 )/At] 42 ^w e sinL 


~ *42+ 
M 


i 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 1 

o' 

[(A<£F G - 3600^I^)/At] 42 -W e sinL 


m 42 “ 


i 

1 

0 

0 

1 

0 

0 

0 

0 

' 0 

0 0 

1 

A${P G /At) 52+ 

= 

m 52+ 

= 

i 

0 

0 

1 

0 

0 

0 

0 

0 

1 o" 
1 

- T52t „ 
w 0 

0 

A$(P G /At) 52 " 


'52- 

M 


i 

0 

0 

1 

0 

0 

0 

0 

0 

1 

1 0 

T52- „ 

w 0 

0 

A$(P G /At) 62+ 


• 62+ 
M 


i 

0 

-1 

0 

0 

1 

0 

0 

0 

T624- A n 

| w 0 0 

0 

_ A«(P G /At) 62 “ 


*62- 
M _ 


_i 

0 

-1 

0 

0 

1 

0 

0 

0 

1 T62- 

r 

0 0 

0_ 


3=1 3=1 


V 

M 


B r;g 

9 

G n g ' 

, 2 
\sS g 
‘> 

c ts g 


JO" 


c os g 

NO A 


NSA 

6ii + 

c 

&y 

c 



3=1 


Notes: 


(1) w T = S 0 P 0 /At (deg/sec) 

(2) c m+ = »(w T+ + w T ") m /360 (deg/hr sec) 


(3) c m s ff (w T+ - w T ) m / 3 6 0 {deg/hr sec) 

(4) Matrix elements shown as .7 are actually equal to 1 //T 
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•62+ *6 
M - M 


*92+ 

M 

M 62 '- M 6 


3- 

CD 

tS3 

(M 42 - M 42 “)/2 


•41 

M . 

(j!i 52 - m S2 ')/2 


,•51 

M 

[(M 42+ + M 42 ") - c 42+ NSA]/2 


•40 

M 

[(M 52+ + M 52 ") - c 52+ NOA]/2 


M 50 

A$(F G /At) 8 - w sinL 
e 


M 6 

A$(P G /At) X + 0 


M 1 

A$(P G /At) 2 + sinL 


M 2 

A$(P G /At) 3 + 0 


M 3 

4$(P G /At) ? + 0 


M 7 

A$(P G /At) 8 - w e (cosL + sinL)/yj 


M 8 

A4>(P G /At) 9 + (w sinL)/^” 
6 


M 9 

A^(P G /At) 10 + 0 


• 10 
M 

A$(P G /At) X1 + iV p (cosL + sixiL)/V2" 


M 11 

A®(P G /At) 42 - w^sinL/Vlf 


M 12 

A${P G Mt) X3 + w e (cosL - sinL)/V§” 


M 13 


i=2 




o 



! 1 
o | ; ; 

1 0 







1 «. 0 ! o 

pc 0 , 

1 0 

0 -1 

0 

0 

0 

0 

0 1 

1 0 

1 0 

0 

1 

0 

0 

o 1 

1 1 

0 0 

1 

0 

0 

0 

1 

0 1 

1 0 

-1 0 

0 

1 

0 

0 

0 1 
1 

1 -1 

0 0 

1 

0 

0 

0 

1 

0 I 

1 0 

0 1 

0 
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0 

0 

0 

1 0 

.7 .7 

0 

.5 

0 

0 

O 

1 .7 

0 -.7 

.5 

0 

0 

-5 

0 . 

1 -.7 

-.7 0 

.5 

.5 

.5 

0 

0 i 

1 0 

-.7 .7 

0 

.5 

0 

0 

1 

-.5 

1 -.7 

0 -.7 

.5 

0 

0 

.5 

o 1 

1 .7 

-.7 0 

.5 

.5 

-.5 

0 

0 

1 .7 

0 .7 

.5 

0 

0 

.5 

0 1 


A^P G /At) 42+ 


• 42+ 
M 


1 

0 

0 

-1 

0 

0 

0 

0 

0 

' 0 

T42t- 

“W 

0 

0 

A$(P G /At) 42 ” 


• 42- 
M 


1 

0 

0 

-1 

0 

0 

0 

0 

0 

\ 

1 0 

T42- 

~W 

0 

0 

A$(F G /At) 52+ 


•52+ 

M 


1 

0 

1 

0 

0 

1 

0 

0 

0 

1 T52f 

r 

0 

0 

0 

A$(P G /At) 52 " 


M 52 * 


1 

0 

1 

0 

0 

1 

0 

0 

0 

I T52- 
1 w 

0 

0 

0 

[A<I> P G - 36OOS 0 F 0 )/At] 62+ - w sinL 


M 62 * 


1 

1 

0 

0 

1 

0 

0 

0 

0 

1 

i 0 

0 

1 

0 

[A$P°- 36OOS 0 F^)/At] 62 '-w sinL 


M 62 - 


1 

1 

0 

0 

1 

0 

0 

0 

0 

! « 

0 

0 

1 


3=2 3=2 


B x g 

V 

B Q g 

G l/ 
„ 2 
C SS g 

r- 2 

C IS g 
,, 2 
C IO S 

C 

^OS= 

NOA 


NSA_ 

* X 

6<i ' 

c 

6 <T 

c _ 

3 = 



J = 2 


Notes: (1) w 2 = S 0 P 0 /At (deg/sec) 

(2) c s )r(w T + \v 2 ) m /360 {deg/hr sec) 


_m T+ T- in — v 

(3) c = j,-( w - w ) /360 (deg/hr sec) 

(4) Matrix elements shown as . 7 are actually equal to t Ifl 
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"m 52+ - M 5 
M 52 " - M 5 


*82+ 

M 

.* 82- 
M 




r 

) 




j o j 

1 0 
0 1 

,,♦,42+ *42-. 

(M - M )/2 


« 41 
M 





J 




| 0 c 42 '! 

o 

(M 62+ - M 62 ")/2 


-01 

M 









I" 0 | 

[(M 42+ +M 42_ ) -c 42+ NSA]/2 


M 40 


1 

0 0-1 

0 

0 

0 

0 

0 

'1 

j 


A#(P G /At) 5 - w sinL 


M 5 


1 

10 0 

■ 1 

0 

0 

0 

0 

1 


[(M 62+ +M 62 )-c 62+ NOA]/2 


M 60 


1 

0-10 

0 

1 

0 

0 

0 

! 


A$(P G /At) 4 + w sinL 
e 


M 1 


1 

-10 0 

1 

0 

0 

0 

0 

1 

1 


A${F G /At) 2 + 0 


M 2 


1 

0 10 

0 

1 

0 

0 

0 

1 


A$(P^/At)^ + w cosL 
e 


M 3 


1 

0 0 1 

0 

0 

0 

0 

0 

1 C\ 


A$(P G /At)^ + w^(cosL - sinL)/-\/2~ 

= 

i k 1 

- 

1 

.7 0 .7 

.5 

0 

0 

.5 

0 

1 \J 
1 


A$(P G /At) 8 + 0 


M 8 


1 

0 -.7 -.7 

0 

.5 

0 

0 

.5 

1 


A$(P G /At) 9 + (w e sinL)/y2‘ 


M 9 


1 

-.7 .7 0 

.5 

.5 

-.5 

0 

0 

1 

1 


A$(P G /At) 49 + w^cosL + sinL)//! - 


•10 

M 


1 

-.7 0 .7 

.5 

0 

0 

-.5 

0 

i 

! 


A$(P G /At) 44 + 0 


*ii 

M 


1 

0 .7 -.7 

0 

.5 

0 

0 

-.5 

1 


A$(P G /At) 12 + (w e sinL)//2" 


M 12 


1 

-.7 -.7 0 

.5 

.5 

.5 

0 

0 

1 

1 


A$(P G /At) 13 + 0 


M 13 _ 


1 

0 -.7 .7 

0 

.5 

0 

0 

-.5 

1 

1 



3=3 


3=3- 


A${P G /At) 42+ - w cosL 
e 


-42+ 

M 


1 

0 

0 

-1 

0 

0 

0 

0 

0 

1 

l 

0 

42+ 

w 

0 

0 

A${P G /At) 42 “ W cosL 
e 


♦42- 

M 


1 

0 

0 

-1 

0 

0 

0 

0 

0 


-0 

CO 

1 

0 

0 

[<A$P G -36OOS 0 F^)/At] 52+ - WgSinL 


•52+ 

M 


1 

1 

0 

0 

1 

0 

0 

0 

0 

1 

1 

0 

0 

1 

0 

[(A$P G - 3600S^ F®)/At] 82 - w sinL 

e* 

= 

2 — 
M 

= 

1 

1 

0 

0 

1 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 

A$(P G /At) 62+ 


*•*62+ 

M 


1 

0 

-1 

0 

0 

1 

0 

0 

0 

1 

1 

S2t 

w 

0 

0 

0 

A®(P°/At) 62 " 


• 62- 
M 


1 

0 

-1 

0 

0 

1 

0 

0 

0 

1 

1 

32- 
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0 

0 

0 


3=3 


V 

V 

E 0 g 

C II g2 

c ss g2 

c is g2 

c ic/ 

NOA 

NSA 


6$ 

c 
6ib 

L c J 


3 = 3 


3 = 3 


3=3 


3=3 


Notes: (1) = S^P^/At (deg/sec) 

(2) c m+ a 5r { w T+ + w T ") m /360 (deg/hr sec) 


m T+ t- m 

(3) c = ff(w - w ) /360 (deg/hr sec) 

(4) Matrix elements shown as . 7 are actually equal to l//f~T 
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APPENDIX C 

CALIBRATION ESTIMATION EQUATIONS 


This appendix summarizes the accelerometer and gyro calibration 
estimation equations for each of the four levels of calibration (minimum,, 
partial, full and scale factor error). 


C. 1 EQUATIONS FOR ESTIMATING CALIBRATION TERMS ASSOCIATED 
WITH THE ACCELEROMETERS 


Tables Ol, 02 and C-3 contain the accelerometer calibration 
estimation equations for the minimum, partial and full calibrations, 
respectively. Each table presents the equations for each accelerometer. 
The M"? 1 terms, which are functions of the accelerometer pulse count 


-.A. 


3 




(P ), runtime (At), accelerometer scale factor (A<£> ) and magnitude 
of gravity (g), are computed as indicated in Tables B~2, B _ 3 and B~4. 

The number of calculations can very likely be reduced for manual data 
reduction by expressing certain calibration terms as functions of other 
calibration terms. 


The accelerometer scale factor error estimation equations are 
contained in Table 04, both in general form and for the data obtained 
in Tables 02 and 03. These equations are applied separately to each 
accelerometer (j). 

A 

The calculations for determing the T ~ matrix are listed in 
Section 9. 2 (Eqs. 9-25 through 9~30). 
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c. 2 EQUATIONS FOR ESTIMATING CALIBRATION TERMS 
ASSOCIATED WITH THE GYROS 

Tables 06, 07 and C~8 contain the gyro calibration estimation 

'equations for the minimum, partial and full calibrations, respectively. 

Each table presents the equations for each gyro. The m* 1 terms, which 

G ^ 

are functions of the gyro pulse count (P ), run time (At), gyro scale 

G $ $ 

factor (A$ ), test table pulse count (P ) and scale factor (S ), earth 

rate (w ) and latitude (L), are computed as indicated in Tables B~5, 

B~6 and B _ 7. The number of calculations can very likely be reduced 

for manual data reduction by expressing certain calibration terms as 

functions of other calibration terms. 


The gyro scale factor error estimation equations are contained 

A . A * — 

in Table 09, in which the 60 and 60 terms are determined as shown 

c c 

AgQ 

in Table 05. The calculations to determine the T matrix elements 

A 

are listed in Section 9. 3 (Eqs. 9-33 through 9~3 8), in whic h the NOA and 

A 

NSA terms are computed as shown in Table C~5. It will be noted from 

A A 

Eq. (9-31) that the NOA, NSA pairs are the estimated misalignments of 
the gyro IAs relative to the bubble level coordinate frame. 
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TABLE C-l ACCELEROMETER MINIMUM CALIBRATION ESTIMATION 


EQUATIONS 


Accelerometer #1 


A 








*4 

K 


0 

. 5 

0 

. 5 

0 


M 

A 








• 5 

-(NOA)g 


0 

5 

1 

5 

0 


M 

(NPA)g 

= 

0 

. 5 

0 

5 

0 


M 6 

tfv +1 


0 

0 

0 

0 

1 


M 1 

c 








6V 1 
c 


1 

0 

0 

0 

0 _ 
-I 

- 1 

. q 

_ M _ 


3=1 3=1 


Accelerometer #2 


A 


— 




" 


• 4 

A 1 


.5 

0 

0 

0 

. 5 


M 

U 

A 








*5 

"(NOA)g 


.5 

-1 

0 

0 

. 5 


M 

A 








• 6 

(NPA)g 


5 

0 

0 

0 

5 


M 

A • +1 
5Y 


0 

0 

1 

0 

0 


M 2 

6V 1 
c 

ft 

_0 

0 

0 

1 

0 _ 

n 

_ m 3 _ 


Accelerometer #3 


A 




' 4 

K 

-(NOA)g 

(NPA)g 


0 0 .5 0 .5 

0 0 -. 5 0 .5 

-1 0 .5 0 .5 


M 

M 5 

M 6 

6y +1 

c 


0 10 0 0 


M 1 

Sv—1 

c J 


0 0 0.1 0_ 

* ft 

_M 2 _ 


• The Mj terms are computed as indicated in Tables B-2, B-3 and B-4. 
*Not shown are compensations for rotary axis tilt error (from vertical). 


C-3 



TABLE C-2 ACCELEROMETER PARTIAL CALIBRATION ESTIMATION 

EQUATIONS 


Accelerometer #1 


A 





A i 

A 


0 

. 5 

0 

- (NOA)g 

A 


0 

5 

1 

- ( NPA)g 
A .+1 

= 

0 

. 5 

0 

6 y 

c 


0 

0 

0 

A ,-l 





5 V 

c 

A £ +. 7 


1 

0 

0 

5 Y 
c' 

A * 7 


0 

- 354 

0 

L « v o J 


0 

. 354 

707 


Accelerometer #2 


A 





A 1 

A b 


.5 

0 

0 

- ( NOA)g 

A 


. 5 

-1 

0 

( NPA)g 


. 5 

0 

0 

6v c 
A .-1 

= 

0 

0 

1 

6 V c 


0 

0 

0 

a .+. 7 





6 V c 
a .-.7 


354 

0 

0 

6 V 

L c J 


354 

. 707 

0 




. 5 

0 

0 

0 


M 4 


.5 

0 

0 

0 


M 5 


. 5 

0 

0 

0 


M 6 


0 

1 

0 

0 


M 1 


0 

0 

0 

0 


M 3 


. 354 

0 

1 

0 


M 7 - 

• 5 « 2c io 

. 354 

0 

0 

1 

=> 1 

m 8 + 

3=1 

0 

. 5 

0 

0 


M 4 

1 

0 

. 5 

0 

0 


m d 


0 

5 

0 

0 


CD 


0 

0 

0 

0 


M 2 


1 

0 

0 

0 


m 3 


0 . 

354 

1 

0 


M 8 - 

•5g 2 Cio 

0 

354 

0 

1 


m 9 + 

• 5g C IE 


3 = 2 j=2 
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TABLE 02 ACCELEROMETER PARTIAL CALIBRATION ESTIMATION 

EQUATIONS {Coat. ) 


Accelerometer #3 


1 

<~<f 

1 


1 

0 

. 5 

0 

. 5- 

0 

0 


M 4 

i 

i 

-(NOA)g' 


0 

0 

5 

0 

. 5 

0 

0 


M 5 

i 

<NPA)g 


-1 

0 

. 5 

0 

. 5 

0 

0 


M 6 


A • + 1 

= 

0 

] 

0 

0 

0 

0 

0 


M 1 


o> > 

<1 - 
0 s I 


0 

0 

0 

1 

0 

0 

0 


M 2 


A .+ .7 
6V C 


.7 

0 

354 

0 

354 

1 

0 


M 7 -. 

c 2 

5g C IO 

a .-.7 
__ 6V C 


0 

0 

. 354 

0 

354 

0 

1 


• q o 

M +. 5g 

r 

3 







y 

--3 


CO 

1 

I ^ 


• m 

»The M-j terms are computed as indicated in Tables B-2, B-3 and B-4. 

* Not shown are compensations for rotary axis tilt error (from vertical). 
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TABLE C-3 ACCELEROMETER FULL CALIBRATION ESTIMATION 

EQUATIONS 

Accelerometer #1 


A 

a* 

t> A 


-( NOA)g 

A 

(NPA)g 


0 

. 517 

-. 058 

.459 

0 

0 

0 

0 

0 

0 

-. 435 

+ . 779 

656 

0 

0 

0 

0 

0 

0 

. 459 

. 139 

-. 401 

0 

0 

0 

0 

0 


£ «+l 
6V c 

£ . -1 

6V c 

a . + . 7 
6 V„ 


0 1 


0 0 


0 0 0 


0 0 0 0 0 0 


0.5 0 


0 0 .0 


0 -. 615 1.101 -.928 0 


0 . 5 


0-101 


M 10 


0 -. 649 -.198 . 568 0 +1 0 -1 0 


279j l— 1 J 
3=1 3=1 


(NPA)g 5 


A «+l 

6V c 


*£ 7 

A .-.7 
6V C 

h 

aIP 


. 125 

0 

0 

. 375 

-. 088 

. 088 

088 

. 088 

. 625 

0 

0 

-. 125 

. 265 ■ 

-. 265 

. 265 

265 

. 125 

0 

0 

. 375 

-. 088 

. 088 

1 

o 

CO 

CO 

CO 

CO 

o 

+ 

0 

1 

0 

0 

.0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

-. 088 

0 

0 ■ 

-. 088 

.437 

. 063 

-. 063 

+. 063 

.442 

0 

0 

. 265 

. 187 

. 313 

. 687 

. 313 

-. 354 

0 

0 ■ 

354 

. 250 

. 750 

750 

-1. 250 

-. 354 

0 

0 ■ 

-.354 

-.750 - 

250 

. 250 

-. 250 


|M J ' 

3=2 3=2 
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TABLE C-3 ACCELEROMETER PULL CALIBRATION ESTIMATION 

EQUATIONS (Cont. ) 


Accelerometer #3 


A 






A b 


058 

0 

. 517 

0 

- (NOA)g 


140 

0 

-.459 

0 

(NPA)g 


-.779 

0 

.435 

0 

A , + l 
5V C 


0 

1 

0 

0 

A 

6V C 


0 

0 

0 

1 

A '+ 7 
6 V 
c 

" 

1. 101 

0 

615 

0 

A • -.7 

6 V c 


0 

0 

0 

0 

A 






C IP 


-.197 

0 

-.649 

0 

A 

fio 


-1. 101 

0 

. 615 

0 


3 = 3 


. 459 

0 

0 

0 

. 082 

0 


' 4 
M * 

.401 

0 

0 

0 

. 197 

0 


M 5 

. 656 

0 

0 

0 

-. 313 

0 


M 6 

0 

0 

0 

0 

0 

0 


M 1 

0 

0 

0 

0 

0 

0 


M 2 

.928 

1 

-. 5 

1 

.442 

-. 5 


M 7 

0 

0 

. 5 

0 

■ 0 

. 5 


M 9 

. 568 

0 

-1 

0 

. 279 

1 


M 10 

. 928 

0 

1 

-2 

-. 442 “ 

1 


M 11 






- 


M 12 






j = 

3 

j 


* YQ 

»The M . terms are computed as indicated in Tables B-2, B-3 and B-4. 

J 

* Not shown are compensations for rotary axis tilt error (from vertical). 


C-7 



TABLE C-4 ACCELEROMETER SCALE FACTOR ERROR ESTIMATION 

EQUATIONS 


The following calculations are performed for each accelerometer(j) 


2 , 


^ ' < F 4 Z o - F 2 Z 2 )/(KF 4 - P 2 1 

L = <*(5 Z 1 - F 4 Z 3 )/( V6 - F 4 2 > 
d 2 - (kz 2 - f 2 z 0 )/(kf 4 - f 2 2 ) 

I T7> >7 _ TP ”7 \l IT? 

2 6 4 


y 


D 3 = (F,Z, - t 4 Z 1 )/fF.,F 0 - F. 2 ) 


General form-for + 
pairs of inputs 


where 


and 


z ° ' I <5 ^ , 

Z 1 * z(a i> k 

Z 2 s f(a.) 2 (' 6 vJ 

^ 1 i ; k c k 


K 


F 2 ^ \ 

F 4 ‘ 

K 6 
F c ^ L (a. ) b 

k 


1 1 


_ K 3 a m 

Z 3 _ ^ ( a A (5 F ^ K = number of data samples 


<6V C > » 5 V" 1 

Uv> - 6V- 7 


A * 

(6V cl 


( 6 VJ 


(S V ) 
c 5 




A* -t. 7 

a v 

c 

A t + 1 

6 V 

c 


(a.) = -g 
1 _ 

(a. ) = -g/ /T 
2 

(a. ) = o 
1 3 u 

( a i> 4 = +gA^"2 

(a-) = +g 
1 5 


For the inputs shown, the above general estimation equations 
may be simplified, to yield: 


" 2 ' 

A b = (2. 5 Zq - 3Zg/ g* 4 )/ 3. 5 

A 9 0 

T> 1 - (2. 25 Z 2 - 2. 5Z 3 /gn/. 5g^ I 

D 2 = (-3 Z 0 + 5 Z 2 lg)l 3.5 g 2 
D, = (-2.5 Z x + 3 Z 3 /g 2 )/. 5g 2 

in which only the Z terms need to be computed. 


For the data obtained in 
Tables C-2 and C-3 
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TABLE C-5 GYRO NONORTHOGONALITY AND RATE ERROR 
ESTIMATION EQUATIONS 


„ m- _ , T+ T- X m , 

Define c - ir{ w - w } /360 


Gyro #1 


N($A 


0 

0 

0 

1/c 61 " 

NSA 


0 

0 

1/c 51 " 

0 

A ; + 
A 

= 

1 

0 

0 

0 

_S^c - 


0 

1 

0 

0 


3=1 


r - C.+* 1 

mV 

M 51 
M 61 


3=1 


3=1 


where = a function of table rate (w ) used = 71, 72, 73, 


Gyro #2 


A 

NOA 
NSA 

A •. + 

Qlp 

c 

Jfc J, 


3 = 2 - 

C, 

i 

Gyro #3 


” 

0 

0 

0 

-1/c 51 " 


CNJ 

a 

i 

0 

0 

-1/c 41 " 

0 


• c 2~ 
M 






• 41 

1 

0 

0 

0 


M 






•5i 

0 

1 

0 

0 


M 

- 



- 

3=2 

L. -J 


3=2 


whex*e = a function of table rate (w x ) used = 91, 92, 93, . . . . 


= 


NOA 


A 


NSA 



= 

^ *c 


a , - 


L 8 h , 



0 

0 


0 

0 


-1 Id 
0 

0 


0 1/c 

41- 


6iq 

o 

o 


o 


t~ Cq “I 

M 

<y 

M J 
41 


3 = 3 


3=3 


M 


M 


61 


3=3 


where = a function of table rate (w 1 ) used = 81, 82, 83, 


A • * A A _ — „ 

The units of and M are deg/ hr. , NOA and NSA are in sec, and 
^ * m 

The M^. terms are computed as indicated in 
Tables B-5, B-6, and B-7. 


T 

w in deg/ sec. 
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TABLE C-6 GYRO MINIMUM CALIBRATION ESTIMATION EQUATIONS 


Gyro #1 



• The M . terms are computed as indicated in Tables B-5, B-6 

J 

and B-7. 


• Not shown are compensations for earth rate coupling through 
NOA and NSA. 
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TABLE C-7 GYRO PARTIAL CALIBRATION ESTIMATION EQUATIONS 


Gyro #1 


r A -j 

R t, 


014 

,497 

-. 015 

A 

gBi 


. 268 

.474 

150 

A 


. 118 

. 021 

379 

A 


020 

.496 

021 

Ai 


. 551 

-1 

-. 036 

2 * 





g c ss 


069 

512 

. 429 


3=1 


. 500 

. 003 

. 004 

. 024 

081 

. 026 

-. 786 

. 249 

-. 004 

. 479 

-. 034 

201 

500 

. 004 

. 006 

. 034 

413 

. 006 

.839 

. 060 

-. 498 

. 512 

. 020 

.118 


3*1 


M' 


50 


M 


M 60 

M 1 


M 

M 7 +. 5g 2 C IQ j 
M 8 +.5g 2 C IS 
3=l" 


Gyro #2 


pA — 



R b 


.414 .003 

A 

gBj 


-.150 -.081 

A ' 



gSg 


.021 .499 

gB o 


-. 379 -. 004 

g 2< c n 


-.465 .085 

2* 

8 C SS 


-.429 .498 

j" 

-2 


Gyro # 

"a 

Rj 

3 

.515 -.083 

b 



gBj 


. 026 • 268 

A 


.004 -.020 

gB 0 


-.479 -.118 

8% 

\ 

-. 506 . 620 

9 a 


-.512 .069 

If u ss 


L- 


083 

015 

. 515 

. 142 

268 

. 474 

. 026 

. 249 

020 

-.496 

-. 004 

034 

118 

. 021 

. 479 

-. 201 

620 

494 

506 

058 

069 

. 512 

512 

-. 118 

003 

-. 015 

.414 

. 142 

081 

.474 

-. 150 

. 249 

499 

. 496 

021 

. 034 

004 

-. 021 

. 379 

. 201 

085 

-.494 

-.465 

-. 058 

.498 

. 512 

-.429 

-. 118 


. 024 


M 40 

M 50 

.786 
. 006 


M 6 
’ 1 

. 034 
. 819 


M 

* 3 

. 020 


M 8 +. 5g 2 C 
• q 2 

3=2 


M 5g C Ig 
r ■!=-? 

. 024~j 


• 40 

M 

• 5 

■. 786 


M 

M 60 

. 006 


M 2 

j^ 3 

. 034 


M* 7 “• C io 

. 819 


. 9 . _ 2„ 

020 


M 5g Cjg 

3 

i=3 

j=3 


1=3 


Notes: (1) The Mj 11 terms are computed as indicated in Tables B-5, B-6 


and B-7. 

(2) Not shown are compensations for earth rate coupling through 


NOA and NSA. 
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TABLE C-8 GYRO FULL CALIBRATION ESTIMATION EQUATIONS 


Gyro #1 


-A 

*b ■ 


025 

. 498 

-. Oil'” 

" .498 

Oil 

. 004 

. 022 

. 004 

. 022 

1 

0 

A 

gBi 


. 217 

. 193 

-. 084 

.193 

-. 084 -. 385 

. 168 

-. 385 

. 168 

0 

g^ s 


0 

0 

-. 5 

0 

. 5 

0 

0 

0 

0 

0 

gB 0 


0 

. 5 

0 

-. 5 

0 

0 

0 

0 

0 

0 

2 A 

S C II 

2\ 

g c ss 
2 a 


.460 

-. 716 

054 

-. 716 

-. 054 . 433 

. 107 

. 433 

. 107 

0 


-. 123 

-. 509 

. 448 

-. 509 

.448 . 018- 

. 105 

. 018 

. 105 

0 

g C IS 


0 

0 

. 707 

0 

-.707 

0 

-1 

0 

1 

0 

2 A 

8 C IO 


0 

. 707 

0 

707 

0 

-1 

0 

1 

0 

0 

2 A 

1? C OSL 
i= 

L 

.174 

. 220 

-1. 133 

-1. 194 

. 281 -. 255 - 

. 148 

-. 255 - 

.148 

1 

3= 


M 

M 50 

M 60 


M 

* 9 
M 

M 7 

M 8 

:M 10 

M 11 

;M 12 


3=1 


Gyro #2 


R b 


. 336 

-. 006 

-. 138 

. 037 

.466 

.168 

-. 304 

0 

.071 -. 

030 

. 097 

A 

gBj 

A 


. 325 

-. 054 

. 199 

. 315 

-. 388 

. 062 

-. 261 

0 

-.504 -. 

261 

. 566 

g£! s 

A 


093 

.474 

. 083 

-.348 

. 050 

. 043 

-. 126 

0 

.101 -. 

126 

-. 058 

gB 0 


-. 379 

. 021 

-.001 

-.122 

.479 

-. 100 

. 101 

0 

-.101 . 

101 

. 001 

2 * 

g c n 


-. 827 

. 085 

.547 

-. 493 

-. 217 

. 045 

.409 

0 

.431 . 

409 

387 

0 A 

§ c ss 


-. 300 

. 525 

.079 

. 356 

-. 501 

-.198 

. 119 

0 

-.142 . 

119 

- 056 

2 A 

g Cjg 


. 460 

-. 076 

.281 

.446 

-.548 

. 088 

. 630 

0 

-.713 4 

.370 

.801 

« 2 =io 


. 364 

-. 059 

.004 

. 346 

-.646 

-.717 

-. 287 

0 

. 286 - 

287 

.997 

_ g2 v 


. 028 

-. 128 

-.315 

. 749 

1. 036 

-. 065 

. 379 

-2 

-.287 . 

379 

.223 


3=2 3=2 


M 40 

M 50 

M 6 


M x 

M 3 

M 8 


M" 


M 


M 


10 


M 


11 

12 


M 


13 


3=2 
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TABLE C-8 
Gyro #3 


GYRO FULL CALIBRATION ESTIMATION EQUATIONS 
(cont. ) 


"A 

R b 


.504 

073 

. 048 

. 004 

. 372 . 

124 

-. 051 

.124 

. 031 

051 

-. 031 

A 

gB I 


001 

.441 

. 189 

. 206 

. 048 -. 

046 

-. 395 

-. 046 

-. 012 

-. 395 

. 012 

gB g 


124 

-.026 

358 

. 376 

. 080 . 

044 

-.018 

. 044 

.188 

018 

-. 188 

gB 0 


371 

087 

074 

. 129 

. 239 . 

132 

-. 055 

.132 

-. 144 

-. 055 

. 144 

2A 

§ C II 

r 

-.498 

. 546 

-.181 

-. 205 

-. 570 . 

068 

.386 

. 068 

. 017 

. 386 

-. 017 

2 A 

g c ss 


505 

. 103 

.433 

. 495 

-. 319 -. 

176 

. 073 

-. 176 

-. 044 

. 073 

. 044 

2 A 

g C IS 


175 

036 

-. 506 

. 532 

.113 . 

062 

-1. 026 

. 062 

. 265 

. 974 

-. 265 

2* 

g C IO 


.002 

-.624 

~. 268 

-. 291 

-. 0681. 

066 

. 559 

-.934 

. 016 

. 559 

-. 016 

! aq 

DO 

o> 

O 


. 502 

043 

. 528 

. 502 

.425 . 

073 

-. 030 

. 073 

-. 982 

-. 030 

-1. 013 


M 


40 


M l 


60 


M 


M 

M 3 


M 


M 

M 10 


M 


11 


M 12 

M 13 


Notes: (1) The terms are computed as indicated in Tables 
B-5 3 B-6 & B-7 

(2) Not shown are compensations for earth rate coupling 
through NOA and NS A. 
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TABLE C-9 GYRO SCALE FACTOR ERROR ESTIMATION EQUATIONS 


2>1 


R b *• < Z o - P 2 Z 2> 1 (KP 4 - F 2 > 


= Z 1 ' F 2 


A 

Qr 


KZ 2 / (KF 4 - f 2 ) 


where 

z o * f <«ic> k 


K 


z i 5 ? < w \ z k 


1 

K T 2 

Z 2 = E (W 1 ), z 
1 


'k k 


„ _ K rp 

Z 3 = E (W 1 ). 


1 general form - for + pairs 
’ of inputs. 


TT 

• rp p 

f 2 ■ r (w ) k 

K 

T 4 

F 4 * Z <W >k 


k = number of data samples 
s (6 i ) “ Z /k 

Ci, 0 


A . 


11 • . . i j ■ j ™ 

and z = R ' for W = 0 is included. The (5 1/) , 6f ) pairlTfor various 


W are obtained from Table C-5. 


If a third order coefficient (Q g ) is to be estimated, the F g terms in 
Eqs. (9-7) and (9-9) must be retained, where 


F 


G 


K rp 6 

E (W X ) k 
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APPENDIX D 


EFFECT OF ENVIRONMENTAL ERROR SOURCES 

The environmental error sources are as specified in Figs. 10-1 
and 10-2 of Section 10. 3. In the case of A0^, the angular motion of the 
test stand base, the effect is to cause cross coupling of (1) earth rate 
components to the IA's of the gyros and (2) components of specific force 
due to gravity to the IA's of the accelerometers. Components of the rate 
of change of A0^ (i. e. , A0 ) in the direction of the gyro IAs are sensed 
directly as erroneous rate errors. In the case of a , the linear vibration 
of the test stand, the effect is significant only at the input to the accelerom- 
eters. 


This appendix presents the transfer functions and resulting effects 
associated with the above error sources. The equations are presented in 
this appendix and plots are provided in Appendices E and F. The analyses 
used in the UNIVAC report (Ref. 4) are utilized, and in several instances 
apparent inconsistencies were noted. * Some of these were resolved; others 
are more involved and require more detailed analyses beyond the scope of 
this study. 

D. 1 EFFECT ON ACCELEROMETERS 

The component of SM due to theA0 and a environmental error 

n v 

sources are given in the time domain as follows: 


These are identified wherever they occur. 


D-l 



, tA. 
< S n ) 


m 


= { e 


A 

n 


m 



,m 

(a v ) ]dt 


(D-l) 


where 

, .A .m r g{ M ) for IA level 

(< W = j 

| g(A0 )v2” for I A at 45° to level 
0 for IA vertical 


(D-2) 

D-3) 

(D-4) 


Since A0 and a^ are given in terms of power spectral densities 
(PSDs), the autocorrelation function and Fourier Transform of Eq. (D-l) 
is computed. This was done in Appendix A of Ref. 4 and the results are 
summarized as follows: 

Accelerometer IA vertical 


a V ) = 


At 


m 


S m 


•At m cos 
(At m 2tt1) 2 


2 7rf) 


{ PSD(a ) 
v 


} df 


(D-5) 


Accelerometer IA at 45° to level 



Accelerometer IA level 


a 


V) . 


n 



2(1-At m cos 2?rf 
(At m 27rf) 2 


{ PSD(a y ) + g 2 PSD(A0 n )3 df 

(D-7) 


where the PSDs are given in Figs. 10-2 and 10-1 and approximated as 
shown in Tables 10-2 and 10-1, for the case when the bubble levels are 
not used. 


D-2 



These equations have been solved (in Appendix A of Ref. 4) and the 

2 ,_A 

equations that relate the PSD's to o (e ) are repeated here as follows: 


n 


cr 2 (e^) due tol 

PSD (a ) 

V 


= 1, 1 x 10 
At 


■15 


+ 


1. 92 x 10 
“ (At) 2 


-8 


(D-8) 


2,. A, , 
o - (e ) due to ; 
n 

g 2 PSD(A0 ) 
n 


-1 5 

= 9. 4 x 10" 9 - 2. 3 x 10" 14 At + 6 * 7 X 10 

/ A J \9 


(At ) 2 


(D-9) 


fX 2 (e j4 ) due to 
u n 


g PSD (A0 ) 
2 n 


2 A. 


= 4. 7 x 10 


-9 


^14 -15 

1.15x10 At + 3. 35 x 10 


Gt)‘ 


(D-10) 


where a (e ) has units of g and t is the run time in minutes. The effect 
is shown plotted in Fig. E-6 of Appendix E (as reproduced from 
Table 1-1 of Ref. 8, Fig. 2-3 of Ref. A, and Eq. (D-8), which is derived 
in Appendix A of Ref. 4).* It will be assumed that the curve from Ref. 8 is 
the correct one since the others are based on it and are inconsistent with 
each other. Even if Fig. 2-3 of Ref. 4 were the correct curve, the effect 
on the measurement errors. would be small (see Table E-2). 


Without bubble level compensation, the effect of A0 as given by 
Eqs. (D-9) and (D-10) is to cause excessively large errors (viz. , 97 and 
68 ug for IA horizontal and at 45°, respectively). However, when bubble 
levels are used, as is the case considered in the error analysis, the terms 
in Eqs. (D-9) and (D-10) are reduced considerably. Upon applying the 

*The curve from Ref. 8 is assumed correct, particularly since the two 
.from Ref. A are apparently inconsistent. 
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two bubble level models discussed in Section 10. 3, the resulting average 
a.ccelerometer outpu,i;jsrrors ,are as shown in Fig. E-5 of Appendix E. 


D. 2 EFFECT ON GYROS 


The component of 6 M due to random gyro drift (e ) and the 
. nr 

A 0 and A 0 environmental error sources are given in the time domain 
n n .s. 

as follows: 



) - <e«) 


m 



dt 

(D-ll) 


where 


• G 

e = gyro internal noise (in terms of rate) 
nr 

. q J^e C 0 S L) ^0 for IA vertical 

A0 n. ]w sinL) . , for IA horizontal 
f e A0n 


(D-12) 

(D-13) 


• G lO for IA vertical 

a . . a for IA horizontal 

A0n A0n 


(D-14) 
(D-l 5) 


Q 

Since e and A0^ are given in terms of power spectral densities (PSDs), 
the autocorrelation function and Fourier Transform for Eq. (D-ll) is 
computed. This was done in Appendix A of Ref. 4 and the results are 
summarized as follows: 


Gyro I A vertical 


2 .G m ( 1 

cr (el =\ 7~mr 


n 


A t 


j 2(l-AC .. co| . 2 ?r f) .,t pgD( . e G , 
o (At 2?rf) 


nr 


+ (w cosL) PSD (A0 )} df 
e n 


(D-16) 
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Gyro IA Horizontal 


2 ,. G 


m 


a ( e ) = ( 

n a .m 

A t 


m 

P 20dH co s 2rf ) Cp SD (^ G ) +( W sinL) 2 PSD (A 0 ) 

g (At m 27rf)2 nr e n 


+ PSD( A0 ) df 
n 


(D-17) 


where 


PSD(e G ) = 
nr 


PSD of e 


G 

nr 


PSD{A0 ) 
n 

PSD(A0 ) 
n 


PSD of angular motion of test stand base, 
about level axes 

PSD of angular vibration of test stand, 
about level axes 

(2 7rf) 2 PSD ( A 0 ) 
n 


where the PSDs are given in Tables 4-6 and Fig. 10-1. The standard deviation 
of Eq. (D-17) is plotted in Fig. F-8 of Appendix F. Each of the terms in 
Eqs. (D - 16) and(D“17)have been solved separately (in Appendix A of Ref 4) 

2 - g 

and the equations that relate the PSD's to a (e ) are repeated here as 
follows: 


r 2 , . G 


due tol 17 x 10 -7 + 2 .6.x 1 0~ 7 + 5J5xlO 


PSD (e G ) 
nr 

tr 2 (e G ) due to 
n 


<i> psd <h- 


At 


(At) 


(D-18) 


1. 8 x 10~ 7 - 2. 4 x 10 12 At - 8 ‘ 5 X 9 10 

(At) 2 


-12 


(D-19) 


D-5 



a 2 (e^ ) due H- 

PSD (A0 ) 
n 


4. 35 x 10 
At 


-2 


+ 


2. 11 x 10 
(At ) 2 


-1 


(D _ 20) 


where the variances have units of (deg/hr. ) and At is in minutes. 

Fig. 2-1 of Ref. 4 is a plot of Eqs. (D-18) and (D-19) and indicates 
G 

o(e n ) is equal to 1. 4 mdh for At> 2. 5 min. Although this is inconsistent 
with the above equations, it is used in the error analysis (see Table F-2) 
since the equations in Ref. 4 appear to be inconsistent dimensionally and 
it is more likely that the final results(in the form of curves) are correct. 


Fig. 2-1 of Ref. 4 also contains a plot of Eq(D-20)*, which is 
reproduced in Fig. F-8 of Appendix F. Again the equation and curve are 
inconsistent and the curve is chosen as being correct. Results using the 
equations indicate very large calibration errors which appear to be 
unrealistic. The results -using the curve are one of the largest error 
sources, as shown in Table F-2, and so further analyses are recommended 
to verify results. Fig, F-8 also shows the effect when the gyro IA is at 
45° to the vertical. 

2 G 

The effect one (e ) in Eq. (D-;19) due to using bubble levels, 
which changes PSD (A 0 ), is not significant since the effect was small in 
the first place. 


Although Fig. 2-1 of Ref. 4 includes Eqs. (D-18) and (D-19), the pre- 
dominant term is that in Eq. (D-20). 
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APPENDIX E 


ACCELEROMETER ERROR EQUATIONS 

The form of the propagation of the accelerometer error sources 
into the accelerometer output errors (£>M) are shown in Table E-l for 
accelerometer #1 (j = 1), Each component of 6M is identified by a dash 
number following the letter "A" and is related to the potentially significant 
error sources of Table 10-3 as shown in Tables E-l and E-2. The sensi- 
tivities of the error propagations are indicated by the matrix elements in 
Table E-l. This information is used, in conjunction with numerical values 
for the error sources as listed in Table 10-3, to determine the components 
of 0M for each group of accelerometer orientations (m) and run 
times (At). Those calculations are shown in Table E-2, along with the 
statistical combination* of the SM components to yield estimates of the 
expected errors in M (viz. , 6M) for each test condition (m). The total 
sum-squared error for each test condition (m) is designated e^. Although 
the error propagation equations of Table E-l are only for accelerometer 
#1, the errors can be grouped according to ISU orientation and run time 
(At) and then extended to the other two accelerometers in accordance with 
like conditions. The results are summarized in Table E-2. Not included 
in the Table are pre -calibration alignment errors, which will be considered 
separately as discussed in Section 10. 2 and 10. 3. 


The errors are assumed to be statistically independent and therefore are 
combined on a sum-squared basis. 
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Plots of the effect on 5M of each potentially significant error 
source are shown in Figures E-l thru E-6, These sensitivity curves are 
useful in providing a "balanced design" such that trade-offs can be per- 
formed to ease equipment requirements whenever possible, without com- 
promising calibration precision unnecessarily. Figure E-l graphically 
shows the estimation precision goals stated in Table 4-3 and provides a 
basis for determining acceptable levels of the various error sources. 
However, to use the curves effectively, the propagation of the 6M m 
errors into errors in the estimates of the calibration terms must be 
considered. This sensitivity is given by the calibration estimation equa- 
tions themselves (Equations C-l thru C-3), in which the weighting factors 
are generally less than unity and in some cases will attenuate the errors 
by more than a factor of 10. 

The curves in Figure E-6,as reproduced from Refs. 4 and 8, show 
inconsistencies. The plot of Equation (D-8) is most likely incorrect since 
the equations in Ref. 4 are dimensionally inconsistent and are numerically 
inconsistent with Figure 2-3 of the same report (Ref. 4). Figure 2-3 pro- 
bably reflects the final result with less chance for typographical errors, 
etc. The difference between Figure 2-3 of Ref. 4 and Table 1-1 of Ref. 8 
is not important since the error source is not major (see Table E-2). 

Since Figure 2-3 of Ref. 4 was based on data from Ref. 8, the curve from 
Ref. 8 is assumed to be the correct one. 

The curves presented in Figure E-5 are based on the use of bubble 

levels and reflect the difference between two models for PSD (&0 ) that 

n 

were considered (see Sections 10. 2, 10. 3 and D. 1). Since this is the criti- 
cal error source, it is recommended that further analyses in this area be 
conducted, as discussed in Section 10. 3. 
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Figure E-2 Sensititivity of 6M to Compliance Cross Coupling 
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Figure E~6 Sensitivity of 6M to Linear Vibration (a^) 
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APPENDIX F 


* GYRO 'ERROR EQUATIONS 


The form of the propagation of the gyro error sources into the gjo c 

- * 

output errors (8M) are shown in Table F-l for gyro #l(j = 1). Each com- 
ponent of 5M is identified by a dash number following the letter "G" and is 
related to the potentially significant error sources of Table 10-5 as shown 
in Tables F-l and F-2. The sensitivities of the error propagations are 
indicated by the matrix elements in Table F-l. This information is used, 
in conjunction with numerical values for the error sources as listed in 

Table 10-5 } to determine the components of 6M for each group of gyro 

T 

orientations(m), run times(At) and table rates(W ). Those calculations 
are shown in Table F-2, along with the statistical combination* of the 

• # t 

OM components to yield estimates of the expected errors in M(viz. , 6M) 

for each test condition (m). The total sum-squared error for each test 

condition(m) is designated e m . Although the error propagation equations 

s s 

of Table F-l are only for Gyro #1, the errors can be grouped according 

T 

to ISU orientation, run time(At) and table rate (W ) and then extended to 
the other two gyros in accordance with like conditions. The results are 
summarized in Table F-2. The effect of pre- calibration alignment 
errors is expected to have little effect on the e m in Table F-2. They are 

A gQ. 

important, however, in the estimation precision of T , which is dis- 
cussed i n Section 10. 5. 2. 

* 

Plots of the effect on M of each potentially significant error source 
are shown in Figs. F-l thru F-8. These sensitivity curves are useful in 

n The errors are assumed to be statistically independent and therefore are 
combined on a sum-squared basis. 
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providing a "balanced design" such that trade-offs can be performed to ease 
equipment requirements whenever possible, without compromising calibra- 
tion precision unnecessarily. Figs. F-l thru F-3 graphically show the 
estimation precision goals stated in Table 4-5 and provides a basis for 
determining acceptable levels of the various error sources. However, 
to use the curves effectively, the propagation of the OM errors into 
errors in the estimates of the calibration terms must be considered. 

This sensitivity is given by the calibration estimation equations them- 
selves (Eqs. C-5 thru C-8), in which the weighting factors are generally 
less than unity and in some cases will attenuate the errors by more than 
a factor of 10. 

The "IA Horiz. " curve in Fig. F-8 was reproduced from Fig. 2-1 
of Ref. 4. Although Eq. (D-2 0) is inconsistent with this curve, the curve 

is assumed correct, as explained in Appendix D-2. Because the angular 

* 

rate vibration (A0 n ) is one of the largest error sources , as shown in 
Table F-2, further analyses are recommended. 
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Figure F-4 Sensitivity of 6M to Centrifugal Acceleration 
and Misalignments of IA and SA 



Figure F-5 Seiisitivity of 5 M to Misalignment Between Gyro IA and Rotary Axis 
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FORMATION OF SYMBOLS 


APPENDIX G 
GLOSSARY 


® Wherever possible symbols will be used which suggest the name of 
the parameter involved. 

® Lower case subscripts are used to modify prime symbols and to index 
over several items of the same kind. 

• The lowercase script m isusedto index over different test conditions. 

° Uppercase superscripts will be used to indicate coordinate frames. 

For example, T is used to identify a transformation matrix and 
lettered superscripts such as ML in T MLl identify the particular 
transformation, as being in this case from the L frame to the M 
frame. 

® Matrices will be denoted by capital letters. 

® An underline will identify a vector. 

® Unit vectors are used to identify lines in space such as instrument 
axes and the axes of all frames of reference. 

• The component of any vector along any axis is indicated by a dot 
product of that vector with the unit vector along the axis of interest. 

• A triple line symbol {=) will be used for definitions. 

® A superior denotes an average of the quantity, but is omitted when 
the meaning is clear. 

v 

° A superior " A " denotes an estimate of a vector or scalar. 

9 A superior "dot" denotes scalar differentiation. 
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® The abbreviation "exp" is used to indicate an exponent (e. g. , 

AexpB = A^), but only when necessary to avoid ambiguity. 

6 The Greek sigma (£) will be used for summations. Where the limits 
of summation are clear from the context, they will not be indicated 
with the symbol. 

® The Greek A is generally used to indicate a change. 

® The Greek 6 is always used to indicate an error. 

© S 0 and C0 are sometimes used to identify the sine and cosine of the 
angle 0. 

® Misalignments expressed in the form M a ^ l± c are inf er P re 'ted as the 
misalignment of about with respect to L , measured about 
the positive L axil in accordance with the riglif-hand rule. 
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LIST OF PRIME SYMBOLS 


Old Symbol 


e 

a 

[3] 


Redefined 

[3] 


D 

o 

[3] 


Symbol 


Definition 


a acceleration or specific force, defined 

as sum of applied acceleration and 
gravitational attraction. 

a acceleration due to vibration of test 

stand base 


A 


—1—3—2 


A b’ A b' 

a T V 


unit vector representing accelerometer 
input axis 

unit vectors representing the accelerometer 
system coordinates as defined by the 
accelerometer IAs. 

accelerometer bias, prime refers to . 
quantity modeled in calibration equation 

accelerations along IA, OA, SA and PA 
of gyros and accelerometers, as 
applicable 


B r B s’ B o 


— r— 2’— 3 


B.* A. 
“i -J 

B . G. 
-c -j 


gyro mass unbalance coefficients due 
to acceleration along I A, SA and OA, 
respectively. B = -MUSA and 
Bg= -MUIA 1 

unit vectors representing the 1SU body 
axes.. as defined by the ISU optical 
cube faces corresponding to the 
nominal directions of the gyro IAs 
(assumed to be orthogonal) 

RA 

elements of T 

EC 1 

elements of T 
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Old Symbol 


Symbol 


Definition 


An ’V n G 

[3] 


BL 

c 

C k- C T" C M 

c, 

ir ss 


Bubble Level 

gyro damping coefficient; constant 
(definedwhereused); abbreviation for cosine 

count of pulses in rebalance loop 

major compliance coefficients of gyro 


C IS 3 C IO’ c os 


cross compliance coefficients of gyro 


C IF C XO 


D 


cross compliance coefficients of 
accelerometer (C^.p is the vibro- 
pendulous coefficient and Cjo is 
the compliance of the pendulous axis). 

gyro damping coefficient (= -C/Hr); 
down 



accelerometer dynamic bias 


D r d 2 - d s 


accelerometer scale factor error 
coefficients 


e 


error 


A G 


e , e 


A G 
e , e 
n n 


errors in accelerometer and gyro 
outputs, as defined in Eqs. (7-2) 
and (7-10) 

error in accelerometer and gyro outputs 
due to noise 


♦A 

e , 


•G 

e 


m 

e 

ss 


E 


A G 

average rate of change of e and e 
over time At 

sum squared error in fiM™, due to all 
independent error sources, for test 
condition m 

East 
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Old Symbol 


Redefined 

[3] 


Symbol Definition 

El ] expected value of quantity in brackets 

E unit vector in -East direction 




f 

f<> 


unit vectors representing the Earth's 

coordinate system (E is up, E 9 is East and 

E is North) 
o 

frequency (H^) 

function of quantity in parenthesis; 
numerical subscripts refer to functions 
defined where used 


VV f 3 1 

F 2 ’ F 4’ F 6J 


functions of inputs to inertial sensors 
used in scale factor error estimation 
equations 


F F F 
—1—2-3 


unit vectors representing SATS frame 
of coordinates of base 


g magnitude of specific force due to 

gravitational attraction 

_g the vector directed up that presents 

the specific force due to gravitational 
attraction. Corresponding to popular 
convention, this is referred to as the 
"gravity vector" 

G unit vector representing gyro input 

axis 


G^ G G unit vectors representing the gyro 

system coordinates, as defined by the 
gyro IAs 

H angular momentum of gyro 
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Old Symbol 

Symbol 

Definition 


H 

row vector of coefficients that relate: 
Y to M 


I ’ 

gyro moments of inertia (I, I , I } 


-ihh 

unit vectors representing the TAF 
inner gimbal coordinate system 


IA 

Input Axis 


ISU 

Inertial Sensing Unit 


j 

V^T 


j 

moment of inertia ratio (= I^_ / H } 

OO r 

Redefined 
[3] ' 

K 

constant (defined where used); total 
number of samples used in estimating 
scale factor error coefficients 


i 

radius arm of proof mass in accelerometer 


L 

latitude 


— r— 2 J ~3 

unit vectors representing the Bubble 
Level coordinate system which is 
nominally defined by the TAF inner 
gimbal coordinate system 


m 

moment arm of gyro from test table 
axis; mass 


mdh 

+3 

millidegrees per hour(= deg/hr x 10 ) 

Redefined 

[3] 

-l'-2'-3 

unit vectors representing the ISU 
mirror (optical cube) coordinate 
system (assumed to- be orthogonal) 

Redefined 

[3] 

A 

MJVI 

adjusted measurements of accelerometer 
and gyro outputs 


M • L. 
— c — t 

elements of 
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Old Symbol 


Symbol 


Definition 


MO'3' 


misalignment between O' and O coordinate 
frames, about O' , representing the 
nonortho go nalityTje tween the TAP inner 
and outer gimbal axes 


MR3' 


MU I A, 
MUSA 


misalignment between R 1 and R coordinate 
frames, about R' , representing the 
nonorthogonality'oetween the rotary and 
outer gimbal axes 


Mass unbalance along IA and SA 
of gyro (MUIA = ~ B,-, and 
MUSA = - B ) " 


MOA, MPA, MSA misalignments of IA about OA, SA and 
PA, relative to ISU optical cube 

MIA misalignment of IA(=y MO A 2 + MSA 2 or 

- / v5oa 2 +mpa 2 ) 


Partially 

redefined 


[3] 


NOA, NPA, NSA nonorthogonalities of IA about OA, PA 
and SA, relative to earth coordinates 
for the accelerometers and relative 
to the test table axis for the gyro 

N number of data sampling periods of 

inertial instrument output that are made 
available to the computer; North 


N 


unit vector in the direction of true North 


NUWT nonuniform test table rate 

A G 

O , O unit vectors representing accelerometer 

and gyro output axes 

unit vectors representing the TAP outer 
gimbal coordinate system. The primed 
coordinate frame is defined in Section 5. 


OA 


Output Axis 



Old Symbol 


Symbol 


Definition 


P 

[3] 



PSD( ) 


-A 







E 

q 

[3] 


PA 
A G 

q » q 


Q. 


is 


Q 


IP 


[1] 

Q e ' [3] 


Q r Q 2 

R b- E b 


-1-2-3 

— 1 — 2—3 


■unit vector representing accelerometer 
pendulous axis 


power spectral density, as a function 
of f, of the quantity in parenthesis 


pulse count out of accelerometer, over 
N a 

time At(= E C, ) 
k 


pulse count out of .gyro over time 
N Q 
At (= E Cp 


pulse count out of test table, over 
time At 

pulse count out of time reference 

Pendulous Axis of accelerometer 

quantization in accelerometer and gyro 
outputs to computer 

anisoinertia coefficient relating w^Wg to 
gyro output rate error ( = (I -I )/H ) 

OO li F 

anisoinertia coefficient relating WjW p to 
accelerometer output acceleration error 

gyro scale factor error coefficients 

gyro bias rate error, prime refers 
to quantity modeled in calibration 
equation 

unit vectors representing rotary axis 
coordinate frame that rotates with test 
table. The primed coordinate frame is 
defined in Section 5. 


s abbreviation sine 
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Symbol 


Definition 


Old Symbol 


S 

S 



SA 
SATS 
SE( ) 



T 

t ba 


q g 

[3] 


m ML 

T 


South 

unit vector representing autocollimator 
line of sight 

unit vector representing gyro spin axis 

scale factor of test table 0^ angle 
output pulses 

scale factor of time reference 

Spin Axis of gyro 

Single -Axis Test Stand 

Standard error of quantity in parenthesis 

times at beginning and end of calibration 
run, over which N sensor readouts are 
provided to the computer 

torque; transformation matrix relating 
one coordinate system to another 

transformation matrix for coordinatizing 
vector components in the A coordinate 
frame (accelerometer IAs) to vector 
components in the B frame (body axes) 

transformation matrix for coordinatizing 
vector components in the G coordinate 
frame (gyro IAs) to vector components 
in the B frame (body axes) 


transformation matrix indicating the 
misalignment of the ISU cube (mirror) 
coordinate system relative to the Bubble 
Level coordinate system. 
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Old Symbol 

Symbol 

Definition 


T T T 
- 1—2—3 

unit vectors representing test table 
coordinate frame 


TAF 

Two~Axis Fixture 


U 

unit vector in the up direction 


U 

Up 


V 

velocity 


* 

V 

c 

computer determined acceleration 
input to accelerometer 

co • 

w 

angular velocity 

[3] 

E 

w 

earth rate 

CO 

e 


[3] 

w„ w , w 

rates with respect to inertial space about 


i 0 s 

IA, OA and SA, respectively 

T 

CO 

T 

w 

angular rate of test table 

[3] 


W 

West 


y 

elements of vector Y 


A d 

y 

vector of accelerometer and gyro 
calibration terms to be estimated 

Redefined 

z, , z thru 
k 0 

functions of outputs of inertial sensors, 

[3] 

used in scale factor error estimation 


Z 3 

equations 


5 

error in quantity following symbol 


A 

change in quantity following symbol 
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Old Symbol 

Symbol 

Definition 

D 1 

* A 
A$ 

accelerometer scale factor in velocity units 

[3] 


gyro scale factor in angle units 

Ad 

[3] 

A0 

n 

change in 0 angles (specifically 

0^j and 0-g ) due to motion of test stand 

base 

e 

to 

[3] 

A0 

n 

angular rate of test stand base 
(=?0 2 + 0 2 ) 

i. 

At 

time over which measurement M is 

obtained (= t ~ t ) 

N o 

computer determined change in accelerometer 
input velocity, using accelerometer output 
pulse count and scale factor (-A^ P^-) 

‘n 

[f a.dt] [1] 

t ind 

o 

AV 

c 

< 9 i>lnd 

[1] 

Alb 

c 

computer determined change in gyro 
input angle, using gyro output pulse 
count and scale factor (= A<3>^P^) 


0 

o 

gyro gimbal-to-case angular misalignment 
about OA 


*rV 

*r*4 

angles of Single-Axis Test Table and 
Two -Axis Fixture that define orientation 
of ISU (0i is about trunnion axis, 02 
about rotary axis, 0g about outer gimbal 
axis and 0^ about inner gimbal axis) 


K 

computer determined rate input to gyro 


cr( ) 

standard deviation of quantity in parenthesis 


a 2 ( ) 

variance of quantity in parenthesis 
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Definition 


£y, E6 

[3] 


N 

E summation of pulses out of inertial sensor 

T data sampling period at input to computer 

to angular frequency (rad/sec) 



E 

g 

i 

3 

k 


1 


accelerometer 


computer; index associated with ISU 
optical cube 

earth; index used to identify error 
sources 

East 

gyro 

input axis; 

inertial sensor designator 

index associated with inertial sensor 
pulse rebalance loops; index associated 
with data used to estimate scale factor 
errors 

index associated with bubble level axes 


M time - modulation rebalance loop 

n noise; nominal value 
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Symbol 


Definition 


N 


0 
P 
q 

r 

reb 

s 

T 

u 

v 

V 

z 

List of Superscripts 

A 

B 

E 

F 

' G 

1 
L 


North; reference to number of inertial 
sensor data samplings made available 
to computer 

output axis; initial value (as used in t ) 

pendulous axis 

quantization 

random; record; gyro rotor 

rebalance 

spin axis 

ternary rebalance loop 

uncertainty 

vibration 

vertical 

value resulting from zeroing process 

accelerometer; coordinate frame 
defined- by accelerometer IAs 

body coordinate frame 

earth coordinate frame 

SATS frame of coordinates relative 
to tilt of rotary axis about North 

gyro; coordinate frame defined by 
gyro LAs 

TAF inner gimbal coordinate frame 
Bubble Level coordinate frame 
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Definition 


index associated with test conditions 
for each calibration run 

mirror coordinate frame defined by 
ISU optical cube 

TAF outer gimbal coordinate frame 

rotary axis coordinate frame that 
rotates with test table coordinate 

coordinate frame of autocollimator(s) 

test table; transpose of matrix or 
vector; test table coordinate frame 
(of base) relative to tilt about East 

reference to test table rotary angle output 


plus input rate into gyro 

minus input rate into, gyro 

+1 g input acceleration into accelerometer 

-1 g input acceleration into accelerometer 

+.707 g input acceleration into accelerometer 

-.707 g input acceleration into 
accelerometer 



